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ABSTRACT 


f  PHAS20  performs  simultaneous  multiple  regression  of  a  mathematical 

model  for  the  functional  relations  among  thermodynamic  quantities  to 
experimental  data  for  a  group  of  chemically  related  species.  The 
thermodynamic  quantities  included  are  heat  capacity,  entropy,  enthalpy, 
free  energy,  equilibrium  constant,  electrochemical  potential,  and 
relative  heat  content.  Without  further  adjustment,  P HAS 20  will  fit  the 
model  to  70  sets  of  data  containing  a  total  of  1200  observations  on  20 
species.  The  data  may  include  any  and  all  of  the  following  types  of 
observations : 


For  individual  species: 

1.  heat  capacity 

2 .  third  law  entropy 

3.  enthalpy  of  formation  from  oxides  and 
and  elements 

4.  free  energy  of  formation  from  oxides 
and  elements 


5.  relative  heat  content  data  such  as  are 
obtained  from  by  drop  calorimetry 


For  reactions  : 

1.  heat  capacity  change 

2.  entropy  change 

3 .  enthalpy  change 

4.  enthalpy  change  of  a 
transition  (e.g.,  heat 
of  fusion,  heat;  of 
vaporization,  etc.) 

5.  free  energy  change 

6.  equilibrium  constant 

7.  electrochemical 
potential 


INTRODUCTION 


The  ideal  compilation  of  thermochemical  data,  as  yet  unrealized, 
would  be  both  compete  and  accurate.  The  approach  to  completion  is  a 
function  only  of  the  assiduity  of  the  experimentalist  and  is  therefore 
beyond  the  reach  of  data  compilers.  Assessment  of  the  accuracy  of  a 
thermo chemical  datum,  however,  requires  knowledge  of  the  precision  of 
its  experimental  determination  as  well  as  its  consistency  with  other 
equivalent  measurements  and  with  thermochemical  data  functionally 
related  to  it.  This  program  automatically  fits  to  the  data  a  mathe¬ 
matical  model  consistent  with  che  exact  func^ ' onal  relations  among  heat 
capacity  (Cp) ,  entropy  (S) ,  enthalpy  (H) ,  free  energy  (G) ,  equilibrium 
constant  (K) ,  and  electrochemical  potential  (E)  fov  a  group  of  chemi¬ 
cally  related  species.  The  program  then  returns  a  "best  set"  of 
species-related  constants  from  which  further  thermodynamic  calculations 
may  be  made  or  compilations  assembled.  These  are  accomplished  by  the 
following  steps: 

1.  All  or  selected  data,  weighted  by  the  reciprocal  of  their 
precision,  are  fitted  to  a  mathematical  model  which  is  consistent  with 
the  thermodynamic  theory  for  a  group  of  chemically  related  species. 

2.  Error  plots  for  each  data  set  (grouped  by  type  and  reference 
source)  are  supplied. 


1 


data 


3.  The  user  examines  the  error  plots  and  removes  the  discordant 
sets. 

4..  The  program  is  rerun  and  a  new  fit  of  the  mathematical  model  to 
the  revised  data  is  obtained. 

"5.  Error  plots  for  each  data  set  are  again  supplied  as  a  check 
upon  che  decisions  made  in  step  3. 

6.  Preliminary  table  of  thermodaynamic  values  for  each  species  in 
the  chemical  system  is  printed.  As  a  second  test,  these  values  can  be 
compared  against  existing  compilation  to  insure  all  thermodynamic 
parameters  have  reasonable  values . 

7.  At  the  option  of  the  user,  step-backward  elimination  cf  non¬ 
significant  parameters  in  the  mathematical  model  will  cause  repeating  of 
steps  4,  5,  and  6  for  each  parameter  deleted.  The  maximum  number  of 
steps  to  be  taken  is  set  by  the  user.  However,  if  the  variance  in  all 
the  parameters  drops  below  j.0  percent,  elimination  will  also  stop. 


DESCRIPTION 


The  program  PHAS20  is  divided  into  4  major  parts;  1)  data  input 
and  editing,  2)  multiple  regression,  3)  output  of  error  plots  and 
associated  tubular  data,  and  4)  output  of  tables  of  preliminary 
thermodynamic  values. 


Data  input  and  editing 

The  units  of  observations  found  in  the  literature  are  varied. 
Successful  analysis  requires  a  standard  set  of  units.  In  PHAS20  the 
accepted  units  are  as  follows: 


Parameters  or 
Thermodynamic  Property 
temperature 
pressure 

heat  capacity 

entropy 

enthalpy 

free  energy 

equilibrium  constant 

electrochemical  potential 


Units 

kelvin 

set  by  user  through  the 
data  selection 
joules  mole-1  kelvin-1 
joules  mole-1  kelvin-"*" 
joules  mole-1 
joules  mole 
logarithm  (base  10) 
volt 


Provision  is  made  in  PHAS20  to  accept  most  data  in  the  units  given 
by  the  reference  source.  Th«  temperature  is  edited  by  the  followinq 
equation: 

T  (K)  =  T(input)  +  T(factor) 

where  T( factor)  is  defined  by  the  user.  In  most  instances  T (factor) 
will  have  a  value  273.15  where  temperatures  are  converted  from  the 
centigrade  scale  to  the  kelvin  scale.  Pressure  is  not  an  input  variable. 
To  change  the  units  ft r  thermodynamic  properties  equation  2  is  used: 
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. . . . ■  J-,--, 


Y ( units  given  in  above  listing) 


Y (input)  x  ¥ (factor)  (2) 


where  Y  is  any  thermodynamic  observation  and  Y (factor)  is  the  conversion 
factor  from  the  published  units  to  the  units  used  in  PHAS20.  This 
procedure  simplifies  input  preparation  because  most  data  can  be  punched 
directly  from  the  source  and  edited  on  input  by  PHAS20. 


Multiple  regression 

Multiple  regression  of  the  mathematical  model  described  below  to 
the  data  is  accomplished  by  the  subroutines  0RGLS2  and  MINV20.  These 
have  been  adopted  from  the  Oak  R.idge  general  least  squares  regression 
program  by  Busing  and  Levy  (1962).  The  changes  in  the  above  program 
were  those  necessary  to  adapt  it  to  the  peculiar  problem  encountered 
here.  Since  all  changes  are  superficial,  the  reader  is  referred  to  the 
report  by  Busing  and  Levy  for  the  description  of  the  routine*,. 


The  Mathematical  Model 


Thermodynamic  theory  provides  the  following  exact  functional 
relations  among  the  heat  capacity  (Cp) ,  entropy  (S) ,  enthalpy  (H)  ,  free 
energy  (G) ,  equilibrium  constant  (K) ,  and  electrochemical  potential  (E)  : 


S.  = 
1 


CPi 


■dT 


(3) 


H.  =  'Cp.  dT 
l  *1 


G.  =  H.  -  T  S. 

li  l 


(4) 

(5) 


-R  T  In  K  = 


I 

i=l 


s . 
i 


G. 

l 


(6) 


-n  F  E  =  £  s .  G . 

i-1  1  1 


C71. 


where  T  is  the  thermodynamic  temperature,  R  and  F  are  the  ideal  gas 
constant  and  the  faraday,  respectively,  k  is  the  number  of  species  in 
the  reaction  and  s  is  the  stoichiometric  coefficient,  positive  for 
products  and  negative  for  reactants.  From  these  it  follows  that  an 
empirical  relation  which  adequately  defines  Cp^  as  a  function  of  T  will 
generate  the  mathematical  description  of  the  other  thermodynamic 
parameters.  For  non-aqueous  species,  equation  (8)  has  been  adopted: 


c . 

Cp  =  a  +  2  b.  T  +  — ij  +  f .  t“ 
i  l  l  T  2  i 


+  g./  Si 
i 


(8) 


This  is  an  extension  of  the  Maier-Kelley  equation  (Maier  and  Kelley, 
1932),  the  last  two  terms  being  added  to  describe  adequately  very 


accurate  data. 

For  aqueous  ions  equation  (9)  has  been  adopted; 

c  2  T  r  2  2 

Cp.  =  a.  +  2b.T+  i  +  f.T  +  -g.  —  a  exp  (8  +  otT)  + 
i,aq  i  i  1  i  e  L 

a2  exp  (8  +  aT)  +  —  exp  (8  +  aT)  +  ^ 

where  e  is  the  dielectric  constant  for  liquid  H2O  and  is  given  by 
equation  (10)  : 

e  =  e  exp  .[-exp  (8  +  aT)  -  (10) 

o  o 

The  terms  e  ,  8,  a,  and  0  are  fitted  constants.  Equations  (9)  and  (10) 
are  adopted°from  Helgeson  (1967)  but  retain  the  original  terms  in  the 
Maier-Kelley  equation  (the  first  three  terms  in  the  product  of  equation 
i9)  )  . 

Table  1  gives  the  mathematical  model  as  a  function  of  temperature 
describing  commonly  available  thermodynamic  observations.  For  further 
discussion  the  reader  is  referred  to  Fisher  and  Haas  (1973)  and 
subsequent  reports  presently  in  preparation. 

Output  of  error  plots  and  associated  tabular  data 

After  the  successful  multiple  regression  of  the  mathematical  model 
to  fit  the  data,  the  refined  parameters  are  used  to  calculate  the  follow 
-ing  statistical  data  for  all  observations  in  each  data  set: 

1)  difference,  Y (obs) -Y (calc) 

2)  percent  error,  100  x  (Y (obs ) -Y (calc) )/Y  (obs) 

3)  weighted  difference,  (Y (obs) -Y (calc) ) /standard  derivation  of 
the  experimental  observation 

4)  arithmetic  mean  for  each  of  the  above  type  of  errors 

5)  standard  error  of  each  of  1  thru  3  about  the  respective  means 
4  for  the  data  set. 

As  a  guide  to  the  quality  of  agreement,  '-he  weighted  difference  is 
best.  The  weighted  difference  is  the  number  of  standard  errors  the 
calculated  value  departs  from  the  experimentally  determined  value.  The 
weighted  difference  is  graphed  against  temperature  on  the  error  plots. 

The  PRPLOT  routine  used  in  P HAS 20  was  written  by  Carnahan  and 
Evans  (1961) . 

Tables  of  pieliminary  thermodynami c  values 

The  last  portion  of  P HAS 20  calculates  and  prints  tables  for  each 
species  in  the  group  of  the  chemical  system.  The  tables  contain  the 
following  (the  definition  are  given  by  the  equations) : 
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I 
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TABLE  1.  THE  MATHEMATICAL  MODEL 


A.  Formulae  for  one  nonaqueous  species : 


c.  g. 

Cp.  =  a.  +  2  b.T  +  —  +  f.  T2  +  -  — 

11  1  m2  1 

rp  f- 


S. 

l 


H. 

l 


c ,  T2 

a.  In  T  +  2b.  T  -  —  +  e.  +  f, — 


2TZ 


c . 


2  gi 
v^T 


a .  T  +  b .  T 
i  l 


1  r- 

—  +  d,  +  f,  —  +  2  g.  /T 


(Al) 

(A2) 

(A3) 


c . 

G  =  a.  (T  -  T  In  T)  -  b.T2 - -  +  d  -  e  T 

11  l  _  i  i 

2T 


T'  r- 

f .  —  v  4  g ,  /T 


(A4) 


B-  Formulae  for  one  aqueous  species : 


c . 

Cpi,aq  =  ai  +  2  b.  T  +  -i-  +  f.  T2  -  g .  ±  [  a2  exp2  (0  +  a  T) 


S. 


9  2c t  1  1 

+  a ^  exp  (6  +  a  T)  +  —  exp  (6  +  a  T)  +  — 

6  e2  - 

C  '  m2 

i  aa  =  ai  ln  T  +  2  bi  T  ‘  ~  +  e.  +  f.  — 

'aq  1  2T2  1  1  2 


(Bl) 


H. 

i,aq 


1  i 

-  g.  —  a  exp(6  +  a  T)  +  — 

1  e  L  e  ^ 

c.  f .  I3 

=  a,  T  +  b.  T2  -  —  +  d  +  — — 

i  l  i 

T  3 


+ 


G. 

i  /  aq 


1  -  T 


aexp{$  +  a  T)  -  — 

0  -1 


=  a.  (T  -  T  ln  T)  -  b.  T5 


c . 
i 


2T 


+  u.  e.  T 

i  l 


-  f . 


i  6 


t3  g. 

+  — 


(B2) 


(B3) 


(B4) 
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TABLE  1 .  Continued 


For  two  species  undergoing  an  isochemical  phase  transition  at  T 


(a  -a  )  (T  -T  In  T  ) 
2  1  tr  tr  tr 


(b  -b  )  T 
2  1  tr 


C2"C1 


2  T 


tr 


+  (d2-dl] 


(e  -e  )  T 
2  1  tr 


T  ' 
2  1  tr 


+  4  (92"gl)  ^tr 


tr 
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heat  capacity, 

Cp  (T ) 

entropy , 

S(T) 

enthalpy, 

H  (T ) 

free  energy, 

G  (T) 

equilibrium  constant, 

log  K 

electrochemical  potential, 

E 

relative  heat  content, 

H 

R 

rT  „ 

=  s  +  S£  dT 

0  Jo  T 

-T 
T 

j  298 

-G  (T ) 

2 . 303RT 

-G  (T) 
nF 


=  G  (T)  +  T  S  (T) 

*•-  fT 


f  ,298 


Cp  dT- 


J  298 


-f 

J  298 


Cp  dT 


For  nonaqueous  species  the  data  are  calculated  in  the  tempp“ature 
range  273.15  to  2000K;  for  aqueous  species  the  range  is  273.15  to 
523. 15K.  The  last  portion  of  each  table  contains  the  final  constants 
from  which  the  table  was  calculated. 

Heat  capacity,  entropy,  and  relative  heat  content  are  equivalent 
to  those  given  in  standard  compilations.  The  enthalpy  of  formation  and 
free  energy  of  formation  of  the  species  can  be  calculated  provided  the 
tables  for  the  stable  elemental  species  from  which  the  species  are 
compounded  are  also  available; 

k 

AG°  (species, T)  =  G  (species,  T)  -  L  s,G.  (elemental  species,  T)  (11) 

f  i=l  1  1 

u  k 

AH°  (species,  T)  =  H  (species,  T)  -  E  s,H.  (elemental  species,  T)  (12) 
f  i=l  1  1 

where  k  is  the  number  of  elements  in  the  species  and  s_  is  the  number  of 
atoms  of  the  _ith  element 

INPUT  REQUIREMENTS  AND  DATA  DESCRIPTION 

The  data  deck  must  contain  the  following  information.  Each  item 
will  be  discussed  in  turn  below. 

Run  data 

Problem  data 

Introductory  data 

Observations 

The  Data  sets 

Trial  Parameters 

Plot  control  data  -  -  ..  -  -  - 

IJJIKUSUKM  ralMI'.i;  min  i  l'<17  >;  -  <'T  IlM 
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Run  Data 


The  introductory  data  in  the  data  deck  consists  of  two  cards 
containing  the  information  tabulated  below.  The  first  column  gives 
the  number  of  cards  needed  to  contain  the  data.  The  second  column 
gives  the  variables  on  the  cards.  The  third  column  gives  the  format 
for  each  card.  And  the  last  column  gives  a  brief  description  of  the 
variables. 


No.  cf  Cards 

Variable 

Format 

Description 

1 

DATE 

A8 

Date  or  other  alphanumeric 
identification  of  the  run. 

1 

NREG 

15 

Number  of  separate  problems  in  the 
data  deck. 

Problem  Data 

The  data  for  each  problem  is  logically  broken  down  into  four 
groups:  1)  introductory  data,  2)  observations,  3)  trial  parameters, 
and  4)  plot  control  data. 

Introductory  Data.  The  introductory  data  contains  from  5  to  7 
cards : 


No.  of  cards  Variable  Format 


TITLE  (20)  20A4 


IDO 


NC 


515 


Description 


Alphanumeric  description  of  the 
problem,  80  characters 


Options  :  1  for  both  regression  and 

error  plots, 

2  for  regression  only, 

3  for  error  plots  only. 

Maximum  number  of  cycles  or 
refinements  permitted  on  trial 
parameters.  Usually  set  to  2  or  3. 


"  '  as  K'SVfl  VI  »'|?J  V[  |  \l  oi  f  1  77  t »  4  V;  ,llH 


No.  of  Cards  Variable  Format 

IW 

IL 

IFMOUT 


1  LIS TP  315 

ICY 

NKOLD 

1-3  PNAME (LISTP)  10A8 

1  ION (LISTP)  2012 
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Description 

Weight  control  options ; 

0  for  standard  errors 
supplied  by  users, 

1  for  uniform  standard  error 
of  1  percent  to  be  assigned 
by  the  program  (The  recipro¬ 
cal  of  the  error  is  used  as 
the  weight  on  each  datum) . 
Print  control  options. 

0  for  calculated  results  at 
plot  time  only, 

1  for  tabular  data  for  each 
cycle  of  NC  cycles. 

Format  control  for  punching 
of  refined  parameters. 

Options : 

0  gives  parameters  in  6D12.5 
format , 

1  gives  parameters  in  7A8 
format. 

Note:  With  IFMOUT=l,  there  is 
no  round  off  between  subse¬ 
quent  reruns  of  the  same  or 
related  problems. 

Number  of  chemical  species  in 
the  problem;  limit  is  20 
Maximum  number  of  nonsignifi¬ 
cant  constants  to  be  dropped 
by  step-backward  elimination. 
Numbers  of  parameters  to  be 
forced  during  step-backward 
elimination  even  though 
TEST20  may  indicate  they  are 
non-significant.  Limit  is  40 
parameters . 

8-character  alphanumeric 
identification  of  each  of  the 
species  in  the  chemical 
system. 

Type  of  species: 

1  indicates  a  reference  ele¬ 
mental  species, 

0  indicates  non  reference  and 
nonaqueous  species, 

-1  indicates  aqueous  species 
or  water  (H2O  liquid) . 
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Observations .  PHAS20,  as  described  herein,  will  intercorrelate 


the  type  of  data  described  by  the  equations  in  sections  A,  B,  C,  and  D 
of  Table  1.  Suggestions  for  fitting  the  mathematical  mode]  to  other 
types  of  data  are  described  under  Use  of  the  Subroutine  UNIQUE  below. 

The  data  may  be  grouped  into  a  maximum  of  70  data  sets  containing 
an  aggregate  of  1200  observations.  These  limits  are  set  by  the  dimen¬ 
sion  statements  and  are  readily  changed.  The  grouping  of  data  must  be 
by  identical  type  such  as  all  the  equilibrium  constants  for  the  same 
chemical  reaction  may  be  placed  in  one  data  set. 


But  for  convenience 
in  sorting  out  discordant  data  it  is  recommended  that  such  data  sets  be 
separated  as  to  the  source  reference.  Thus  the  equilibrium  constants 
for  the  same  reaction  but  from  different  sources  should  be  separated 
into  separate  data  sets  and  labeled  completely  on  the  data  set  identi¬ 
fication  card  described  below. 

The  Data  Sets.  The  cards  in  the  observation  deck  are  as  follows 


(Repeat  beginning  after  the  first  card  for  each  data  set) : 


No.  of  Cards 


Variables  & 
Dimension  Control 


Format 


Description 


NSETS 


15 


Number  of  separate  data  sets 
in  the  problem.  Maximum  set 
at  70. 


REF  (10,  NSETS)  10A8 

NPHASE (NSETS) 


615 


80-character  alphanumeric 
description  of  the  data  in 
the  data  set. 

Number  of  species  described 
by  this  set.  Maximum  set  at 


IKOUNT (NSETS) 
IGO (NSETS) 


6.  For  the  heat  capacity  of 
a  single  phase,  set  to  1;  for 
the  reaction  among  4  species, 
set  to  4,  etc. 

Number  of  observations  in  the 
data  set. 

Data  type  options  (see 
Attachment  B  for  examples): 

1  for  Cp  and  ACp, 
for  S  and  AS, 
for  H  and  AH, 
for  G  and  AG, 
for  log  k, 
for  E, 


for  Ht  -Ht 
2  1' 


8  thru  14  optional  for  situa¬ 
tions  described  under  Uses  of 


the  Subroutine  UNIQUE. 
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Format 


of  Cards 


— ~  Variable  &  Format  Description 

Dimension  Control  - - ' - 

ITFACT  Flag  for  editing  Qf  input 

temperatures; 

0  for  editing, 

1  for  no  editing.  See  TFACT 
below. 

IPARA  Flag  for  editing  of  the 

dependent  variable 
YO(I)  and  SIGYO(I) : 

0  for  editing, 

1  for  no  editing. See  PARA 
below. 

ISIG  Flag  for  indicating  the 

type  of  precision  supplied 
for  the  observation: 

0  for  precision  expressed 
as  a  ratio  of  the  observed 
(e.g.  0.01  x  YO) , 

1  for  precision  expressed 
as  an  absolute  range  (e.g., 
0.02  J  mol-lK"l) 

foUoJing0cfar“se:NPHflSE  ^  the  S6t' 


IPARA 


No.  of  Cards 


Variable  & 
Dimension  Control 
ANAME 


Format 


Description 


COEF (NPHASE ,NSETS) 


IS TATE (NPHASE ,NSETS) 


(A8,2X,  8-character  alphanumeric 

DIO . 3  ,  10X,  identifier  of  each  species 
2I5>  referred  to  in  the  data 

set.  ANAME  must  be  found 
in  the  species  list  PNAME, 
character  for  character 
and  space  for  space 
Stoichiometric  coefficient: 

1)  for  species  in  reaction, 
positive  for  products, 
negative  for  reactants;  or 

2)  number  of  moles  upon 
which  observation  was  made 
(generally  1  mole) . 

^  Flag.  If  the  species  is  a 

reference  species  with  AG| 
at  298.15K50  joules  per 
moles,  set  to  1 .  f  ion 
or  H20(1),  set  to  -.1, 
otherwise  set  to  0. 


$4'  'BnWItoNMa  toi 


Mo.  of  Cards 


Variable  & 
Dimension  Control 


Format 


Description 


NINVER (NPHASE ,  NSETS) 


Number  of  inversions 
(equals  the  number  of 
lower  temperature  modifi¬ 
cations  being  considered 
for  the  component) . 
Maximum  is  4. 


If  the  number  of  inversions  is  zero  continue  with  the  next  phase 
in  turn  or  go  to  the  next  section  beginning  after  INVPH  below. 

If  inversions  are  considered,  insert  after  the  phase  identification 
card  the  next  two  cards : 


No.  of  Cards 


Variable  & 
Dimension  Control 
TINV (NPHASfi , 
NSETS,  KINVER) 

INSTAT  (NPHASE, 
NSETS) 


INVSC (NPHASE 
NSETS) 


INVPH (NPHASE, 
NSETS , (KINVER 

+  D) 


Format  Description 

(4D12.5)  Inversion  temperature  at 

which  AG (inversion)  is 
set  at  0  joules 

(15, 13,  If  a  stable  phase  at 

5A8)  298.15  K  is  a  reference 

species,  set  to  1,  other¬ 
wise  set  to  0. 

On  rare  occasions  there 
may  be  a  stoichiometry 
change  at  a  phase  trans¬ 
formation.  INVSC  then 
contains  the  index  of  the 
net  stoichiometry  change 
which  is  found  in  the 
vector  STCOEF (10)  in 
BLOCK  DATA.  The  reaction 
is:  STCOEF  (INVSC)  PHASE1 
=PHASE2 . 

Beginning  with  the  phase 
stable  at  298. 15K,  list 
the  succession  of  phases 
up  temperature  to  and 
including  the  phase  ANAME. 
(All  names  must  be  the 
same  as  found  in  PNAME) 


After  the  phases  and  associated  inversion  data  have  been  identified, 
give  the  data  for  the  set: 


No.  of  Cards  Variable  & 

Dimension  Control 
IKOUNT  X ( 1 , NO) 


Format 


TFACT 


YO  (NO) 
PARA 


Description 

(6D12.5)  First  independent  varia¬ 
ble,  the  temperature  of 
observation. 

Conversion  factor  to 
convert  X(1,N0)  (centi¬ 
grade  scale)  to  X(1,N0) 
(kelvin  scale)  by 
addition 

Dependent  variable 
Conversion  factor  to 
convert  YO(NO)  (any 
units)  to  YO(NO)  in 
J  mol-lx  ,  J  mol-l, 
log  K  (base  10) ,  or 
volts  by  multiplication 
Estimate  of  precision  of 
the  observation  YO 
Second  independent  vari¬ 
able.  For  relative  heat 
content,  X(2,  NO)  is  the 
reference  temperature. 
The  default  is  298. 15K. 
See  Uses  of  the  Subrou¬ 
tine  UNIQUE  for  further 
discussion. 

Repeat,  beginning  with  REF  (10,NSETS)  until  all  data  sets  are  entered. 
This  then  concludes  the  Observation  portion  of  the  input  deck. 

?£lal  Parameters .  The  third  portions  of  the  data  deck  contains 
the  trial  parameters  ; 


SIGYO (NO) 
X  ( 2 ,NO) 


No.  of  Cards  Variable  and 

Dimension  Control 
1  IFMIN  (LISTP) 


Format  Description 

(8011)  Species -related  flag  to  in¬ 
dicate  the  format  used  for 
each  of  LISTP  species: 

0  indicates  (6D12  5/D12.5) 

1  indicates  (7A8) 

Because  there  is  a  format 
flag  for  each  species,  the 
format  for  the  7  trial 
parameters  of  one  species 
may  differ  from  that  of  the 


| 


1 


s,*******im»m*mmu^ 


Plot  Control  Data .  The  plot  control  cards  are  as  follows: 


No.  of  Cards  Variable  and  Format 

Dimension  Control 


NHL 


(415) 


NSBH 


NVL 


NSBV 

1  XMAX  (3E10.3) 

XMIN 

YMAX 


Description 

Number  of  horizontal  lines  in 
the  plot  less  1,  generally 
set  to  4. 

Number  of  spaces  between 
horizontal  lines  plus  1, 
generally  set  to  10. 

Number  of  vertical  lines  in 
the  plot  less  1.  If  the 
temperature  range  is  250K  to 
2000K,  set  to  7.  This  gives 
250  degress  per  division. 
Number  of  spaces  between 
vertical  lines  plus  1.  Gener¬ 
ally  set  to  10. 

REAL*4  variable  giving  the 
upper  bound  on  the  x-coordin- 
ate. 

REAL* 4  variable  giving  the 
lower  bound  on  the  x-coordin- 
ate. 

REAL*4  giving  the  upper  bound 
on  the  y-coordinate.  (YMIN= 
YMAX  for  a  symmetrical 
graph) .  Generally  set  YMAX  to 
10.0EO.  This  allows  weighted 
differences  of  +_  10  standard 
errors  to  be  plotted,  but  re¬ 
tains  enough  resolution  near 
+1  standard  error  to  be 
useful. 


This  concludes  the  deck  for  the  plot  control  cards  and  for  the  problem. 
If  more  than  one  problem  is  contained  in  the  run,  the  TITLE  card  for 
the  next  problem  follows  the  last  card  above. 


Choosing  Trial  Parameters 

Equations  8  and  9  contain  5  parameters  to  describe  the  heat 
capacity  of  the  ith  species.  In  almost  all  situations  these  are  more 
than  adequate,  and  in  most  problems,  using  all  5  parameters  will  cause 
a  singularity  in  the  matrix  inversion.  For  these  reasons  it  is  often 
best  to  graph  the  heat  capacity  data  (or  relative  heat  content  data)  to 
determine  the  parameters  in  the  heat  capacity  equation  (or  the  relative 
heat  content  equation,  Section  C,  Table  1)  which  probably  are  adequate 
to  describe  the  data.  If  no  heat  capacity  or  relative  heat  content 
data  are  available,  then,  at  most,  only  "a",  "a",  and  '‘a"  may  have 
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non-zero  values.  Some  overfitting  can  be  removed  by  step-backward 
elimination. 


Because  successive  refinements  of  the  parameters  is  a  feature  of 
ORGLS  2,  order-of-magnitude  estimation  of  trial  parameters  is  usually 
suffici-ent.  The  following  are  acceptable  first  estimates  for  the  seven 
parameter!. . 


Parameter 

Type  of 

Species 

Non-aqueous 

Aqueous 

a 

10 

10  3 

b 

10“  3 

10 

c 

105 

10  7 

d 

io4  y 

105 

e 

10 

104 

f 

10-7 

10-2 

g 

102 

105 

Constraints 

There  are 

two  preprogrammed  constraints.  The  first 

deals  with 


reference  species  and  associated  reference  state,  and  the  second  with 


phase  transformations. 


Constraint  derived  from  a  reference  state 

The  most  important  constraint  is  that  the  free  energy  of  the 
elemental  species  stable  at  the  reference  temperature  and  reference 
pressure  is  zero.  (Refer  to  Fisher  and  Haas,  1973 

for  a  discussion  of  the  reference  state  con¬ 
vention  adopted  here.)  For  PHAS20  TREF,  the  reference  temperature, 
is  set  at  298.15  in  BLOCK  DATA.  The  reference  pressure  is  never 
explicitly  defined  in  PHAS20  except  for  aqueous  species.  Therefore, 
for  non-aqueous  species ,  the  user  defines  the  reference  pressure  by 
data  selection.  For  aqueous  species  the  algorithm  to  calculate  the 
dielectric  constant  for  water  is  applicable  to  vapor-saturated  water 
only.  Indirectly  then,  the  reference  pressure  is  a  sliding  reference 
set  for  ions  at  the  coexistance  curve  of  water. 

Equation  (El)  of  Table  1  gives  the  algebraic  implications  of  a 
zero  free  energy  at  a  TREF  of  298.15  K.  Equation  (13)  is  derived  from 
(El). 


1/  For  reference  species  or  species  related  to  others  by  phase 
transformation,  "d"  is  set  to  zero  and  not  varied  during  the  regression. 
See  Constraints  below. 
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Mi\  *;■?  rfi is  u?  n<  I  »  o  *:  : 


i 

4 


de  =  ae(298-15-ln  (298.15)  -  298.15) 


+  b  ■ (298.15) 
e 


c 

e 


2*298.15 


+  e  • 298.15  + 
e 


f  (298.15) 3 
e  6 


-  4  g  (298.15)  i/  2 
e 


(13) 


In  the  subroutine  EAFW20,  for  those  phases  specified  as  references 
(i.e.  I STATE  or  INSTAT  =  1) ,  is  calculated  from  TREF  and  the  current 
value  of  tlie  other  parameters  tor  the  reference  species.  Therefore,  in 

selecting  trial  parameters,  SET  De  EQUAL  TO  0.0  AND  SET  KI  CORRESPONDING 
TO  De  TO  0. 


Constraint  derived  from  a  phase  transformation 

At  the  temperature  Ttr  of  a  phase  transformation  the  difference  in 
the  free  energy  between  the  coexiting  phases  is  zero.  This  is  stated 
algebraically  by  equation  (E2)  of  Table  1.  Equation  (E2)  can  be 
rewritten  so  that  d2,  the  "d "-constant  of  the  high  temperature  phase 
may  be  calculated  from  the  remaining  13  constants: 


d  =  (a2-ai )  •  (T  «ln (T  ) 
2  *  x  tr  tr 


(Vci>  „ 

2  T  1 

tr 


T  ) 
tr 


lW 


+  (V6i>  Ttr  +  <VV 


Tr* 

6 


-  4 


(V%> '  V 


1/2 


(14) 


In  the  subroutine  EAFW20 ,  for  those  phases  specified  as  having 
lower  temperature  modifications  (i.e.  NINVER  /  0) ,  d2  is  calculated 
from  the  inversion  temperature  Ttr  and  the  current  value  of  the  other 
six  parameters  for  the  high  temperature  phase  and  the  seven  parameters 
icr  the  lower  temperature  phase.  If  more  than  one  transformation  is 
considered,  then  dp  is  calculated  as  was  d2  using  the  data  for 
successively  lower  temperature  phases.  DIMENSION  statements  limit  the 
number  of  allowed  transformations  to  4. 
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In  selecting  trial  parameters,  SET  "D"  OF  ALL  HIGHER  TEMPERATURE 
POLYMORPHS  EQUAL  TO  0.0  AND  THE  CORRESPONDING  KI  TO  0. 

Use  of  the  Subroutine  PRELIM 

PRELIM  iS  3  feature  retained  from  the  program  of  Busing  and  Levy 
(1962).  In  those  situations  1)  where  the  editing  feature  at  read-in 
is  not  sufficient  to  convert  units  or  prepare  the  data  properly  for 
regression  or  2)  where  a  temporary  adjustment  in  weights  is  desired, 
PRELIM  can  be  used.  .  PRELIM  is  called  from  the  subroutine  0RGLS2  just 
prior  to  the  beginning  of  the  regression  procedure.  Therefore  PRELIM 
presents  a  last  chance  to  make  an  adjustment  in  the  data  used  by  0RGLS2 
or  any  subsequent  routine. 

As  supplied,  PRELIM  is  a  dummy  subroutine  and  returns  control  to 
0RGLS2  without  any  change  in  the  data. 


Use  of  the  Subroutine  UNIQUE 

PHAS20  is  designed  to  allow  simultaneous  regression  of  the 
mathematical  model  given  on  Table  1  to  a  body  of  experimental  observa¬ 
tions  for  each  and  all  species  of  a  chemical  system.  In  some  situations, 
these  are  observations  which  do  not  fall  into  any  of  the  types  giver,  on 
Table  1.  A  regression  result,  to  be  valid,  must  also  permit  accurate 
calculations  of  the  latter  data.  The  Subroutine  UNIQUE  is  designed  to 
allow  the  user  to  fit  the  model  to  these  data  also. 

By  way  of  explanation,  let  us  consider  the  following  problem.  In 
the  System  Sulfur  at  pressures  less  than  one  bar,  the  sum  of  the 
partial  pressures  of  the  gaseous  species  equals  the  total  pressure  over 
the  condensed  phase,  provided  the  vapor  obeys  the  ideal  gas  law: 


8 

p(total)  =.E,  p. 

D=1  3 


(15) 


where  indicates  the  number  of  sulphur  atoms  in  the  eight  polymers 
existing  in  the  sulphur  gas.2-7 


2/ 

-  The  species  3  ,  S  ,  and  S,  „ 

3  1U  1  2 


are  neglected  as  very  minor. 


Equilibrium  between  the  condensed 
be  written 

j  S  (condensed)  • 
and  the  equilibrium  constant  is 


phase  and  the  j_Lh  gaseous 
;  S_.  (ideal  gas) 


species  may 
(16) 


K. 

3 


(17) 


because  the  activity  of  the  condensed  phase  at  vapor  saturation 
Substitution  of  (17)  into  (15)  gives 


is  unity. 
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p  (total)  =  E  K. 


Vapor  pressure-temperature  data  will,  therefore,  provide  addition¬ 
al  constraints  on  the  parameters  aj  thru  gj  which  are  implicit  in  (18)  . 

For  each  datum  the  multiple  regression  routine  0RGLS2  requires 
1)  a  calculated  value  of  the  dependent  variable  obtained  from  the 
appropriate  thermodynamic  model  function  (i.e.,  the  function  given  on 
Table  1),  the  independent  variables,  and  the  current  values  of  the 
adjustable  parameters  aj  thru  gj  and  2)  the  derivatives  of  the 

dependent  variable  with  respect  to  each  of  the  parameters  being  refined 
in  the  regression. 

In  cases  like  sulphur  gas  equilibria  where  the  programmed  code 
(give  algebraically  on  Table  1)  does  not  apply,  the  user  should  take 
the  following  actions : 

^  the  data  are  described  under  Observations  above  except 

IGO  is  given  a  value  of  8  (Also  9  thru  14  are  available  for  other 
constraints  not  included  on  Table  1.  IGO  is  passed  to  UNIQUE  in  common 
statements  and  the  programmer  can  branch  in  UNIQUE  to  the  proper  code 
to  supply  the  requirement  of  2  below)  . 

2)  The  dummy  routine  UNIQUE  is  replaced  by  a  user-coded  routine 
UNIQUE  which  supplies  the  following: 

a)  YC  =  function  (X,  P)  (19) 

dYC 

b)  DCi  =  dP~  =  function  ?)  (20) 

i 

where  YC,  X,  P,  and  DC  are,  respectively,  the  calculated  dependent 
vanabie,  the  independent  variables,  the  parameters,  and  the  derivatives 
o  the  calculated  dependent  variable  with  respect  to  each  of  the 
parameters.  The  subscript  i  above  varies  from  1  thru  NP,  the  number  of 
parameters  in  the  vector  P.  Use  of  UNIQUE  for  vapor  pressure  data  in 
the  System  Sulphur  is  given  in  Attachment  A. 

For  the  mathematical  model  given  on  Table  1,  the  data  cited  by  2a 
and  2b  are  already  programmed  in  EAFW20  and  the  associated  routines 
YDERIV ,  DDERIV,  DIE,  and  FN. 


PROGRAM  RUN  PREPARATION 

The  following  deck  set-up  includes  OS/360  control  cards  as  well  as 
the  problem  deck  requirements.  The  program  is  stored  on  LOADLIB;  the 
card  reader/tape  (SYSIN  file)  ,  the  line  printer/tape  (SYSOUT  file)  and 
punch  (SYSOUT  =  B)  are  the  required  .input/output  devices.  The 
OS/360  control  cards  are  identified  by  //'s  in  columns  1-2.  Each 
control  card  must  be  punched  on  an  IBM  029  keypunch  or  equivalent  card 
codes  must  be  used.  The  symbol  'V  denotes  a  blank  card  column. 

1  User's  JOB  card 

2  / /step„EXEC/vPGM=G466,REGION=250K,TIME=4 

3  //STEP LIB „DD„DSN=SYS1 . LOADLIB ,DISP=SHR 


4  //FT06F001aDDaSYS0UT=A 

5  //FT07F001ADDASYS0UT=B 

6  //FT05F001ADDA  * 

7  data  deck  in  the  following  order: 

run  data  (2  cards) 

First  Problem  title  (1  card) 
option  codes  (1  card) 
number  of  species,  etc.  (1  card) 
species  names  (1  to  3  cards) 

species  types  (1  card) 

number  of  data  sets  (1  card) 

data  sets  (minimum  4  cards  per  set;  1  reference  card,  1 
option  card,  1  to  6  phase  cards  each  followed  by  2 
transition  data  cards  if  applicable,  and  a  minimum  of  1 
observation  card.  The  maximum  number  of  sets  is  70) 
option  card  for  trial  parameter-input  format  (1  card) 
parameters  (1  or  2  cards  per  phase  depending  on 
option  code . ) 

flags  for  parameters  to  be  varied  (1  or  2  cards) 
indices  of  parameters  forced  during  elimination  (0  or  1 
card) 

plot  controls  (2  cards) 

Second  Problem  title  (1  card) 
option  codes  (1  card) 


plot  controls  (2  cards)  for  last  problem. 

8  /* 

9  // 

OUTPUT 

The  printed  output 

The  printed  output  from  PHAS20  can  logically  be  divided  into  four 
parts,  each  of  which  will  be  discussed  below. 

1)  Listing  of  input  data  sets. 

2)  Regression  results. 

)  Printer  plots 

4)  Tables  of  preliminary  thermodynamic  data  and  refined  parameters 
For  example  of  the  printed  output,  refer  to  Attachment  B. 

Listing  of  input  data  sets .  The  first  portion  of  the  printed 
output  contains  an  interpreted  listing  of  the  input  data  sets.  The 
data  for  each  set  has  been  edited  and  should  be  in  units  consistent 
with  the  gas  constant  R  and  the  faraday  F.  If  this  is  not  found  to  be 
true  then  the  results  of  the  regression  will  be  of  no  value.  The 
input  data  sets  are  printed  by  the  main  program  PHAS20 . 

Regression  results .  The  second  part  of  the  output,  the  results  of 
regression  begin  with  a  listing  of  the  regression  options  and  the  trial 


parameters  and  ends  with  a  listing  of  the  correlation  matrix.  The 
regression  results  are  printed  from  statements  in  0RGLS2  except  that 
each  reference  (or  data  set  identification)  is  printed  from  statements 
in  EAFW20.  If  IL  =  1,  the  results  of  all  intermediate  results  will  be 
printed.  If  IL  =  0,  only  the  references,  the  agreement  factors,  and 
the  parameters  will  be  printed  for  each  refinement.  The  tables  are 
self  explanatory. 

Printer  Plots.  After  successful  regression  of  the  model  to  the 
data,  the  program  forms  and  prints  plots  of  the  weighted  difference 
(YO(I) -YC) /SIGYO (I)  as  a  function  of  the  first  independent  variable 
X(1,I) .  For  convenience,  tabular  listings  of  the  data  and  the 
difference  (YO(I)-YC),  the  percent  error  (100.0  * (YO(I) -YC)/YO(I) )  and 
the  weighted  difference  are  printed  just  before  each  error  plot.  The 
mean  of  the  errors  and  the  standard  error  about  the  mean  are  also  given. 
The  error  plots  plus  these  simple  statistics  are  useful  in  sorting  out 
discordant  data  sets . 

The  tabular  data  are  printed  from  commands  in  the  main  program, 

PH AS 20 .  The  reference  for  each  data  set  is  printed  by  EAFW20 .  The 
printer  plots  are  generated  and  printed  by  the  entry  PLOT 4  of  PRPLOT.' 

Preliminary  tables  of  thermochemical  data  and  the  refined 
parameters .  The  last  part  of  each  problem  contains  tabular  listings 
of  Cpi,  S^ ,  Hi,  G^ ,  log  ,  Ei  and  HR  as  calculated  from  equations  Al, 

A2 ,  A3,  A4,  D5,  D6 ,  and  Cl  respectively.  The  terms  log  Ki ,  and  Ei,  are 
respectively  the  equilibrium  constant  for  the  formation  of  the  ith 
phase  at  T  from  the  elements  at  TREF  (298. 15K)  and  the  electrochemical 
potential  of  the  ith  phase  at  T  relative  to  the  elements  at  TREF. 

The  preliminary  tables  are  generated  and  printed  from  statements 
in  PH AS 20. 

The  punched  output 

The  punched  output  is  designed  for  use  of  the  refined  parameters 
in  subsequent  problems  and  in  other  programs  where  thermodynamic  data 
are  needed.  The  punched  deck  contains  the  following: 

1  TITLE  card  for  problem. 

2  Number  of  species  in  problem. 

3  Species  name  card(s) .  The  8-character  alphanumeric  name  of 
each  species  in  the  problem. 

4  The  refined  parameters  P(LISTP*7)  punched  seven  at  a  time  using 
the  format  indicated  by  IFMOUT. 

5  The  complete  refined  parameters  AA(LISTP*7)  punched  seven  at  a 
time  using  the  format  indicated  by  IFMOUT. 

6  Repeat  items  4  and  5  for  each  constant  removed  thru  step- 
backward  elimination. 

7  Repeat  items  1  thru  6  for  each  problem  in  the  computer  run. 

In  the  refined  parameters  P(LISTP*7),  the  fourth  parameter  for 

each  species  will  be  zero  if  the  species  is  a  reference  species 
(ISTATE  =1)  or  a  high  temperature  polymorph  (NINVER  ^  0) .  To  use 
such  a  set  of  parameters,  the  fourth  or  "d"  constant  must  be  calculated 
according  to  equations  (13)  and  (14) ,  respectively.  This  set  of 


parameters  can  be  used  as  input  into  further  refinements  or  into  new 
regressions. 

In  the  complete  refined  parameters  AA(LISTP*7)  ,  the  fourth 
parameter  is  zero  only  thru  elimination.  In  this  set,  as  opposed  to 
the  set  above,  the  fourth  or  »d"  parameter  has  been  calculated  from 
equations  13  and  14  where  applicable.  This  set  of  parameters  can  be 
use  to  calculate  the  thermodynamic  values  on  Table  1,  sections  A  thru 
D  without  consideration  of  the  reference  state  constraint  or  of  the 
parameters  of  other  polymorphs. 

EXPLANATION  OF  ERROR  MESSAGES 

.If  on  lnput/  the  name  (8-character  alphanumeric  label)  of  a  species 
is  mispunched  and  cannot  be  found  in  the  list  of  names  PNAME,  the 
computer  prints  "I,  YOUR  FRIENDLY  COMPUTER,  DO  NOW  QUIT.  IN  YOUR  nth 
DATA  SET  YOU  HAVE  MISTYPED  THE  PHASE  NAME".  A  comparison  of  thujas 
n  the  nth  data  set  with  the  species  name  cards  should  locate  the 
problem.  Execution  of  the  run  is  terminated. 

OF  THO^Tv^TvnS..A-ZER°  °IAG0NAL  ELEMENT  CORRESPONDING  TO  PARAMETER  n 
OF  THOSE  VARIED.  is  generated  when  a  diagonal  element  of  the 

regression  matrix  is  zero.  Find  the  nth  element  by  counting  the  l's  in 
the  list  of  KI  (NP)  at  the  beginning  of  the  regression  results  until  n 
s  reacied.  Because  the  programmed  derivatives  in  EAFW20  will  not 
normally  be  zero,  is  wrong  information  stored  in  KI(NP)?  Is  a 
stoichiometric  coefficient  zero  by  error?  Is  the  user-supplied  UNIQUE 

error"1115  *  deriVatlVe  DC(n)  of  zero?  Th^  job  is  terminated  by  this 

"SINGULARITY  RETURNED  BY  THE  MATRIX  INVERTER"  indicates  the  matrix 

2  !hln9Ular'  T °ftSn  tUS  iS  9P=ss  overfitt^.  chS 

whether  you  are  likely  to  need  all  parameters.  Most  often,  if  "Ci"  and 
fi  are  allowed  to  vary  simultaneously  for  the  ith  species,  a  singular 
matrix  is  generated.  The  parameter  "fi"  should  be  used  only  when  there 
e  mite  curvature  in  the  heat  capacity  plot  at  the  high  temperature 
extreme.  An  example  is  found  on  approaching  the  Curie  pernt  or  lamda 
pernt  with  rising  temperature,  with  this  error  message  S  probes 
terminated  and  the  next  problem  is  run. 

n  'SUBROUTINE  TEST  INDICATES  THAT  JOB  IS  TO  BE  TERMINATED  FOR  REASON 
1.  indicates  that  the  change  in  subsequent  refinements  of  the 
parame  ers  P(l)  is  less  than  10  .  Continued  refinement  is  generally 

sound  statistically.  The  control  is  returned  from  ORGLS 2  to  PHAS20 
for  continuation  of  the  problem. 

WATr  "***WARNING***  THE  TERM  SIG/(NO-NV)  IS  NEGATIVE.  THE  ABSOLUTE 
VALUE  IS  TAKEN  AND  THE  REGRESSION  CONTINUES."  is  self  explanatory 
Before  you  accept  the  results  of  the  regression,  make  sure  they  make 

sense.  The  term  SIG/(NO-NV)  is  used  in  calculating  the  statistics  on 
the  parameters  P(NP). 

"IN  EAFW20  I  IS  GREATER  THAN  NO.  THE  NUMBER  OF  OBSERVATIONS. 
THEREFORE,  I  (THE  COMPUTER)  AM  QUITTING.  THE  PROBLEM  IS  YOURS.'1' 
normally  cannot  be  generated  unless  NO  is  redefined  by  a  user-supplied 
subroutine.  Check  common  blocks  in  user-supplied  PRELIM  or  UNIQUE  for 


NO.  Make  sure  it  was  not  redefined  thru  oversight. 

"YOUR  n-TH  DATA  SET  CALLED  UNIQUE  FROM  EAFW20,  IGO(n)  HAS  A  VALUE 
OF  m  BUT  UNIQUE  IS  CURRENTLY  UNPROGRAMMED.  "  is  most  likely  generated 
when  IGO(n)  has  a  number  greater  than  7  entered  thru  error  at  input. 
Check  you  data  set  option  card.  If  you  have  programmed  UNIQUE,  then 
you  may  still  have  the  dummy  version  in  your  source  deck  or  object 
deck . 

GLOSSARY  OF  MAIN  FORTRAN  SYMBOLS 
NAME  (MINIMUM  DIMENSION)  ,  Definition  or  description 

A (nPHASE;  7)  ,  In  PUTOUT,  local  matrix  containing  the  parameters  for  the 
data  set  being  plotted. 

A,  Elsewhere,  variable  name  assigned  to  the  constant  in  the  equation 
for  the  dielectric  constant  for  water. 

AA(NP) ,  Vector  containing  the  complete  set  of  constants  for  all 
species  in  the  problem. 

ACOEF  (NPHASE ,  NSETS),  REALM  variable  used  for  output  of  the  coeffi-r 
cients  of  a  species. 

AN,  The  number  cf  items  in  a  data  set. 

ANAME,  The  name  (or  identification)  of  a  species,  the  constants  of 
which  ate  constrained  by  the  associated  data  set.  ANAME  must  appear 
in  PNAME  (LISTP ) ,  character  for  character,  blank  for  blank  and  is 
used  to  find  the  index  of  the  species  in  PNAME.  The  index  is  then 
used  in  all  subsequent  calculations  and  is  assigned  to  the  matrix 
IPHASE  (NPHASE,  LISTP). 

ASTAR,  Decoration  for  output. 

AVAL ( 7 ) ,  REALM  variable  used  to  print  preliminary  table  of  thermody¬ 
namic  values . 

B,  Variable  name  assigned  to  the  constant  3  in  the  equation  for  the 
dielectric  constant  for  water. 

BCD,  L0GICAL*1  variable  containing  the  alphanumeric  character  used  on 
the  printer  plot. 

COEF (NPHASE ,  NSETS),  Matrix  of  stoichiometric  coefficients  for  react¬ 
ing  species  in  the  LISTP  data  sets.  The  coefficients  are  negative 
for  reactants  and  positive  for  products  by  convention. 

CPRIME,  When  the  constants  for  the  reaction  are  printed  at  plot  time, 
the  constant  g  is  separated  into  g(non  aqueous)  and  g (aqueous). 

CPRIME  contains  the  g (aqueous)  summation  while  A ( (NPHASE+1) ,7)  con¬ 
tains  the  g(non  aqueous)  summation. 

D,  In  PUTOUT  D  contains  the  result  di  =  ai  (TlnT-T)  +bjT2+Ci/2T+eil4- 
fiT3/6-4g^Tl/2  where  T=298.15,  the  reference  temperature.  The  above 
relation  is  derived  from  the  constraint (E-l)  of  Table  1. 

DATE,  8-character  alphanumeric  label  to  identify  the  run.  The  date  is 
suggested. 

DC(NP),  Vector  containing  the  derivatives  of  YC  with  respect  to  the 
parameter  P (NP) . 

DDDC(7) ,  Component  derivative  of  YC  with  respect  to  the  parameters 
P  (NP)  generated  by  the  constraints  given  on  Table  1,  section  E. 

DDER1V ,  Subroutine  used  to  calculate  DDDC. 
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NAME  (MINIMUM  DIMENSION) ,  Definition  or  description 

DERRP ,  Local  variable  used  in  TEST20. 

DIE,  Functions  which  returns  the  value  for  the  dielectric  constant  of 
water  at  temperature  T  and  vapor  saturation. 

DIEO,  eQ  of  the  equation  for  the  dielectric  constant  of  water. 

DYDC(7),  Component  of  derivative  of  YC  with  respect  to  the  parameters 
P(NP)  at  the  temperature  of  the  observation  (or  reference  temperature 
for  relative  heat  content  data) . 

E(l),  The  weighted  difference  (YO-YC)/SIGYO  passed  to  the  printer  plot 
routine  PRPLOT. 

EAFW20,  Subroutine  which  calculates  YC  and  DC(NP) . 

EBAR(3) ,  Mean  of  the  difference,  the  percent  error,  and  the  weighted 
difference  calculated  at  plot  time  for  each  data  set. 

ERR,  The  difference  YO-YC,  REAL*4  for  use  in  output. 

ERRDP ,  The  difference  YO-YC  used  for  calculations. 

ERRP(NP),  The  error  in  P(NP)  calculated  in  0RGLS2  and  used  in  signifi¬ 
cance  test  in  TEST20. 

F,  The  faraday  (96487.0  C  mol-l) . 

FN,  Function  returning  the  result  FN=exp(8+aT)  from  the  equation  for 
the  dielectric  constant  of  water. 

FNT ,  Result  of  calling  FN (T) . 

FOUR,  The  value  4.0D0. 

I,  Local  index. 

ICY,  Number  of  steps  in  step-backward  elimination.  Default  =  0. 

IDO,  Flag  to  indicate  operations  to  be  executed.  Options:  1  for 
regression  and  plots,  2  for  regression  only,  3  for  plots  only. 

IFMIN ,  Flag  to  indicate  input  format  for  parameters  P.  Options:  0  for 
(6D12.5/D12.5)  ,  1  for  (7A8) . 

IFMOUT ,  Flag  to  indicate  output  format  for  parameters  P.  Options  same 
as  IFMIN. 

IGO,  Flag  to  indicate  type  of  data  in  data  set.  For  the  types  Cp,  S, 
H,  G,  log  K,  E,  and  HT2"HTif  use  I  thru  7,  respectively.  For  types 
programmed  by  the  user  in  UNIQUE,  use  the  unassigned  flags  8  thru  14. 

IGOES ,  Local  variable  having  the  value  IGO+7  and  used  exclusively  in 
call  to  YDERIV  from  EAFW20  or  PHAS20. 

IHOLD (NHOLD)  ,  Index  of  the  parameters  in  P  which  are  to  remain  in  the 
regression  during  step-backward  elimination  without  regard  to 
significance . 

II,  Local  index. 

I  ICY ,  Flag  to  indicate  the  number  of  parameters  dropped  thru  step- 
backward  elimination. 

IINVR,  Local  index. 

IKOUNT(NSETS)  ,  At  input,  the  number  of  observations  in  a  data  set. 
During  run,  IKOUNT  contains  the  index  of  the  last  item  in  each  data 
set. 

IL,  Flag  to  reduce  print.  Options:  IL=0  for  calculated  results  after 
cycle  3  or  plot  time  only,  IL=1  for  calculated  results  after  each 
refinement  cycle. 
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NAME (MINIMUM  DIMENSION) ,  Definition  or  description 

IMAGE (5000),  LOGICAL* 1  vector  containing  the  plot  image  in  vector  form. 

IN,  Local  index. 

INDEX,  Local  index. 

INDEX2,  Local  index. 

INSTAT (NPHASE,  NSETS) ,  Flag  to  indicate  whether  the  species  stable  at 
298. 15K  is  a  reference  species  (=1)  or  not  (=0). 

INVPH (NPHASE,  NSETS,  NINVER) ,  Index  of  the  species  in  PNAME  (LISTP) 
which  are  involved  in  polymorphic  transitions. 

INVSC  (NPHASE,  NSETS) ,  If  the  polymorphic  phase  transformation  at  the 
highest  temperature  involves  a  stoichiometry  change,  find  the 
appropriate  reactant  coefficient  in  the  vector  STCOEF  tnd  assign  the 
index  to  INVSC: 

Index  STCOEF 

1  1.0D0 

2  2.0D0 

3  3.0D0 

4  4.0D0 

5  5.0D0 

6  1/2 . 0D0 

7  1/2. 5D0 

8  1/3. 0D0 

9  1/4 . 0D0 

10  1/5. 0D0 

For  example,  consider  the  reaction:  2  FeCl3 (liquid)  =  Fe2Cl6  (ideal 

gas) .  Here  the  reactant  coefficient  is  2.0  and  the  index  2  should  be 
assigned  to  INVSC. 

ID,  Local  index. 

ION  (LISTP),  Flag  to  indicate  the  type  of  species:  +1  for  reference 
species ,  -1  for  aqueous  species  and  H20(liq) ,  0  for  all  other. 

IPARA,  Flag  to  indicate  whether  input  YO  and  SIGYO  are  to  be  edited 
IPARA=0  executes  YO=YO*PARA,  IPARA=1  executes  no  change  in  YO. 

IPHASE (NPHASE, NSETS) ,  Index  of  the  names  of  the  NPHASE  species  in  the 
vector  PNAME.  Determined  by  the  computer  from  comparison  of  ANAME 
with  PNAME  (LISTP) . 

IREDO,  Index  of  parameter  having  least  significance. 

IREG,  Decal  index. 

ISIG,  Flag  to  indicate  whether  SIGYO  at  input  is  an  absolute  measure 
(-1)  or  relative  measure  (=0)  of  the  precision  of  the  value  of  YO. 

ISING,  Flag  to  indicate  singularity  was  found  by  MINV20,  the  matrix 
inverter. 

ISTATE (NPHASE, NSETS),  If  the  species  is  a  reference  species,  set  to  1 
if  the  species  is  an  ion  or  H20(liq)  set  to  -1,  otherwise  set  to  0 

ISTAY,  Flag  generated  in  TEST20  to  indicate  that  at  least  one  parameter 
has  a  percent  change  from  the  previous  cycle  greater  than  1  x  10“6 

ISTOP,  Flag  generated  in  TEST20  to  indicate  that  all  changes  in 
parameters  from  those  of  the  previous  cycle  are  less  than  1  x  10“® 
percent.  Refinement  will  terminate. 
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NAME (MINIMUM  DIMENSION),  Definition  or  description 

ITFACT,  Flag  to  indicate  that  the  independent  variables  should  be 
edited:  X(J,I)  =  X(J, I)  +  TFACT . 

IW,  Flag:  options  for  weighting;  0  for  user  supplied  weights,  1  for 
weights  uniformly  set  at  1  percent  of  observed. 

IWRITE,  Local  flag  used  in  EAFW20  to  indicate  whether  a  new  data  set 
is  being  calculated.  This  permits  output  of  the  data  set  reference. 

J,  Local  variable. 

JDFLAG ,  Flag  to  indicate  whether  call  to  UNIQUE  is  from  0RGLS2  thru 
EAFW20  ( JDFLAG=0)  or  from  another  routine  (JDFLAG^O) .  If  from  0RGLS2, 
then  the  derivatives  DC(NP)  must  be  supplied  in  addition  to  YC. 

JN,  Local  index.  Just  before  the  call  to  0RGIG2,  JN  is  set  to  zero 
and  becomes  JDFLAG  in  COMMON  or  in  the  arguments  of  the  subroutines. 

JO,  Local  index. 

K,  Local  index. 

KI(NP) ,  Flag  to  indicate  whether  parameters  are  to  refined  (KI (I ) =1) , 
or  held  constant  (KI(I)=0). 

KINVER,  Local  index. 

KO,  Local  index. 

KOUNT,  Local  index. 

L,  Local  index. 

LABEL (NL),  L0GICAL*1  Variable  containing  the  alphanumeric  legend  for 
the  ordinate  of  the  printer  plots. 

LAST,  Local  index 

LISTP ,  The  number  of  species  in  the  vector  PNAME .  Limit  is  set  at  20. 

LL,  Local  index.  '  ” 

LLL ,  Local  i ndex . 

LST,  Local  index. 

LSTINV ,  Local  index. 

LSTPHA,  Local  index. 

M,  Local  index. 

NC,  The  number  of  refinement  cycles  to  be  taken  by  0RGLS2.  Generally 
set  to  2  or  3. 

NDATA,  The  number  of  items  in  a  data  set. 

NHL,  The  number  of  horizontal  lines  on  the  printer  plot  less  1. 

NHOLD,  The  number  of  parameters  to  be  forced  to  remain  in  the  regres¬ 
sion  without  regard  to  the  results  of  TEST20 . 

NINVER (NPHASE ,NSETS) ,  The  number  of  phase  transitions  which  occur  in 
the  chemical  component  between  the  temperature  of  the  observation 
X(1,I)  and  the  reference  temperature  for  the  elements  at  298. 15K. 

NL,  The  number  of  characters  in  the  legend  LABEL. 

NO,  The  total  number  of  observations  in  the  problem.  Limit  is  1200 

observations  total.  ' 

NOYES,  Local  index. 

NP,  The  numbers  of  parameters.  The  limit  is  set  at  140. 

NPHAS ,  Local  index. 

NPHASE (NSETS) ,  The  number  of  species  in  the  jth  set  of  data. 

NREG,  The  number  of  problems  in  the  computer  run.  No  limit. 

2  7 
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*escripti  on 

Printer  plot.  P“CS"  betwe6n  the  horizontal  lines  plus  i  on  the 

"er'pl^r  "  —  b—  vertical  lines  plus  1  on  tte  print- 

NSCALE  (5)  ,  Arguments  used  to  set  ud  i-ho  • 

Here  all  values  are  set  to  zero  P  lma9e  by  the  routine  PRPLOT. 

"T;  p^rerS  °f  data  “*  “  a  U-it  is  70  data  sets 

S:  The  "Se“  °J  “Xbr  hTied  (i-e-  “«>-U  • 

ONE,  The  constant  1.0D0.  tables.  Limit  set  at  2 - 

P  (LISTP*7^  f°p  °"e0pe^^0ne^s°ofGpHA^LBUSing  LeVy  Q962) 

Pa£  in™*  Parameters  the  species. 

™(NV),  ?a^^a~  ^2  f6  dependent  variable  Y0. 

changed.  ges  for  those  parameters  which  are 

'"Jnvirs'i^^enL^'^iTwrTSeT'e  ^  ^  PbaSe  ib 

to  include  the  phase  stLle  at  *he  phaSe  stable  298. 15K  and  progress 

PLOTg '  E'™  POintS  into  0,6  -ti  Sot”  °f  obseb''ati°h- 

PLOT3 

PLOT4 

"“bl-”' '  8-charact“  alphanumeric  labels  for  each  species  in  the 
PUTOUT,  Routine  called  at  Dint  F-i 

constants,  including  "dj  "  if  constra  i  °  calculate  1:116  complete  set  of 
R/  The  gas  constant  (8  31  d?  7  ?tS  are  applicable. 

REF(10,NSETS)  ,  80  character  T^h  *  }  timeS  ln<10)  =  19-14437900 

HELHRH,  Local  varlSS^^Sr^  ^  **  eaCh  data  3et 
value  relative  to  "the  oblervefdatr0611"  °f  eaCh  calc^ated 

si«  s?  -  r— — 

SERRSQ  (3)  ,  The  ITumTof  ^rrorfsq* 

SGN (2) ,  Local  variable  used  d  ■  £  3  data  set- 

derived  from  the  ^hL  p^eSrl  CalcUlation  of  'V  which  is 
SIGYO(NO)  ,  The  standard  error  of  an  T^6  °f  COnstraints • 

is  used  as  the  weight.  observation  Y0(I) .  The  reciprocal 

SIX,  The  constant  6.0D0. 

STCOEF  (19)  ,  Vector  of  possible  stni  v,  • 

Change  in  stoichiometry  at  '  t“c  “KtIK  coefficients  for  a 
above.  7  at  a  pbase  transformation.  Refer  to  ihvsc 

STDEV  (3)  ,  The  standard  error  n-F  Fh 

a  data  set.  6  ata  ab°ut  the  mean  error  EBAR  for 

T,  Local  variable. 
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NAME (MINIMUM  DIMENSION),  Definition  or  description 

TEST20 ,  Routine  which  does  the  following:  1)  Signals  the  end  of 

refinement  cycling  if  the  percent  change  in  a  parameter  is  less  than 
10  2)  Locates  the  least  significant  parameter  being  refined  but 

not  constrained  by  IHOLD  to  stay  in  step-backward  elimination, 

3)  Signals  the  end  of  step-backward  elimination  if  the  error  on  all 
parameters  is  less  than  10  percent. 

TFACT ,  Local  input  variable  used  in  editing  X(J,I)  by  addition. 

THETA,  The  constant  0  of  the  equation  for  the  dielectric  constant  of 
water . 

THREE,  The  constant  3.0D0. 

TINV (NPHASE ,NSETS ,NINVER)  ,  The  inversion  temperatures  beginning  with 
the  closest  temperature  to  *.98.15K. 

TITLE (20)  ,  80  character  alphanumeric  label  for  the  problem. 

TK(  19 ,2)  ,  Array  of  temperatures  used  to  calculate  the  preliminary 
tables  of  thermodynamic  constants  for  a  species. 

TO,  Local  variable. 

TREF,  The  constant  298.15  which  is  the  reference  temperature  from  the 
reference  state  convention. 

TWO,  The  constant  2.0D0. 

TYPE ,  Legends  for  output  of  data  type . 

UNIQUE,  Subroutine  to  be  programmed  by  usex*  if  data  not  considered  on 
Table  1  are  available  and  used  in  refinement.  Refer  to  Use  of  the 
Subroutine  UNIQUE.  ~~ 

UNITS,  8  character  alphanumeric  label  for  the  type  of  energy  data 
consistent  with  R  and  F.  As  supplied,  UNITS  has  the  label  'JOULES'. 

VAL ( 7) ,  Local  variables. 

X (2 ,N0) ,  The  independent  variable.  X(1,N0)  is  the  absolute  tempera¬ 
ture.  X(2,NO)  is  an  optional  second  independent  variable  used  in 
the  current  code  for  the  reference  temperature  of  relative  heat 
content  data. 

XI  (1)  ,  REAL*4  variable  containing  the  x-coordinate  of  the  datum  being 
passed  to  the  routine  PRPLOT. 

XMIN,  The  lowest  temperature  for  the  error  plot. 

XMAX,  The  highest  temperature  for  the  error  plot 

YC,  The  value  of  the  dependent  variable  calculated  from  the  current 
set  of  parameters. 

YDERIV,  Routine  which  calculates  the  temperature-dependent  part  of  the 
derivative  d  YC/d  P^ 

YESN0(2)  ,  The  labels  'YES'  and  'NO'  used  in  printing  input  data. 

YMAX,  The  upper  bound  of  the  y-axis  on  the  error  plot.  YMAX=  -YMIN . 

YMIN,  The  lower  bound  of  the  y-axis  on  the  error  plot.  YMIN=  -YMAX. 

YO(NO),  The  dependent  variable. 

ZERO,  The  constant  0.0D0. 
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ATTACHMENT  A 

Example  of  the  Use  of  the  Subroutine  UNIQUE 

In  the  text,  there  is  a  discussion  of  the  use  of  the  Subroutine 
UNIQUE  with  the  System  Sulfur  given  as  an  example.  The  attached 
routines  UNIQUE,  VPRESS,  and  KJCALC  will  provide  values  for  YC  and 
DC(NP)  for  an  IGO  of  8. 

The  subroutine  UNIQUE  calls  the  appropriate  subroutine  VPRESS  when 
IGO  for  the  data  set  is  8.  If  IGO  of  the  data  set  is  9 ,  10,  11  or  12, 
subroutines  SNCOEF,  ACTSUL,  INFLO,  and  ININT  are  called.  These  are  not 
given,  but  would  contain  other  constraints  on  the  parameters  for  the 
gaseous  species  in  the  System  Sulfur. 

The  subroutine  KJCALC  is  a  modified  version  of  EAFW20  and  returns 
the  eight  constants  KJ  for  reaction  16  of  the  text.  KJCALC  also  gen¬ 
erates  the  data  for  the  matrix  TDC(8,NP) ,  the  derivatives  of  aach 
equilibrium  constant  with  respect  to  the  NP  parameters  of  the  vector  P. 

The  subroutine  VPRESS  calculates  YC  and  DN (NP)  by  summing  KJ(8) 
and  TDC(8,NP)  for  the  eight  reactions. 


SUeROUTINF  UNIQUE  I YC.N , JO,  JDFL AG  I 
IMPLICIT  PEAL*8IA-H,C-Z! 

DI-LNSION  CQE*16,70I  ,PNAREt20)  , TI NV C 6 , 70 , 4 1  ,  I FHASF 1 6 , 70 1 , 

1  NPHASt(70I.IK3UMC  70),  I  GUI  70)  ,  ESTATE!  6,70)  ,MNVFR(6,70)  , 

2  H  STATE  t  ,7C»,INVFHI6,  70,01  ,II,VSCC6, 701 

C0K1CN  /tAKTH/  CCEF  ,  FNAPF  ,  T  INV  ,  I  PHA  SE  *NPMSE  »  I KCUNT  ,  I GC  ,NSFTS« 
1  IS  TAT  c  ,  N  INVER  ,  INST  AT  ,  INVPH,INVSC,  LISTP 
M=  IGC ( JO  I -7 


GO  rr  l  IOC, 200,  30C,  ACC, 500)  , 
100  CONTINUE 

CALL  VPRfcJS  ETC , N , JC , JOFLAG ) 
RETURN 


200  CONTINUE 

CALL  SNCOEF  (»C  ,N  ,  JC  ,  JOFLAG  I 


RETURN 
300  CONTINUE 

CALL  ACTSLL  t»C,N . JO, JCFLAGI 


RETURN 
400  CJNTINUE 

CALL  INFLC  I  YC  ,N  ,  JO  ,  JCFLAGI 
RETURN 

500  CONTINUE 

CALL  ININT  IYC.N, JO, JCFLAGI 

RETURN 
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SUfiKOUT INF  V PRESS  < Y C .  I , JO, JCFL AG  I  VPR  10 

IMPLICIT  FEAL*8<A-H,C-Z»  VPR  20 

REAL**  T I TLE ( 20  I  VPR  30 

RFAL»8  LOOK<8I,KJ<8)  VPP  AO 

DIMENSION  TCCIS.IAOI  VPR  50 

DIMENSION  XI2, 120  0  ,VCt  12001,  S1GYCI  12001  ,PI  IACI  ,KI(HOI ,  VPR  60 

1  OCUAJ)  .PCUAOI  VPR  70 

CCMMCN  /AIR/  X,  P,  OC ,  TITLE,  YC ,  SIGYC,  FC,  KI,  NC,  NV  ,  NX,  IU,  VPR  80 

1  NP,  NO,  I S I N  G ,  1STCP,  1L  V®R  90 

CC*MCN  /NATGR/ZERC,  CN  =  ,  TwC,  THREE,  FOUR,  SIX,  R,  F  VPR  100 

COMMON  /6PI M/LOGK ,KJ , TDC  VPR  110 

YC*ZERO  VPR  120 

T* X(  I ,  I  I  VPP  I  30 

CALL  KJCALC  IT.JOI  VPP  1A3 

CO  1  0 J  J  =  2,8  VPR  15C 

YC*YC*KJ( Jl  VPP  160 

100  CCNTINJL  VPR  170 

IF  ( JCF L AC « Nc « 0 j  RETURN  VPR  180 

0 J  2 CO  K»I,NP  VPP190 

CO  200  J-2,8  VPR  200 

CC(KIOC(KI»TOC(  J.KI  VPR  210 

200  CONTINUE  VPR  220 

RETURN  VPR  230 

SNO  VPR  2A0- 
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scbrcutine  kjcalc  <t.j)  kjc  10 

INFLlCtT  PSaL»8<A-F,  C-Z )  KJC  20 

C ( “i N  S ( CN  SGM2)  ,  CC  (  1AO)  KJC  3Q 

01  PE  NSC. IN  C3EF(6«7CI,PNAPE<20I,TINV(6,73,A1,1PHA<E(6,70),  KJC  A3 

1  NPHASL<701tIKCUM(70»,ir,Q(7JI,ISTATH6,70»,Nl,VER<6,70),  KJC  50 

2  1NSTATI6 ,70) ,INVFhC6, 70,51, 1NV$C(6,70)  KJC  60 

ctAL*A  riTLE<20)  *  KJC  70 

CIPLNSIUN  XI2.120CJ ,YC( 120J|,$1GY0( (2001  ,F(  1ACI.K1 (  1AC),  KJC  80 

1  ll ( 1AJ)  ,P3< LAO)  KJC  90 

Cl“i.NSICN  P5F<10,73),FRPP<1A0I  KJC  130 

C1“LNS1CN  CLCCI7I ,CYCC<71  KJC  110 

OCPJNSIJN  SCINV<2I,STCCEF( 101  KJC  120 

9fAL*3  LCCK(3I,KJ(81  KJC  130 

CIK£H1  CN  TDC  (  8  «  1  A3  1  KJC  liO 

CC“NU  /EARTH/  CCf  F  ,FNA“E,T1NV,  1PHA  SE.NPFASE  ,  1KCLNT,  IGC.NSFTS,  KJC  150 

1  I  STATE,  MWE7,  (NSTAT,  (NVPF,  1NVSC.L1STP  KJC  160 

CC“NCN  /  A  T  F  /  x,  P,  ll,  TITLE,  YC,  S1GY0,  FC,  K  1 ,  NC,  MV,  NX,  1 V, ,  KJC  173 
1  NP,  Nl,  1 S ! NG ,  lST'-P,  1L,  JCPLAG  KJC  180 

CO“MCN  /FI  PI/  RcF,  EPFP,  I  »p IT" »  ICY,  llCY,  KCCC  KJC  190 

C0“MC1  /1,/TFP/itHC,  CAS ,  T  NC ,  ThPFL,  Ft  L  F  ,  SIX,  P,  F,  KJC  200 

1  SCINV,  TP  "  F ,  STCCFF  KJC  210 

CCNNCN  /SPAC*/  CCCC,  OYOC,  SC,  TC  KJC  220 

CC““CN  /JP1M/LQGK,KJ,TCC  KJC  230 

CC““CN  /„C“AN/LMC,LCCC  KJC  2A0 

CALL  Y  0  r.  A  I V  (T,5I  KJC  250 

IF  < (G3< J I.NC.12I  GC  TC  10J  KJC  260 

*1=1 F1ASS (3, J)  KJC  270 

*2=1  phase i  3 ,  j  j  kjc  2ao 

CC  TC  230  KJC  290 

100  CCNT1NUS  KJC  300 

“1  =  2  KJC  310 

“2=8  KJC  320 

230  CONTINUE  KJC  330 

CO  1 8 J J  M=Ml,H2  KJC  3A0 

CO  3C0  K=l,NP  KJC  350 

CC(K|  =  7ERC  KJC  360 

TCCI “,<l =7ER0  KJC  370 

300  CONTINUE  KJC  380 

LOGK <“l=ZERO  KJC  390 

NTSMF  =  :iPHASE<  Jl  KJC  AOO 

NPFASc(J) =2  KJC  A10 

(PHASEI2 , Jl=“  KJC  A20 

C0FF(1,JI«-F  KJC  A30 

C0EFI2, Jl =1  KJC  AAO 

ISTATc(2, J1=0  KJC  A50 

NINVERI2, J1=0  KJC  A60 

LAST<NPHA<E< J>  KJC  A70 

CO  1501  L  =  l ,  LAST  "JC  A80 

lNCtX  =  l*7«<IPHASE<L,Jl-n  KJC  A90 

SC  =  CCEF(L  ,JI  KJC  500 

IF  <  lSTATEd, Jl.NE.lt  GO  TO  500  KJC  510 

CALL  03t  R  IV  (TREF)  KJC  520 

CO  AOO  K *  1 , 7  KJC  530 

CC<K*INDEX-ll»CC<KMNCcX-ll»SC*OCOC<KI*OYCC<AI  KJC  5A0 

AOO  CONTINUE  KJC  550 

GC  TC  130C  KJC  560 

500  CONTINUE  KJC  570 

IF  <MNVEML,J).LE.O)  GC  TO  1300  KJC  580 

CO  6 CO  L L L= 1 , 2  KJC  590 

SGN(LLU  =  SCINV<LLL)  KJC  600 
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6jo  continue  kjc  6io 

IF  (  INV  SC  (  L  ,  J  I  .  L:C  . 0  I  GC  TO  700  KJC  620 

SGM  1  l  =  SGN<  I  >  *  S  T  CC  =  r  (  INVSCI  L,  Jl  I  KJC  630 

SGN(i)=SGN(2l«STCCFF(INVSC(L,JII  KJC  660 

700  CONTINUE  KJC  650 

IF  (  INSTAT(L,J1.EC.0I  GC  TO  900  KJC  660 

IINV»*l*7»(INVPH(L,J, 11-11  KJC  670 

CALL  OJcHIV  (TREF)  KJC  68C 

CO  30J  K=l,7  KJC  690 

rC  IK  *  IIN'/o-l  )  =  CC(  MI  IN'/R-ll  »ODOC<K  HSCNI  2  1*  SC  »0  YCC  (  6 )  KJC  700 

8J0  CONTINUE  KJC  710 

900  CONTINUE  KJC  720 

LSTINV=NINVER(L,J1  KJC  730 

DO  1200  LL=1,LSTINV  KJC  760 

CALL  COw^IV  (TINV(LtJiLL)l  KJC  750 

IF  (  INVSCIL,  Jl  .CC.OI  GC  TO  1000  KJC  760 

IF  (LL.N'.LSTINVI  GO  TC  UoO  KJC  770 

SGM  2I  =  SGM2I/STCC-F(ISVSC(L,JI  I  KJC  700 

1J00  CCMIVJ"  KJC  750 

DO  11CJ  L  LL  =  1  *  2  KJC  800 

IISVP  =  U7»(INVPH(L»J,(LLFLLL-1)1-11  KJC  810 

00  110J  K»l,7  KJC  820 

OC  (K  +  IINVP— 1)  =  0C(  MI  INVP-ll  «000C  (K  I  ♦  SGNI  L  LL  1  *  SC  «CY  OC  I  6  )  KJC  830 

1100  CONTINUE  KJC  860 

12J0  CCNT  INUE  KJC  85C 

1300  COM  T  INUE  KJC  860 

CO  1 6  J J  K»l,7  KJC  870 

nC(K«IN^O-ll=DC(K«IKCFX-lMSC*CYCC(K>  KJC  880 

1 6 j o  continue  kjc  ooo 

1500  CCNT -NUE  KJC  900 

CO  160 J  K=1,NP  KJC  910 

L0GMM)*LC0K(M|4P<KMCC(K1  kjc  920 

1600  CONTINUE  KJC  930 

IF  (lOGKINJ.Gr. 7. sesscil  GO  TO  1900  KJC  960 

KJ(rt|-1.0Cl**LCGK(HI  KJC  950 

OC  1700  K  =  1  *  NP  KJC  960 

TOCIF.K  )«TOC(M,K  I«KJ(H*IN10*0C(K)  KJC  970 

1700  CONTINUE  KJC  980 

1800  CONTINUE  KJC  990 

NPFASEU1-NTEMP  KJC1000 

COEF  (  1 .  J  I  *CN6  KJC1010 

RETURN  K  JC  1020 

1900  CONTINUE  KJC  1030 

FRITE  (6,20001  N , LOCK (N )  KJC1060 

STCP  K  JC 1050 

2000  FCPMAT  ( • CFOR  YOU  F  •  ,  1 2,  '-TH  SPECIES,  YGU  FAVE  TPIFC  TO  CALCULATE  1KJC1060 
lO.O**',oi2.5,,.,/<  I  JLST  CANNOT  HACK  SUCH  LAPGE  NUPBE0  5 •  IT  IS  TKJC1070 
2IRIN0  ON  TEE  DRAIN. »l  KJC1080 

RNO  KJC 1090— 
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ATTACHMENT  B 

EXAMPLES  OF  AN  EXECUTION  OF  PHAS20 
Problem  1.  The  System  H2O  and  the  self-ionization  of  water.  . 
Probelm  2.  The  System  Ni-0  and  the  thermodynamics  of 

Bunsenite  (NiO) . 
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Problem  1 .  The  System  H2O  and  the  self- ionization  of  water 

Consider  the  System  H2O  including  the  self-ionization  of  liquid 
H2O (water)  .  There  are  four  species: 

H20(ideal  gas) 

H2O (liquid,  i.c.  water) 

H+ (aqueous  hydrated  ion) 

OH- (aqueous  hydrated  ion) 

In  considering  this  chemical  system,  the  constants  for  H2O (ideal  gas) 
are  based  upon  an  earlier  run  and  are  not  further  refined  here;  the 
constants  for  H20(liquid)  ,  H+  (aq)  ,  and  OH“(aq)  are  sought.  Only  the 
following  portion  of  the  data  available  are  used  in  the  refinement  to 
keep  the  problem  in  manageable  bounds  for  publication: 

1.  The  equilibrium  constant  reaction  H20(liquid)  =  l^CKideal  gas) 

2.  The  relative  heat  content  data  for  liquid  H2O. 

3.  ACp  for  the  self-ionization  of  liquid  H2O. 

4.  AH  for  the  self-ionization  of  liquid  H2O. 

5.  Selected  equilibrium  constants  for  the  self-ionization  of 
liquid  H20. 

The  following  pages  give  the  listing  of  the  input  and  output  for  the 
regression  of  the  model  to  these  data. 

Listing  of  input  deck.  The  next  four  pages  contain  a  listing  of 
the  input  deck.  Deleted  from  the  input  are  the  constants  for  H20(G) 
which  were  in  the  format  (7A8)  . 
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0 

0 

(L  )H» 


1 


1 


DH- 


H2CILI  * 

H20(  G  1  .LOG 

K. 

HAAS, 

1970 

2  36 

5  0 

1  1 

H20  CGI 

1. 

CCOD  00 

0 

0 

H20  (LI 

-1  . 

OOOO  00 

- 

1 

0 

0.0100 

00 

27  3.  1500 

00-2. 21  97  CO 

00 

o.oobood 

00 

10.0000 

CO 

273.1500 

00-1.917200 

00 

0. 004000 

00 

20.0000 

CO 

273.  15uU 

OC-1. 637600 

00 

0.004  00D 

00 

30.0000 

03 

2  7  3.  1500 

CO-1.  3788C0 

00 

0.  C04000 

00 

40.  0000 

00 

2  7  3.  1500 

00-1. 139000 

00 

0.004000 

00 

50.0000 

00 

273.  1500 

00-0.  9161CD 

00 

0. 004  OOD 

00 

60.0000 

oc 

2  73.1500 

00-0.708500 

00 

0.0J400D 

00 

70.0000 

00 

273.  150C 

00-0.  5150C0 

00 

0. J0400D 

00 

80.0000 

00 

273.1500 

00-0.334100 

00 

0.C04C00 

00 

40.0000 

00 

273.  1500 

CO-G.  1649C0 

00 

0.0J4JJ0 

00 

100.0000 

00 

273. 15C0 

CC-O.  006300 

00 

0.  C0400D 

00 

110.0000 

00 

2  73.  1500 

00  0.142400 

00 

0.  GO 5 OOD 

00 

120.0000 

oc 

2  73. 15 CO 

00  0.2822CD 

00 

0. 00600D 

00 

130.  0000 

00 

2  73.  1500 

00  0.413800 

00 

0.007000 

00 

140.0000 

GJ 

273. 15C0 

CO  0.53  7600 

00 

0.  00800D 

00 

15C.OOCO 

00 

273. 1500 

00  0.654400 

00 

0.009000 

00 

160. GOOD 

00 

273.  15CD 

CO  0. 764 7 CO 

00 

0.010000 

00 

170.0000 

00 

273.1500 

00  0.868900 

00 

0. Cl  1  COD 

00 

lBJ.OCCJ 

CO 

273.  15C9 

CO  0. 9674C0 

00 

0.012000 

00 

190.0000 

GO 

273. 1500 

CO  1.C606CD 

00 

0.012000 

00 

2C0.  OCOO 

GO 

273.  1500 

00  1.148900 

oc 

0.0 130 JO 

00 

2 10.0003 

00 

273. 1500 

00  1.  2327C0 

oc 

0.  014000 

00 

22O.CO0D 

00 

273. 1500 

00  1.312200 

00 

0.01 5000 

00 

230. OCOO 

CO 

273. 1500 

OC  1. 3 8 86 CO 

oc 

0.  J15000 

00 

240.0000 

00 

273.1500 

00  1.459400 

00 

0.C16C0D 

00 

250. oocn 

00 

273. 1500 

OC  1.52  7600 

00 

0.016000 

00 

260. OGLU 

00 

273.1500 

OC  1.586000 

oc 

0.  C 1  76  JO 

00 

2  70.  OCOO 

00 

2  73.  1500 

CO  1.654200 

oo 

0.017000 

00 

2a:.ooco 

CO 

273.1500 

00  1.  713100 

00 

0.  018000 

00 

240.  0000 

00 

273.1500 

00  1  .  76  8  2  00 

00 

0.018000 

00 

300.0CCJ 

CO 

273. 1500 

CC  1.8227C0 

00 

0. J19C00 

00 

310.0000 

00 

273.1500 

00  1.873H00 

00 

0.  019 COD 

00 

320.0000 

Go 

273.  15C0 

CC  1.S227C0 

00 

0.  J19000 

00 

330.000U 

00 

273.1500 

00  1.969400 

00 

0.  020060 

00 

340.0000 

00 

273.  15C0 

00  2.014900 

00 

0.020000 

00 

350.0000 

1.0 

273.  15C0 

OC  2.C570C0 

00 

0.  C20000 

00 

H20  (LI. 

H(  T I-HC273. 

161 

OSBORNE, 

STINSON, 

AND 

2  73 

7  0 

0  0 

H20  (LI 

1 

.0000  00 

-1 

0 

H20  (LI 

-1 

.  OCCO  00 

-1 

0 

5.C000 

00 

273.1500 

00  5.0280 

00 

1.8C0660 

01 

0.  000500 

00 

10. OCOO 

00 

2  7  3.  1500 

00  1C.04C0 

JO 

1  H.J660 

01 

O.CJo!  JO 

00 

15. OCOO 

CO 

273.1500 

00  15.0430 

00 

L • 80066 J 

01 

C.  00  C  5  JO 

00 

20.0000 

00 

273.15C0 

00  2  0.0400 

oc 

1  .600660 

01 

0  Ch-  j v  00 

00 

25.0CC0 

00 

273.15CO 

OC  25.0330 

00 

1 . 600660 

01 

0. CJ050D 

00 

30.0009 

00 

273.  1500 

00  30.0240 

JO 

1  .8J06b0 

01 

J.0JC5 JO 

JO 

35.0000 

00 

273.  1500 

OC  35.0150 

00 

1. 800660 

J1 

0.0005  00 

JO 

40.0000 

00 

273.  1500 

00  40.0060 

00 

1.800660 

01 

C.00J500 

00 

ANO  G INNINGS,  1939 


O.OIOOCO  00 
0.010000  00 
O.OIOOCO  00 
0.0100CD  00 
0.010000  00 
O.CIOOCD  00 
0. 010000  00 
0.0  1  00  00  0  0 
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5  . 

.  DC  OD 

00 

5  J  > 

.0000 

00 

55, 

,OuOr' 

00 

6  J. 

1  DuOO 

0  0 

65, 

.CJOn 

OD 

70. 

,00CD 

00 

75. 

0  COD 

02 

83. 

,  JCD  J 

JD 

85. 

,  ODD  J 

JD 

90. 

,  c  -CD 

CO 

95. 

.  0  OC  0 

00 

100. 

,  jCOD 

CC 

105. 

.GOOD 

CO 

11C. 

00C0 

JD 

1 1  6  . 

1 0  JCD 

00 

123. 

oOC? 

DO 

125. 

OCCD 

CO 

13  0. 

,DDO" 

00 

135. 

OCOO 

JU 

150. 

ocon 

3j 

155. 

ojIO 

00 

150 

3  JC  J 

00 

155. 

coco 

DO 

160. 

occu 

jD 

1'65. 

0002 

00 

17D. 

GOOD 

CO 

175. 

OCOO 

00 

lfaj. 

JoOD 

00 

185. 

OCCD 

C  3 

190. 

0  CDIi 

0  3 

195. 

ODCD 

CC 

2  CO. 

OCOO 

OC 

205. 

00  00 

CO 

210. 

C20D 

CO 

215. 

0  *00 

OD 

220. 

occc 

30 

225. 

0000 

00 

230. 

OODD 

OD 

2  35. 

uCCO 

<-C 

25C. 

ooco 

00 

25  5. 

GOOD 

00 

250. 

3  CD" 

00 

265. 

COO" 

DO 

260. 

OCOO 

OD 

26  5. 

ODOC 

JO 

2  7C. 

OCCD 

CO 

2  75. 

ODCD 

00 

280. 

ODOO 

JO 

285. 

OCOO 

JC 

290. 

O  ' JD 

DC 

295. 

j  COD 

OC 

300. 

OCOP 

DO 

3C5. 

OODD 

DO 

310. 

DDCD 

CC 

315. 

ocon 

JC 

32C. 

DCOD 

OD 

325. 

DOOD 

DO 

330. 

000  J 

CC. 

335. 

DOOD 

DO 

350. 

0000 

00 

355. 

OODD 

DC 

273.1520  cc 
273.15.jD  DO 
273.  15CD  CO 
273. 15JD  00 
2  73.  1 5 CD  DO 
273. 15u0  CO 
273.15CD  00 
273.1 5uD  CD 
273.  1500  00 
2  7  3.  1 5 CD  OC 
273.15CD  DO 
2  73.  15  CD  CC 

273.1500  00 
27  3.  1  5 CD  DO 
27  3.  15CD  CO 

273. 1500  DO 
273.  1500  00 
273.150D  00 
273.  1500  DC 
2  73.15CD  Oj 
273.  1500  Ou 
273. 1 5uD  OC 
273.  151.0  DO 
27  3.  150D  00 
273.  1  5 jD  00 
273.15C0  DC 

273. 1500  DO 
273. 15CD  00 

273.1500  DO 
2  7  3.  1500  00 

273. 1500  DO 
273.  1500  00 
273.150D  CC 
273. 15C0  DO 
273.  1 5CD  CO 

273.1500  00 
273.  15 CO  DC 
273.15CD  DO 
273.  1500  02 
273.151.0  DO 
273.1 50D  00 
2  73.  1 5 CD  00 
273.150U  DO 
273.  15DD  00 

273.1500  DO 
273.  150D  DD 
273.  1  500  DO 
273.15CD  00 
273.  15CD  DO 
273.  150r  DO 
273.  15oD  DO 
273.15DD  00 
2  73.  1500  OC 
273.  15CD  OC 
2  7  3.  1  5 OD  00 
273.  15JD  00 
273.150D  00 
2  7  3.  1 5  CD  OC 
273.150D  DO 
273.  1500  OC 

273.1500  OD 


'.5.09  70  3  0 
59.9900  Oj 
55.  9050  OC 

59.90  10  00 
65.  90  ID  OD 
6S.S05D  00 
75.991D  30 
00.  0020  OC 
05.0101)  OD 
90.J39D  00 
95.067D  00 
IOC. 102D  CC 

10  5.  15  0''  OC 
113.  1900  00 

11  5.2503  OC 
12 0. 32 JD  00 

12  5.5003  JO 
1 3D. 59 CD  DO 
135  .  60  00  JD 
15G.72CD  JD 
155.  95 CD  JO 
151. COCO  OC 
156.16  ,D  OC 
161.  35 CO  OC 
Ufa. 5500  30 
1/1.77C0  00 
177.0100  OC 
10  2.  2HCD  00 
10  7.  5 7 CO  30 
1 92. 00 CD  00 
190.2300  DC 
2C3.59jD  33 
2C9.C3C3  OD 
215.5390  OD 
219.  93 CD  00 
22  5.5C3D  00 
23C.95cO  30 
236.  53C0  OC 
2*  2.  16011  30 
25  7.  05  CD  00 
263. 5733  00 
2  59.  35C3  JO 
265. 1900  OC 
271.09CD  00 
277. 0600  DC 
203.C90D  00 
269. 21  CD  DC 
295.51  3D  30 
301. 70 CD  00 
300. 10C0  OD 
315.60 CO  00 
321.220"  OD 
32  7.  9SCD  00 
335.69CC  30 
351.97CD  OC 
359.  25 CD  3D 
356.7300  Oo 
365.50CD  00 
372.  55 3D  00 
38C.93DD  00 
369. 7600  DD 


1.000660  D 1 
1.O0066D  01 
1 . 000660  01 
1.603660  Cl 
1.00066D  D1 
1. 000660  D 1 
1  .630660  01 
1.UjD66D  Cl 
1.80066D  01 

1.600660  01 
1  «  60  0660  01 
1.00D66D  Cl 

1 . 60  0  fa  60  01 
1.HJJ66D  01 

1 . 90  0660  01 
1.9G066D  D  1 
1.00066D  01 
1.60066D  01 
1.HC066D  01 
1.B3D66D  01 
1 *6u 066D  J 1 

1.600660  01 
l.B00u6O  01 
1.00066D  01 
1.60066"  Cl 
1.000663  D 1 
1.03066D  01 

1.90  0660  01 
1.80J6fc9  D 1 
l.HJGfaoO  D 1 
1.00C66D  D 1 
1 .000660  0 1 
1 . 000660  0 1 

1.90  0660  Cl 
1 . HO  06  6D  0 1 
1.H3D66D  Cl 
1 . 6 J0660  D 1 

1.600660  Cl 
1.60J66D  Cl 

1.900660  01 

1 .900660  D 1 
1.9C066D  D 1 
1  .6006  60  D  1 

1.900660  01 
1.900663  01 
l.OJQfafaD  01 
1.0J066D  01 
1  . fa 0 3  fa oD  01 
1 • 6 JC  660  D 1 
1  .800660  0  1 
1.0jO66D  01 
1.9JJ66D  01 

1.900660  D 1 
1.030660  D 1 
1.800660  01 
1.60066D  D 1 
1  .60066D  D  1 
1.  900660  D1 
1.30D66D  D 1 
1.00066D  01 
l.faODfafaD  D 1 


0.  9  32  5  3D  00 
D. 030500  00 
D.30C5uD  00 
C.  JJ050D  JC 
u •  000501)  00 

C.  3JC5CD  DC 
0.2005JD  00 
0. CJ050D  00 
0. 2005 JD  DO 
0.  00C5JD  JO 
0. C0C5CD  CD 
0.00050D  00 
0.CO050"  JO 

D. 0305JD  00 
C.  00  j  500  OD 
0.00050D  00 
C.  0C..5  0D  OD 
J.3J050D  00 
0. jj j50D  OD 
C.3CC50D  CO 
J.00O50D  JO 
C. CJ050D  OJ 
0. 32 0500  OC 
c.jooson  jo 
C. J3  3600  00 
2. u0«5  30  30 
C. Oj35jD  00 
0.3jj5JD  03 
C.  232500  OC 
O..Jj5jD  30 
0. 000500  00 
0.03050D  30 
O.33O530  30 
J.0JJ533  30 
o.o:o5or  30 
C. 0 Jo50P  20 
0 . 002  5  OP  JD 
0.  J0C5CP  00 
0.2JC5JD  30 

C.  3305  JD  JD 
J.0035JD  00 

D. C00502  00 
O.OOD5CO  DO 
0.030500  JO 

C.  CC050D  00 
0. 02050D  03 

D. CD050P  3D 
0.20050P  00 
D. 300500  OD 
D.0J05CD  00 
0.30C5JD  JO 
C.C0C5CD  DD 
0.JJC5OD  00 
0.  OOC50D  00 
O.OCJ5UD  DO 
0.  00050"  00 
J.00050D  DC 
0. OC  050C  JC 

C. 33C50 0  00 
0 . 000 500  3D 

D.  uOC5CD  00 


0 •  01  GOOD  DO 
D.D10DCD  00 
0.0100CD  DO 
0.  Cl  00 CD  00 
D  •  0 10000  00 

C. D10DCD  00 
0.C1000D  00 

D. 010C0D  00 
D.0100CD  00 
D. 01 DOOD  DD 
D.0100CD  00 
0.01000D  00 

C.  D 1  00 CO  00 
0.010000  00 

D.  D10DCD  00 
D.D100CD  00 
D.D1000D  00 
D»  Cl  DO  CD  00 
D.D1000D  00 
D.0100CP  00 
C.D1000D  OD 
0.  C1DDCD  00 
C  .  0 100 JD  00 

C. OICODD  00 
O.UlOCCD  DO 

D. CluGDD  00 
O.OIDOOD  00 
C.C1000D  CO 
O.CIODCD  00 

C. D10CCD  00 

D. DICOCD  00 
O.CIODOD  DO 
O.OIJOOD  00 
O.OIDOOD  00 
D.01C0CD  00 
0.  Cl  GOOD  00 

C. 01000D  00 

D. C10CCD  00 
D  •  0  1  GOOD  00 
D.OIDOOD  00 

C. D10DCD  00 

D. 01JCDD  00 

C. 0100CD  00 
O.D1COOD  00 
O.DIDDCD  00 

D. D1000D  00 
O.D1COCD  00 
D.D1D0DD  00 
D.OIDOOD  00 
0.D100C0  00 
D.OIODOD  00 
D.C1000D  00 
D.010DCD  00 
0.  C10DCD  00 
0.0  10D00  00 
D. 01 ODOO  00 
D.01000D  00 
0.010000  00 
O.CIODOD  00 
C.D100DD  00 
0.0100CD  00 


350. COO!)  CO  273.  1500  CC  395.17:0  00  1 

355.0000  00  273.15C0  00  5C9. 36CD  00  1 

360.0000  00  273.  15C0  CO  520.7500  00  1 

365.0000  00  273.  1500  00  535.12C0  00  1 

H*  ♦  OH-  -  H20(U.  OFLTA  CP.  FACTION. 

3  5  10  11 

H20  (LI  1.0000  00  -1 

H*  -l.OOOD  00  -1 

OH-  -1.0000  00  -1 

lO.OCCO  00  273.  1500  OC  65.2C0C0  00 
30.0000  00  273.1500  CC  50.10000  00 
50.00C0  00  273.  1500  CO  51.80000  30 
7J.000J  00  273.  1530  CO  52.300CO  00 
H20  =>  OH—  ♦  H  ♦  ,  HEAT  CF  I  ON  I  7  AT  ION 


.1)00060  Jl  0.000500  00  0.010000  00 
•  80J660  01  C. 000500  00  0.01C0C0  00 
.8OO60D  01  0.000500  00  0.010000  00 
.830660  01  0. C005 wO  00  O.CIOOCO  00 
ACKFR«ANN,  1958. 


3.000000  00 
3.000000  00 
3.330CC0  00 
3.000000  00 

OLOFS SON  ANJ  DLOFSSON,  1973. 


3  6  3  1  0  1 

H20  (LI  -1.0000  00  “1 

H*  1.0000  00  -1 

OH-  1  .0000  00  -1 

298.  1500  CO  5,5.9  2000  0  32. 

323.1500  00  5C.92CCD  032. 

35  7.  5500  00  56.670./D  0  32. 

373.5500  00  52.3  7039  032. 

398.5500  00  37.35000  03?. 

517.7500  00  33.550CD  032. 

LOG  KW 

3  25  50  1  1 

H20  (LI  -1.0000  00  -1 

Ht  1.0000  00  -1 

OH-  1.0000  00  -1 

51.00C0  00  273.1500  00-13. 251C0  00 

60.9CCD  00  273.1500  00-12.99190  00 

70.50CC  00  273. 15 CO  00-12.80160  00 

80.2C0D  00  273.1500  00-12.651  70  00 

89.9C00  00  273.15CO  CO-12. 590CD  00 

99.30C3  00  2  73.1500  00-  12. 23'- 30  00 
1C3.6C00  00  273.  1500  00-12.07550  00 
118.0000  00  273.1500  00-11.96650  00 
127.3C00  00  273.1500  00-11.86990  00 
136.700  0  00  273.  1500  00-11.76720  00 
156.00C0  00  2  73.1500  00-11.72030  OC 
155.5CCC  00  273.1500  00-11. 60760  00 

165.7000  00  273.1500  00-11.55770  00 
175.20CO  00  273.  15LC  OC-1  1.57280  OC 
183. 60C'  00  273.1500  CO-11.39290  00 
193.000)  00  273.  1500  OC-1 1 . 332  6  0  00 

202.5000  00  273.1500  CO-ll. 279CO  00 
212. 1  COO  00  273. 150n  00-11.25860  00 

221.7000  OC  273.  1500  00-11.22830  00 

231.5000  00  273.1500  00-11.21550  00 
251.2C0J  00  273.  1500  CC-11. 20210  00 
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Il-P  INDEX  CF  Th„  LAST  ITEM  OF  THE  ABOVE  DATA  SETS  IN  THE  VECTORS  XII, II.  YCII),  ANC  SlGVOtll 
IS  AS  FCLLCwS : 


NljPBEP  CF  lNOtPcNJcllT  V&r.I4eLcS  PER  OBSERVATION  IS  2 
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AGREEMENT  FACTI.RS  BASIC  ON  PARAMETERS  BEFORE  CYCLE 


SYiTfcM  H2C  AN  J  THE  SELF-  I OM  t  A 1 1  ON  CF  WATER 
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SYSTLM  H2C  AND  The  SELF- I  CM/ AT  ICN  OF  WATER 


Tf-ESE  PESU.7S  wEAt  J8TAINFC  IN  A  BIN  CN  10/31/73 


THESE  SESLl  TS  »eRi  UBTAINEC  IN  A  (HA  CN  10/31/73 


THESE  RFSllTS  WERE  OfaT R I  NEC  If* 


Problem  2.  The  System  Ni-0  and  the  thermodynamics  of  Bunsenite  (NiO) ■ 

As  a  second  example,  consider  the  selected  calorimetric  data  for 
the  following  phases  for  the  System  Ni-O: 

T<631  K  rie  pernt) 

T>631  K 

T<525  1  ^lurie  pernt) 

525  K<T< 565  K  (Transition  point) 

T> 565  K 

ideal  gas,  all  T 

The  constants  for  nickel  and  oxygen  are  derived  from  an  earlier 
regression  and  not  refined  in  this  problem.  Only  refined  are  the 
constants  for  Bunsenite,  NiO(B),  and  NiO(C) . 

The  data  are  the  results  of  studies  by  E.G.  King  and  his  associ 
ates  and  include: 

1.  Cp  (Bunsenite)  from  low  temperature  calorimetry 

2.  S2g8 (Bunsenite)  from  third  law  evaluation  of  low 
temperature  calorimetry. 

3.  H (T) -H ( 298)  of  Bunsenite  from  drop  calorimetry. 

4.  H (T)  of  NiO (B)  -  H ( 298)  of  Bunsenite  from  drop  calorimetry. 

5.  H (T)  of  NiO (C) -H (298)  of  Bunsenite  from  drop  calorimetry. 

6.  Ah°  of  Bunsenite  from  combustion  calorimetry. 

The  input  and  results  for  this  problem  are  given  below. 

Listing  of  thput  deck.  The  next  two  pages  contain  a  listing  of 
the  input  deck.  Deleted  from  the  input  are  the  constants  for  Ni (A) , 

Ni  (B)  ,  and  0^ (G)  which  were  in  the  format  (7A8) . 


—  There  is  no  first-order  structural  change  at  the  Curie  pernt. 
However,  because  Cpi  is  discontinuous,  the  mathematical  model  requires 
separate  algebraic  expressions  for  the  temperature  regions  above  and 
below  the  Curie  pernt.  The  Curie  pernt  (and  also  all  similar  phenome¬ 
non)  are  treated  here  as  "first"  order  transitions  with  mathematically 
derived  AH  and  AS  of  "inversion"  to  describe  the  local  anomaly. 
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THE  SfSTEM  NI-Q  AND  ThE  T 
0  0  1 
0 

(BIBUNSENI  TNIO  (8JN10 


10/31/73 

PROBLEM  2,  THE  SfSTEM  M 

13  0  0 

6  0  0 

N1  (4INI  (BIBUNSENITt 
1  0  0  0  0  1 

BUNSENITE,  Ht  AT  CAPA  CITY 

1  10  1  1 

BUNSENIT  1.0000  00 
206. 13D  00 
216.190  00 
225.850  CO 
236.070  00 
245.630  00 
256.290  00 

266.150  00 
276. 01D  00 
286.430  00 
295.940  00 

BUNSENITE,  ENTPCPr 

112  1 
bOnsenite  l.OOCD  00 

298.150  00 

BUNS  tN IT  E  >  RELATIVE  HEA 
2  5  7  1 

BUNScNITr  -l.OOCD  00 
BUNSENITE  1.0000  00 

364.900  00 
433.000  00 

502.900  00 

513.200  00 

524.200  00 


HEBMuDTNARICS  OF  BUNSENITE. 


<8)910  I C 1 02  10) 


8.13400  00 
8.45700  00 
8.76200  00 
9.05000  00 
9 . 306 CD  00 
9.59700  00 
9.84500  00 
10.0  80  00  OC 
10.32000  30 
10.55000  00 


9.0  n  JOQ  00 

lathe  heat  content. 


0.7500D  03 
1.58000  03 
2.5850D  03 
2.74000  03 
2.92500  03 


NICKEL  C  X  IDE  t  RELATIVE  HtAT  CONTENT 
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Printed  output .  The  following  pages  contain  a  complete  printed 
output  from  the  execution  of  PHAS20  using  the  preceding  data  deck. 

By  setting  IL  to  0,  the  results  of  the  intermediate  refinements  were 
deleted.  Constant  18,  corresponding  to  the  "d"  constant  for  bunsenite, 
is  the  least  significant  of  the  constants  included  in  the  regression. 

If  ICY  had  been  set  equal  to  or  greater  than  1,  this  constant  would 
have  automatically  been  set  to  0.0  and  the  regression  rerun. 
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OOOriOOOOOOOOOO 


PHAS20 


THIS  VERSION  *ILL  HANOLE: 

20  PHASES 
70  OATA  SETS 

ANC  1200  INDEPENDENT  C8 SER VAT  IONS . 
VERS  I CN  TESTEO  8/15/73  8V  HAAS. 


- PHA 

PHA 

PHA 

PHA 

PHA 

PHA 

PHA 

PHA 

PHA 

PhA 

PHA 


IMPLICIT  REAL*8< A-H.O-Z  )  PHA 

LOGICAL*!  LABEL, IMAGE, 8C0  PHA 

REAL*'.  XMAX,XMIN,YMAX,YMIN,XI  ,E  ,  T I TLE  ,  ACGEF  ,  T,  TO,  A  VAL  ,  Rt  L  ER  R  ,  SERR,  PHA 
1  ScRRSU,C6AR,STCEV,ERR,AN  PHA 

01  PENSION  AVAL  17 1  ,SSRP(3! ,SERRSC(3) ,E8AR (31 , STOcVI 3 » , IMAGE!  5000) ,  PHA 
1  XK1I.E  (ll.PHINVtfel  ,VAL(7),IdN(20),IF«IN(20i  PHA 

01  PENSION  CO HP (6, 7  0) , PNAMt (20  I, TINVt  6 , 70,41 , IPHASE ( 6, 70) ,  PHA 

1  NPHASEI 701 , IKOLNT! 70) ,IG0(70) , ISTATC! 6 , 70 ) ,NI NVER ! 6, 70 ) ,  PHA 

2  INSTAT! t, 701, INVPHI6, 70,51, INVSCC6, 70)  PHA 

oi Pension  x <2.1200) ,vc(  12001 .sigtoi 1200) ,pi iaoi.kii iaoi,  pha 

1  OC  (  I'iO  ) , POI 1  AO  I ,  f  I  TLE  ( 20 1  PHA 

01  PENSION  REF! 10, 70) , ERRP! 1AGI , I HOL  Q( 40 )  PHA 

01  PENSION  SCINV(2),STCDEF( 101 , YCSN 3! 2) , TK C 19 ,2 1 ,  TYPE!  14 1 , NSC  ALE  I  5) PH A 
1  .LABEL ( 50 )  PHA 

01  Mc  NSI  ON  or  cc  (  7  )  ,  OV OC  (  7 )  PHA 

Cl PENSI ON  A/I140)  PHA 


COMMON  BLOCKS 
NAPE 
FARTH 
AIR 


ROUTINES 

MAIN  PROGRAM,  0RGLS2,  EAF-20,  PUTDUT 

MAIN  PROGRAM,  0RGLS2,  TEST20,  FAFH20 

PUTDUT  ANO  UNICUE. 

MAIN  PROGRAM,  OPGLSZ,  TEST20,  EAFWZO 

MAIN  PROGRAM,  eafWZO,  VOERIV,  ODER  IV, 

PUTOUT,  FN,  OIE,  ANO  PLCCK  OATA. 

MAIN  PROGRAM,  PLTOUT 

MAIN  PROGRAM,  EAFW20,  YOERIV,  OOERIV, 
ANO  PUTPUT 

MAIN  PROGRAM ,  PUTOUT 


COMMON  /EARTH/  COEF.PNAPE.TINV, IPHASE, NPHASE, IKOUNT, I  GO, NSET S , 

1  I  STATE,  MOVER,  IN  ST  AT,  INVPH,  INVSC  ,LISTP 

COPMCN  /AIR/  X,  P,  OC,  TITLE,  YO,  SIGYO,  PO,  Kl,  NC ,  NV,  NX,  IW, 
l  NP,  NO,  ISING,  I  STOP ,  IL,  ON 

COMMON  /FIRE/  REF,  ERRP,  IWRITE,  ICY,  IICY,  IRECC,  IREG,  NHOLO , 

1  I  HOLD 

COMMON  / HATER/ ZERC ,  ONE,  TWC,  THREE,  FOUR,  SIX,  R,  F, 

1  SCINV,  TREF,  STCOEF,  OIEO,  ADIE,  BCIE,  THETA,  YESNO,  TK, 

2  ASTAR,  TYPE,  NL ,  NSC  ALE,  LABEL ,  8C0 
COPMCN  /TI“E/  OATE 

COPMCN  /SPACE/  OCCC,  OYOC,  SC,  TLOW 

COPMCN  /PAN/  AA 

EXTERNAL  EAFW20 

CALL  ERRSET  (208,256,0,2) 


REAO  IN  ANO  STORE  INPUT  OATA — 

OATE  --  TOOAY'S  OATc 

NREG  --  NUMBER  OF  SEPARATE  PROBLEMS  IN  THE  DATA  SET 


oooonnoooonoooooor-oor  o  r>  r 


RE  AO  (5.44C0)  OATE 
RE  AO  (5,45CU)  NREG 
00  4300  IPFGM.NREG 

INITIAUI2E  CONSTANTS 
WRITE  OAT f 


nCY-o 
K0UNT*J 
IWRlTE*l 
WRITE  (6,5300) 

WRITE  (6,5400)  OATE 


READ  IN  AND  STCRE  — 

TITLc  (20)  —  title  for  correlation 
IDO  --  1  FOR  REGRESSION  AMO  ERROR  PLOTS 

2  FOR  REGRESSION  ONLY 

3  FOR  ERROR  PLOTS  ONLY 

NC  —  NUMBER  OF  CYCLES.  GcNFRALLY  2  IS  SUFFICIENT. 

SUBROUTINE  TFST  WILL  TERMINATE  REGRESSION 
IF  THt  RELATIVE  CHANGE  IN  SUCCESSIVE 
PARAMETERS  IS  LESS  THAN  THAN  1.00-8. 

IW  —  0  FOR  WEIGHTEO  OATA 

1  FOR  UNWEIGHTEO  OATA 

IL  „  q  for  CALCULATED  RESULTS  AT  PLOT  TIME  ONLY 
1  FllR  CALCULAfEO  RESULTS  with  EACH  CYCLE 
IFMOLT  —  PUNCH  FORMAT  for  OUTPUT  OF  PARAMETERS— 

0  GIVES  PARAMETERS  IN  6012.5  FORMAT 
1  GIVES  PARAMETERS  IN  6A8  FORMAT. 

LISTP  —  NUMEEP  OF  PHASES  IN  CORRELATION 

ICY  —  NUMBER  OF  STEPS  IN  STFP-8ACKWARC  ELIMINATION. 

DEFAULT  IS  0. 

NHOLC  —  NUMeER  OF  PARAMETERS  TC  BE  KEPT  WITHIN  THE 

REGRESSION  EVEN  THOUGH  'TEST'  WILL  INOICATE  THEY  ARE 
NON-SIGNIFICANT. 

PNAMEILISTPI  —  S-CHARACTIR  LABEL  FCR  EACH  PHASE  IN 
CORRELATION 

I  ONI  LI STP I  —  FLAG  WHICH  INDICATES  THE  SPECIES  IS  A 

REFERENCE  ELEMENT  (II,  A  COMPOUND  (0),  OR  AN  ION  (-1). 
NSETS  --  NUMBER  CF  OATA  SETS  IN  CDRFELATICN 


REAO  ( 

READ  ( 

PE  AO  ( 

REAO  ( 

READ  ( 

REAO  < 

WRITE 

WRITE 

WRITE 

WRITE 

WRITF 

°UNCH 

PUNCH 

PUNCH 


5,4600  I 
5.4EC0) 
5.45C0)  l 
5 , 4  4C0 I 
5,47001 
5.45C0) 
(6,55001 
(6,5600) 
(6,53001 
<6,57001 
(6, 5800) 
460C,  (T 
450C,  LI 
440C ,  (P 


(TITLE! 1 1,1*1 ,201 
IOC,NC,IW,IL,IFMOUT 
LISTP, ICY.NHOLO 
(PNA"E(M|,M=1, LISTP) 

(  I CN ( I  1, 1*1, LISTP) 

NSETS 

(TITLE! 1 1  *  1  =  1 *201 
( PNA«E(MI ,M*1, LISTP) 

I 

i  (TITLE!  II, 1*1, 201, OATE 

I 

riTLEl 11,1*1,201 

[  STP 

>NAME(H),M*1, LISTP) 


REAO  IN  ANO  STORE  INPUT  OATA  FUR  EACH  OATA  SET  — 

RtF  (10, J)  —  REFERENCE  FOR  'J'TH  DATA  SET  (OATA  SET 
IDENTIFICATION) 

N PHASE ! J)  —  NUM8ER  OF  PHASES  IN  REACTION,  SET  TO  l  WHEN 


-PHA  610 
PHA  620 
PHA  630 
PHA  640 
-PHA  650 
PHA  660 
PHA  670 
-PHA  680 
PHA  690 
PHA  700 
PHA  710 
pha  no 
PHA  730 
-PHA  740 
PHA  750 
PHA  760 
PHA  770 
PHA  780 
PHA  790 
PHA  800 
PHA  310 
PHA  820 
PHA  830 
PHA  840 
PHA  850 
PHA  360 
PHA  870 
PHA  880 
PHA  890 
PHA  900 
PHA  910 
PHA  920 
PHA  930 
PHA  940 
PHA  950 
PHA  960 
PHA  970 
PHA  980 
PHA  990 
PHA1D00 
PHA 10 10 
— PHAIOZO 
PH  A 1030 
PHA 1040 
PHA 1050 
PHA 1060 
PHA 1070 
PHA1080 
PHA 1090 
PHA 1100 
PHA1110 
PHA 1120 
FHA1130 
PH  A 1 140 
PHA1150 
PHA1160 
— PHA l 1 70 
PHA 1 1 80 
PHA 1 190 
PHA1200 
PHA  1210 
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DATA  FOP  CNE  PHASC  ONLY. 

IK3UNTIJI  —  NUMBER  OF  OBSERVATIONS  IN  DATA  SET 
IGOIJI  --  DATA  TYPE » 

FOR  HEAT  CAPACITY,  IGOIJI  IS  I 
ENTROPY  2 

ENTHALPY  3 

GIBBS  ENERGY  4 

LOG! 101  K  5 

CELL  OATA  6 

DIFFERENTIAL  HEAT, 

HIT  )-  HI  TO  )  7 

FOR  ANY  DATA  NOT  INCLUDED  IN  TYPES  1 
THROUGH  7,  SET  IGOIJI  TO  ANY  INTEGFR 
8  THRU  14.  SUPPLY  SUBROUTINE  UNIQUE 
TU  CALCULATE  YC  (ANC  THE  DERIVATIVES 
OCII)  IF  JN  I  OR  JDFLAG)  =  ZERO  FROM 
THE  CURRENT  VALLES  Cp  THE  PARAMETERS. 

IF  MORE  THAN  ONE  TYPE  OF  DATA  AFE  TC  BE 
CONSIDERED,  THESE  ARF  KEYED  BY  01 FFERFNT 
VALUES  OF  IGOIJI. 

ITFACT  —  TFE  VALUES  IS  1  IF  TEMPERATURE  IS  IN'  KELVINS, 

0  IF  THE  TFMPERATURE  IS  IN  CENTIGRAOE  DEGREES. 

IPARA  --  PROGRAM  IS  SET  F JR  UNITS  IN  JCULESI*!,  VOLTS,  CR 
LCGI10IK  IN  MOLAL  UNI TS. 

IF  SO,  SET  IPARA  TO  1 

IF  NIT,  T  IPARA  TO  J.  PROGRAM  MULTIPLIES  YOU)  ANO 
S I  GY  J 1 1  I  BY  PARA  TO  BE  R  FAD  LATER. 

I S I G  —  D  IF  ^FIGHTS  Arc  RELATIVE. 

1  IF  WEIGHTS  ARE  ABSOLUTE. 

•  ABSOLUTE  WEIGHTS  ARE  USEC  IN  CALCULATIONS.  I 

(*l  FUR  CATA  IN  CALORIES,  AOJUST  R  ANO  F  IN  BLOCK  DATA.  SEE 
COMMENTS  IN  FLOCK  OATA  wHICH  CESCRIBE  TFE  NECESSAPY 
ADJUSTMENTS. 


CD  laOJ  J*1,NSETS 

READ  4  5,4400)  IREFI I , J) ,I>1 , 10  I 

WRITE  16,5900) 

WRITE  16,40301  ASTAR 
WRITE  16,59001 

WRITp  I  6 , 6 100 )  J, IREFI  I  ,  J  I  ,  I  =1,101 

READ  15,45001  NPHASE 1 J I , IKUUNT I J  I ,  IGOI J  I  ,  ITFACT, IPARA,  I  SI C 
WRirt  (6,6200)  NPFASEIJI, IK  CUNT! J)  .TYPE!  IGOIJI I 
NPFA S*NPHASE(J I 

FOP  EACH  CATA  SET,  REAO  ANO  STORE  F CR  CNE  PHASE  AT  A  TIME, 

ANAME  —  8-CHARACTER  PHASE  LABEL.  MUST  BE  THE  SAME  AS 
GIVEN  IN  PNAWEIIISTP). 

COEF  1 1 ,  J I  —  COEFFICIENT  IN  REACTICN  tU.000  IF  NOT 
REACTION) 

I STATE1 I , J I  --  FOR  ELEMENTS  IN  LOWEST  TEMPERATURE  MODIFICA¬ 
TION,  SET  TO  1 

FOR  COMPJUNOS  ANO  ELEMENTS  IN  HIGHER  TEMPER' 
ATURE  MODIFICATIONS,  SFT  TO  0 

NINVERI I ,  J  I  —  NUMBER  OF  I NVER S I CNS I  EQUALS  NUMBER  OF  LTWFR 

TEMPERATURE  MCC I F ICAT I CNS  BEING  CONSIDERED 
FOR  COMPONENT.) 


NINVERI I , J )  — 


00  700  I «1,NPHAS 

REAO  I5.4FC0)  ANAME, COEFII, J ) , I  STATE  1 1 , J)  ,NINVER(I,J) 


PHA  1220 
PFA 12  30 
PHA 1240 
PHA 12  50 
PHA1260 
PHA 1270 
PH A 12  80 
PFA 1290 
PHA 1300 
PHA 13 10 
PHA 1320 
PHA1330 
PFA 1340 
PFA 1350 
PFA1360 
PFA 13 70 
PHA 1330 
PFA  1 390 
PHA  1400 
PFA 14 10 
PFA 1420 
PHA 1430 
PHA 1440 
PFA1450 
PFA 1460 
PF  A  14  70 
PHA 1480 
P  FA  1490 
PFA1500 
PFA  1510 
PHA 1520 
PFA 1530 
PFA 1540 
PFA1550 
PHA1560 
•PHA1570 
P^  AI580 
MF41590 
PFA1600 
PF  A 161 0 
PHA 1620 
PHA1630 
PHA 1640 
PhA 1650 
PFA166P 
■  PHA 16  70 
PHA 1680 
PHA1690 
PHA 1 700 
PF41710 
PHA 1720 
PFA 1 730 
PFA1740 
-  PFA  1750 
PFA 1760 
PFA1770 
PFA17B0 
PHA 1790 
-PHA1B00 
PHA 1810 
PFA 1820 


c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 


c 

c 

c 

c 

c 


LOCATE  AN  APE  IN  LIST 
IPHASFCI, J). 


pnameilistpi  and  asign  the  index  to 


CC  ICO  K*I,LISTP 

IF  ( ANAME.EO.PNAPEIK It  GO  TO  200 
100  CONTINUE 
GJ  TC  A  IOC 
200  I  PHASE (I ,  Jl  *K 

NOYE  $= ( I S  7ATc(  I,JI+3)/2 
ACrEF=COCF(I,J( 

WPITP  (  6  «  A  300)  PNAMF (  I  PHASE!  I  ,  J  )  I ,  ACCEF  ,  YE  SNC 1  NCYES  » t  N INVE1’  (  I ,  J) 
IF  ( NINVERI  I,J).EC.O)  GC  TO  700 

IF  (MNVFRII  y  J I  .GT.O  I  RFAD  ANO  STORE, 

INST  AT ( I y  J I  --  FOR  ELEMENTS, ASIGN  l 

FOR  COHPOIJNCS,  ASIGN  0 

INVSClIyJ)  —  IF  UPON  INVERSION  THFRF  IS  A  STC ICH I CMSTO Y 
CHANGE  (EG—  2  F6CL3  =  FE2CL6),  THEN  FIND  THE  APPRO¬ 
PRIATE  REACTANT  COEFFICIENT  (HEPfc  »2*>  IN  THE  VECTOR 
STCOFF  ANO  ENTER  HERE.  DEFAULT  GIVES  A  COEFFICIENT  OF 

•1’. 

TINVININVER!  I,  JM  —  INVERSION  TEPPEPATUPES,  KELVINS 
PHINV1N  INVER  I 1 y  J I  *  1 1  --  PHASE  NAME  FOR  ALL  POD  I F ICATI HNS , 

BEGINNING  WITH  LOWEST  TEMPERATURE  PODIFICA 
TIJN  FIRST.  NAPE  POST  CCNFORT  WITH  NAMES 
IN  “NAME ( L I STP ) . 

kinver=n:nvfr(  i.ji 

READ  (5, (.SCO)  ITINV! I , J ,K I ,K=1 yKINVFR ! 

KlNvcR=MNvrP»l 

REAd  (5.5CCO)  INSTAT ( I  , J| , INVSC I  I , J I , I  PHI NV(K) ,K  =  1 yKINVSR ) 

IF  (  INVSClIyJ). EC. 0)  INVSC!  I,  J  I  ■=  1 

LOCATE  PFINVININVER!  I , J 1*1)  IN  LIST  PNAME(LISTP)  ANO  ASIGN  THE 
INDEX  TO  INVPH(lyJ). 

CO  500  K* I yK INVER 
CC  3CJ  LM.LISTP 

IF  ( PHINVIK I . EO.PNAMEtL  ) )  GO  TO  400 
300  CONTINUE 
GO  TO  A  IOC 
400  :NVPH( I, J,K)*L 
500  CONTINUE 

KINVER=K1NVER-1 
WRITE  (6,64001 
ACCE  F*  1. 0 
CC  6C0  K-I.KINVER 

IF  ( K.EQ. KINVEP )  ACOEF»STCOFF! INVSC! lyJII 
T»TIHV(I#J,K) 

WRITE  (6,6500)  ACCEF, PNAME( 1NVPH! I , J, K ) ) , PNAPE ( INVPH! I , J ,K* 1 ) ) ,T 
600  CONTINUE 

IF  (  INSTAT(tyJ).EC.l)  WRITE  (6,66001  PN'APEI  IN'VPHII  ,  J,  1 ) ) 

703  CONTINUE 

INITIALIZE  KO,  KCUNT , ANO  IKOUNT(I).  AFTER  UATA  IS  STORLO, 

I KUUNT ( I )  CONTAINS  THE  LOCATION  OF  THE  LAST  ITEM  IN  THE 
' J’TH  OATA  SET. 

K0=KCUNT*1 

KCUNT»KOUNTMKOUNT(  J) 


-  PHA 18  30 
PHAI8 40 
PHA 1850 

-PHA1860 
PFA1870 
PHA 1880 
PHA 1890 
PHA 1900 
PHA  1910 
PFA1920 
PHA 1930 
PF  A 19  40 
PHA 1950 

-  PFA 1960 
PFA  1970 
PFA1980 
PFA  1990 
PFA2J00 
PHA2310 
PHA2020 
PFA2030 
PHA  2040 
PHA2050 
PFA2060 

-PFA2070 
PFA2080 
PHA2090 
-PFA2100 
PFA2110 
PHA2120 
PFA  21 30 
PFA2140 
PFA  2150 
-PFA2160 
PHA  21 70 
PFA2180 
-PFA2190 
PHA2200 
PFA2210 
PFA2220 
PHA 2230 
PFA2240 
PHA2250 
PHA2260 
PFA2270 
PHA2280 
PFA2290 
PFA2300 
PHA  2310 
PHA2320 
PHA  23  30 
PHA2340 
PHA2350 
PHA  2360 
•PHA2370 
PHA2380 
PFA2390 
PHA2400 
•PHA2410 
PHA2420 
PHA2430 
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lKCUNTl  J)  *KCUM  PFA  2440 

WR  1 T  C-  (6ft  7 JO )  PFA2450 

WRITE  (6,68001  TYPE(IGOU))  PHA2460 

C  - PHA2470 

L  RFAO  ANO  STORE  TFE  DATA  (M  THE  *J'TH  CATA  SET,  PFA2480 

C  Xll.II  —  TEMPF  F ATUP  fc  (SEE  TFACT  RELCWI  PFA2490 

C  TFACT  —  O.OCO  IF  TEMPERATURE  IN  KELVINS  PFA2500 

C  273.1500  IF  TEMPERATURE  IN  CLNTIGRA05  CEGRFES  PFA2510 

C  VO  ( I  I  --  Ct  FENCANT  00SF.RVAT  1 CN  PF42520 

C  PARA  —  CONVERSION  FACTOR  TC  CCNVCRT  DATA  TO  CALn=I£S  C°  PHA2530 

C  V  US  (1.002  FJR  CATA  IN  K ( LOC ALORI t S ,  PFA2540 

C  1.00-3  FOR  CATA  IN  MILLIVOLTS,  —OR  THF.  PF425S0 

C  APPROP I  AT  2  CCNVeRS  ION  FACTOR  FG5  CATA  IN  PFA2560 

C  JOULES,  OTU/LB,  CUAR T»STCNL S/ACRE ,F TC . I  PH  A  28  70 

C  SIGYCIII  —  SEC  PROGRAM  0:SC»IPTICN  F0=  (SIGNING  -RIGHTS.  PFA  2580 

C  X ( 2 ,  1 1  —  FCR  CATA  TYPES  ( I  GO ( J I »  FRCM  l  THRU  8,  PFA2590 

C  X ( 2 , 1  1  »  7HRC.  PHA  26 JO 

t  FOR  RELATIVE  HE  AT  CONTENT  CATA,  XI2.II  IS  THE  EASE  PFA261J 

C  Tc  WP. RAT  URL .  Cr  r (LL  T  IS  298.15  <•  FFA2620 

C  FCR  DATA  TYPES  TO  PC  PPOGRAP"Fr  BY  THE  JSlR.  PFA2630 

C  XI2.II  IS  AS  OPTIONAL  SECCNC  I NOEPENCENT  PF72640 

C  VAR(A6LL.  PHA2650 

C  - - - PFA2660 

CO  I7JJ  1=KP,KGUNT  PH  A  26  7  0 

PEAU  (5, 4S00 )  X(  1,  I  i,  TFACT, YC( I » ,PARA,SIGYC(  (I ,X(2, 1 1  PFA2600 

IF  ( (1«.EC.3I.AN0.(  SIGYul  (  I .EU.O.UDOI.ANC .( I SIG.cC.J)  1  S I GYC  (  I  I  * ( . F F A 2o90 
IOCO  PFA27G0 

IF  ( ITFACT.tC.ll  GC  TC  1000  PFA2710 

c  - PMA2  720 

C  IF  (TFACT  *  Q  X  (  1 ,  I  )  -  X(1,I(  ♦  TFACT  PHA2730 

c  - PFA2740 

X(  1 , I  I « X ( 1 , I  I »TF ACT  OFA2750 

IF  (IGK  Jl.rc./I  CO  TO  800  PFA2760 

XL  2,11*0.  COO  PFA27  70 

GO  TC  900  PHA2780 

800  X(  2,  n  =  XC,I»*TFtCT  PFA2790 

900  CONTINUE  PFA2803 

1000  CONTINUE  PHA  2810 

IF  (IPARA.EC.il  GC  TO  11Q0  PHA2820 

C  - PFA2830 

C  IF  IP4RA  »  G  VO ( I  I  *  YQ ( I  I *PA R A  PHA2840 

C  - P  FA  28  50 

YC(I l*YO(  I»*PARA  PHA 2860 

1100  CONTINUE  PHA2870 

IF  (IW.EQ.l)  GO  TC  1300  PFA2880 

IF  (ISIG.RQ.il  GO  TO  1200  PHA2890 

C  - - PH  A  290  0 

C  IF  ISIG  *  Q  SIGYUII ) «  0A8S  ( Y0( I  II » S ( GYO(I)  PFA2910 

C  - PHA2920 

SIGYOII l*CABS(  SIGYOI I )*YO(I  > I  PFA29  30 

GC  TC  140C  PHA2940 

1200  CONTINUE  PFA2950 

C  - PHA2960 

C  IF  IPAKA  *  0  SIGYOd  »  =  SIGYO(  I  1*PARA  PHA2970 

C  - PFA2980 

IF  (IPARA.EQ.OI  SIGYOI I l«0A0S( S I GYOI I )*PARA)  PFA2990 

GC  TC  14GC  PFA  3000 

1300  CONTINUE  PHA30I0 

SI  GYC  ( I  1  =  C  A8S  (  YO  (  I  Ml.  00-2)  PHA  3020 

1400  CONTINUE  PFA3030 

IF  (SIGYOII). EO. 0.000)  S I G Y0( I  I » 1 .000  PFA3040 
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ooor-oor.  nooooooo 


IF  1 1 IGOI  JI.EQ.7I.ANC.IXI2,  1 1 .EQ.O.ODOI  I  X  1 2, 1)  -  25  8. 1  500 

T'XIl.II 

TC-XI2.il 

IF  IT0.NE.0.0I  GO  TO  1500 
WRITE  16,69001  I,7,Y0III,S1GY0I2) 

GO  TC  1600 
1500  CONTINUE 

WRITE  16,  70001  I,T,TO,  YOU  I  ,SIG»OI  II 
1600  CGNTINUE 
170G  CONTINUE 
laOO  CONTINUE 

WRITE  16,71001  I J , I KCUNTI J I  ,J-1,NSETSI 


C  THE  DATA  HAS  NOW  BEEN  STORED.  INITIALIZE  CONSTANTS  NEEDEO  FDR 
C  CRGLS2. 

NC-KCUNT 
NP=LISTP*7 
NX-2 
I  *> 0 


PHA3050 
PHA  3060 
PHA3070 
PWA3Q80 
PH  A  3090 
PHA3100 
PHA3110 
PFA3120 
PHA3130 
PHA  3150 
PHA  3 1 50 
PHA3160 
PHA31 70 
PHA  3 180 
PHA3190 
PHJ3200 
PHA  3210 
PHA3220 
PHA  32  30 
PHA  3250 


READ  A.NC  STORE, 

IFMIN  —  INPIT  FORMAT  CFR  PILISTP*7I 
0  INCICATFS  I6D12.5/D12.5I 
1  INCICATCS  (7A8I 

THE  LATTCP  FORMAT  HAS  NC  RCIJNCCFF  IN  TATA  TRAN'SFr’  . 
PILISTP*6I  —  TRIAL  PARAMETERS.  IF  VARIED,  “AY  NOT  BE  0.000 

IF  NOT  VAMRC,  0*  ICO  OR  A 


PHA3250 
PHA326G 
PHA  32  70 
PHA3280 
PHA  3290 
PHA33JO 
PHA  33 10 
PH  A  3  3  20 


PREDETERMINED  value  may  BE  USCDPHA-330 
KIILISTP.6I  -  0  IF  PARAMETER  IS  HELC  CONSTANT.  PHA-3A0 

I  IF  PARAMETLP  IS  TC  Be  VARIFC  PHA-35U 

IHOLC  —  INDICES  DF  PARAMETERS  IN  PIUSTP*7I  WHICH  ARF  PHA-360 

FORCED  TO  STAY  IN  REGRESS ICN  CESPITT  •TnST'  CETERMINA-  PFA3370 
TIDN  THAT  THeSE  PARAMETERS  MAY  Ee  NON-SIGNIFICANT.  PHA3380 


READ  15,51001  ( IFMINII|,I,1,LISTPI 
DD  2100  I-l.LISTP 
JO-1* 7*11-11 
JN*JC*6 

IF  IIFMIMII.Eq.il  GO  TO  1900 
READ  15,95001  (PI Jl , J-JD, JNI 
GD  TD  200C 
1900  CONTINUE 

READ  I5.55G0I  IPI JI.J.JO, JNI 
2000  CONTINUE 

DD  2100  J-JO.JN 
AAIJI-PI J| 

2100  CONTINUE 

READ  I5.51G0I  IK  1 1 1 1  ,  I »1,NP) 

IF  (NHJLD.NE.OI  READ  15,57001  I  I  HOLD  1 1 1  ,  I  =  I ,  NHCLOI 


C  BYPASS  RFCRESS ICN  IF  ICC  IS  3. 


IF  I  IDD.EC.3I  GO  TO  2600 
2200  CONTINUE 
NV-0 

CD  2300  K-l.NP 
IF  (KIIKI.EO.OI  GC  TO  2300 
NV-NV*1 
2300  CONTINUE 
JN-0 


■PHA3390 
PHA3500 
PH  A  25 10 
PHA3520 
PHA  35  30 
PHA3550 
PHA  39  50 
PHA3560 
PHA  3570 
PH  A  35  80 
PHA3590 
PHA  3500 
PHA3510 
PHA3520 
P.-A3530 
PHA3550 
PHA3550 
PHA3560 
PHA3570 
PHA3580 
PHA3590 
PHA3600 
PHA3610 
PHA3620 
PHA3630 
PHA3650 
PHA2650 


1 
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OOOOOOOOOOOO 


ENT;R  3RGLS2  AND  REFINE  PARAMETERS 


CALL  DRGLS2  (EAFW20) 

IF  1  ISING.NC-.OI  GC  TO  2600 


PUNCH  REFINED  PARAMETERS  ON  CARO.  CNC  PHASE  AT  A  T I  ME  AND  INCLUDE 
THE  PHASE  NAME  ON  THc  CARO  WITH  THE  RELATED  PARAMETERS. 

CO  2  5  30  l  3 1 ,  L I STP 
JC31»7*(  I- 1 1 

IF  (  IFMDUT.60.il  GC  TC  2400 
JN=JCh5 

PUNCH  820C.  (P(JI»J*JC,JN|, PNAM6 (II 
J  N  3  J  N  ,  l 

PUNCH  8  30C ,  P( JN  I  .DATE (PNAMEl 1 1 
GC  TO  25CC 
2400  CONTINUE 
JN=JC»6 

PUNCH  840C.  (  PI  J  I  I J- JC.  JNI  .DATE  »  PNAMEl  I  I 
2500  CONTINUE 

IF  I  ID0.cC.2l  GO  TC  4000 
2600  CONTINUE 

IF  ( IICY.CE.il  GO  TO  2700  _ 

SET  LP  ANC  PRINT  ERRCR  PLOTS. 

READ  IN  ANC  STORE  , 

NUL  —  NUMbER  OF  HORIZONTAL  DIVISIONS  CN  OPCINATE  LESS  I. 
NiaH  —  NUMEEF.  OF  HORIZONTAL  L I  .tS  PER  DIVISION. 

NVL  —  NUMBER  OF  VERTICAL  DIVISIONS  CN  ABSCISSA  LESS  1. 

Nsav  —  NUMBER  OF  VLRTICAL  LINES  PER  DIVISION. 

XMAX  MAXIMUM  tlrperature  ON  ABSCISSA. 

XMIA  —  MINIMUM  TCRPERATU.’C  ON  ABSCISSA. 
yMAA  —  MAXIMUM  ERROP  CN  ORDINATE. 

YMIN  *  -  YMAX.  _  _ 

READ  (5.45001  NHL  .NSBH , NVL ,NSB V 
RE  AO  (5.5  2001  XMAX.XMIN.YMAX 
YMIN=-YMAX 

2700  CONTINUE  .  . 

IF  ( 1  IS ING.NE.OI .AND.! I0C.NE.3I I  GO  TO  42C0 

IN*0 

CC  3303  J 3 1, NS6TS 
WRITE  (6,53001 

WR  I  T  T  (6,57001  (TITLEd  I  ,1*1,201  »OATE 


WRITE  (  6,  72001  TYFSdGOUll 

c 

c 

INITIALIZE  THE  PLCT  IMAGE. 

c 

CALL  PLDTI  (NSCALE. NHL, NSBH, NVL, NS8VI 
CALL  PLOT  2  (  IMAGE  ,  XMAX.XMIN.YMAX,  YMIM 

c 

c 

INITIALIZE  CONSTANTS  ANC  COUNTERS. 

c 

2  800 

IC  =  INU 

I N= I KOUNT ( J 1 

CC  2800  I 3 1 , 3 

SERR  (I  l  =  ZERO 

SERRSQ1 I 1 *ZERO 

CONTINUE 

c 

PHA3660 
PHA3670 
PHA3680 
PHA3690 
PHA3700 
•PHA37  10 
PHA3720 
PHA3730 
■PHA3740 
PHA3750 
PHA3760 
PHA3770 
PHA3780 
PHA  3790 
PHA3800 
PHA3810 
PHA  3820 
PHA  38  30 
PHA3840 
PHA3850 
PHA3860 
PHA  38  70 
PHA3B80 
PHA3890 
-PHA  3900 
PHA3910 
PHA3920 
P  HA  3930 
PHA3940 
PHA3950 
PHA3960 
PHA3970 
PHA  3980 
PHA3990 
PHA4000 
-PHA4010 
PHA4020 
PHA4030 
PHA4040 
PHA4050 
PHA4060 
PHA4070 
PHA4080 
PHA4090 
PHA4100 
PHA4110 
— PHA4120 
PHA4130 
— PHA4140 
PHA4150 
PHA4160 
— PHA4170 
PHA4180 
— PHA4190 
PHA4200 
PHA4210 
PHA4220 
PHA4230 
PHA4240 
PHA4250 
— PHA4260 
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WH«W!» W^.l .  'A .  > -  '«J( W I . ■  M,1  'A  -.'.  i A  ^  V« ,' 


»  3*’<-  .'  J  .  -!■*•+.»  L'Jfl  J  l.y.NJ 


C  CALCULATE  ANC  PRINT  YC.  CALCULATE  AND  SLM  The  CIFFFRENCE  (  FR  P  I  »  PHAA270 

C  PERCENT  CRRCP  IRELEkRI,  ANU  The  WEIGHTED  CIFFCRcNCE  (F|.  PHAA280 

C  - PHAA290 

□C  310J  1*10, IN  PHAA300 

CALL  EAFW20  (YC.II  PMA310 

FRRJP=YU( II-YC  PHAA320 

ERR=ERROP  PhA  A3  30 

IF  <YJIII.EQ.7ERC>  GO  TO  2900  PhAABAO 

ReLcRR  =  cRRCP*l.0C2/0ABS(Y0( 1 1 1  PhA A  350 

GO  TO  300C  PHAA360 

2900  CONTINUE  PHA A3  70 

PfcLEFR=ZERC  PHAA380 

3000  CCNTINUE  PHAA390 

F< 1 1 =eRR/ SIGYO ( I  I  PhAAAOO 

XI ( 1 1  * X ( 1 »I I  ph AAA  10 

TC*  X  (  2,  I  |  P  h  A  A  A  20 

WRITE  <6,73001  X I < 11 ,TG, YO < 1 1 , YC , fcPR  ,RFLEFk ,  H  1 »  PhA AA30 

C  - PhfAAAO 

C  ENTLR  THE  POINT  CN  ThE  PLOT  IMAGE.  FHAAA5C 

C  SUM  THE  ERRORS.  PhAAA60 

C  - PhAAA  70 

CALL  PLOT  3  <I1CC,XI,E,1I  PHAAA80 

S  =  RRI1I*SFRR(1MERR  Ph»AA90 

S  ERR  1 2 1  =  S  ERR  1 2  >  *REL6RR  PhAAbOO 

ScRnl3l*StRR(3l*Elll  PhAASIO 

SeRRSOUl  *S.RP.SJI  1 1  ♦  ERR*ER.P  PhAA520 

SERkSJ(2 l*SCRRS0<2l»RELCRR*RCLERR  PHAA530 

S:°RS9I  II  *S-PRS  )<3ME<  11*1  <  II  PhAAbAO 

C,  - PHAA550 

C  CALCULATE  ANO  WRITE  THE  STANCAPQ  ERROR  CF  ESTIMATE  FOR  THE  PHAA560 

C  DIFFERENCE,  ThE  PERCENT  ERROR,  ANC  THE  WEIGHTLD  01 F  FH»SNCE .  PHAA570 

c  - PhAA580 

3100  CCNTINUE  PHAA590 

NO  AT  A* 1*( IN- 10)  PHAA600 

AN  =NCAT  A  °h  A  A6 1 0 

00  3200  1*1,3  PHAA620 

CBAK<II*SERRU  I/AN  PhA A6 30 

STDe  VI I  I  *  SORT  <  SERRSC  <  I  I /AN-EUAR < 1 1 *EEAR (II)  PHAA6A0 

3200  CCNTINUE  PhAA650 

WRITE  (6,7/UUi  P  h  A  A  66  0 

WRITE  (6  7  ,001  ( E  EAR  < 11,1*1 ,3),(ST0EV1 I  7,  1*1, 3), NO AT A  PhAA670 

C  - PhAA680 

C  CALL  PUTOI.T  TO  CALCULATE  ANO  WRITE  THE  COMPLETE  SET  OF  CONSTANTS  PHAA690 

C  FOR  TH-  PHASES  ANO  FCR  THE  LATA  SET.  EYPASS  THE  LATTER  SET  PHAA700 

C  OF  CONSTANTS  IF  IGOUI  IS  GREATER  THAN  OP  FOUAL  TO  7  OR  IF  Ph AA 7 1 0 

C  THE  J-TH  OATA  SET  HAS  ONLY  ONE  PHASE.  PHAA720 

C  - - - Ph  A  A7  30 

CALL  PUTULT  (J)  PHAA7A0 

C  - -  44750 

C  PRINT  THE  PLOT  IMAGE  FOR  THfc  J-TH  DATA  SET.  PHAA760 

C  - PHAA770 

WRITE  (6,53001  Ph  AA  780 

WRITE  (6,5700)  <  T I TL  E  (I) ,  I  *  1 ,20 1  ,OATE  PHAA790 

WRITE  16,75001  < REF < I , J I , I  *  1 , 10 1  PHAA800 

CALL  PLOT  A  < NL , LABEL!  PHAA810 

WRITE  (6,7600)  PHAA820 

3300  CONTINUE  PHAA830 

C  - - - PHAA8A0 

C  PUT  CUT  TABLES  FOR  EACH  SPECIES  IN  THE  LIST  PNAME ( L I STP ) •  PHAA850 

C  - - - PHAA860 

SC*1 .000  PHAA870 
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ooo  00000 


00  3S0J  L*1,LISTP 
WR  ITE  (  6 ,  ‘4  00)  OATS 
WRITE  (6,55001  ( T ITL £ ( I » , I  *  1 , 20  I 
WRITS  (6,77001  PNAME(L) 

WRITE  ( 6 , 78001  ( TYPE ( I t , I • 1 , 7  I 
I  NOE  X* 1  *  7* ( L-l ) 

KOLNT* 19 

IF  (  IONIL  I.EQ.-l)  KCLNT*15 
M*(ICN(LI«4l/2 
CO  J83J  K*  l,  KCUNT 
T*TK(K,M) 

00  3600  0*1,7 
VAL(  JI*O.COO 
ICOES*J 

if  (  rm  l  i.rc.-i  i  igcfs=igofs*7 

JF  (  IGOES.  SC.  141  GO  TO  3600 
CA'.L  YOlRIV  (TK(K,M  ,  IGCES  ) 

00  J4JJ  1*1,7 

VAL(  J)*VAL(J>*OYCCm*AA(INOfcX*I-tl 

3400  CONTINUE 

IF  ( J.Nc.71  GO  TO  J60C 
CALL  YOERIV  (29fl.15C0.JI 
OC  3  50 J  1*1,7 

VAU  JI=VAL( JI-OYCCI I  l*AA(  IN  CEX ♦ I-  1 1 
35J0  CONTINUE 
3600  CONTINUE 

CO  3  703  1*1, 7 
AVAL)  I  l=VAL(  I  I 
3700  CONTINUE 

WRITE  (0,7930)  T,  ( A  VAL  ( 1 1 ,  I  *  1 , 7  I 
3000  CONTINUE 
JO*I  NOEX 
JN  *  J  C  *6 

WRITE  (6, {7001 

WRIT:  (6,80331  (AA(JI,J*JO,JNI 
PUNCH  050C,  ( AA( J ) , J* JO , JN I , PNAME(L I ,OATE 
3900  CONTINUE 

IF  ( ICO. EC. 31  GO  TO  4200 


IF  STEP-BACKWARO  ELIMINATION  OF  NON-SI  CM  FI  CAN  T  PARAMETERS  IS 
IN  EFFECT  (ICY.NE.UI,  SET  LEAST  SIGNIFICANT  PARAMETER  TO 
ZERO  ANC  REPEAT  REGRESSION. 


PHA4880 
PHA4890 
PHA4900 
PHA4910 
PHA492C 
PH A 4930 
PHA4940 
PHA4950 
PHA4960 
PH  A  49  70 
PHA4980 
PHA4990 
PHA5000 
PHA5010 
PHA5020 
PH  A  50 30 
PHA5040 
PH  A  5 J  50 
PHA5060 
PHA5070 
PHA5080 
PHA  5090 
PHA5100 
PHA5110 
PHA  5120 
PHA5130 
PH  A  5 1 4  0 
PHA5 150 
PHA5160 
PHA  5170 
PHA5180 
PHA5190 
PHA  5200 
PH  A  52 1 0 
PHA5220 
PHA5230 
PHA5240 
PHA5250 
PHA5260 
PHA5270 
PHA5280 
PHA5290 


4000  IF  (  ( IREDC.CQ.OI .CR. ( I ICY.GT. ICY  I  I  GO  TO  4200 
Kill RcOO 1*0 
P( IREOOI =C.OCO 
1 1 C  Y*  1 1  C  Y  *1 
GO  TC  220C 
4100  CONTINUE 

WRITE  (6, £1001  J 
4200  CONTINUE 
4300  CONTINUE 
STOP 


FCPMAT  STATEMENTS 


PHA5300 
PHA5310 
PHA5320 
PHA  5330 
PHA  5340 
PHA  53  50 
P  NA5360 
PHw  53  70 
PHA  53  80 
PHA  5390 
PHA5400 
PHA5410 
PHA  54  20 


4400  FCRMAT  (1CA8I 
4500  FORMAT  ( 1(15) 

4600  FORMAT  (2CA4I 
4  700  FORMAT  ( 4CI2) 

4800  FORMAT  ( A  8, 2X , 010 .3, 10X, 21  5  I 


PHA5430 
PHA5440 
PHA  54  50 
PHA  5460 
PH  A  54  70 
PHA5480 
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4903  FORMAT  (6C12.5I 
5000  FORMAT  (  15,13, 9A8/(  1GA8II 
51J0  FORMAT  (  OC  11  I 
5200  FORMAT  (8110.31 

5400  FORMAT  '|42hI  THESE  RESULTS  WERE  3ET  A I  NEC  IN  A  RUN  ON  ,A8I 
5500  FORMAT  (  IHO,  20A4/  IHO 1 _  ..  r._„,cc,rkl  luu 


PH65490 

PHA5500 

PHA5510 

PHA5520 

PFA5530 

PFA5540 

PHA5550 


5600  FORMAT  (  6CH0  °  PHASFS  CONSIDERED  IN  THIS  REGRESSION  ARE  AS  FCLLCPFA5560 


>  ^  J  - 

IWS-“"“"/IHOflOX*5A2C/(lH  ,10X,5A20II 
5700  FORMAT  (IF  i20A4*40XfA8/lHGI 


PHA5570 

PFA5580 


lilt  5°0::S?  ( 54F0THE  FClLCWlNG  DATA  SETS  HAVE  BEEN  READ  IN  TO  STQRAGE:/PHa5590 


>uw  4  w  ■»'<  i  *  - -  P  h  A  S6  0  0 

1 1 1  HO  SET  NLM8ER  ,1 C X i 9HPcFt RENCE  I  PHA5610 

5900  FORMAT  ( 1F0I  PHA5620 

6000  FORMAT  (1H  ,A8)  PHA5630 

6200  FORMAT  !  2 'HOTHF ‘nUMBER^QF  PHASES  lS«,lo/3lH  THE  NU-8ER  ^SFRVAT  PF  A564  0 
20  IKNS  Is:'l5/21H  THE  TYPE  OF  DATA  I S  : .  3X  .  A  E  /  1  1H0FH  A  S£  N  AMF  .  1  CX,  UHC  PF  650 
20FFF IC 1  ENT  »9X* 10FPEF .  STATE , 1  OX , 1  OH  1  NVC RS ICN S7 1 H  I  .  7Q 

6300  FORMAT  (IF  ,  IX  ,  A  8 . 1 1  X ,  F  8.  3  i  I5X ,  A  3 . 1  4X  ,  1 1 ) 

04 JO  FORMAT  (lFO,lOX,2tHlNV-»SIONS  ARE  AS  FCLLCwS:/lH  I  '  „Q 

6503  FORMAT  (1H  ,  F 1 5. 2 ,3 X *A a »5X, IH= , 5X , A 8 ,F20.  31  _  _ 

6603  FORMAT  (1F3,10X,A8,5X,58HIS  A  REFERENCE  PHASE  AND  HAS  ZERO  GI8BS  E^^703 


1NCRGY  AT  298.15  K/1H0I 
6700  FORMAT  (IF  1 


PFA57IO 

PFA5720 


6  800  FORMAT  i  \l  !  1 5 X,5HNCM.2 IX ,  1 1HTEHPER  ATUR E.  1 8X ,  A  8 ,  1 3X,  5HFR> C>  /  IH  .PM5730 


1 3 bX • 2HT 1 1  1 8X v 2H t  1 / 1M  )  ci  PHA5750 

69j3  FCRMAT  (IF  ,120.F2(j.3.16X11H-,3X,1P2C20.5I  a'760 

HOG  FORMAT  ( 9 ’HI  THE  *  INC  f  X^OF^TH^LAST  ITEM  OF  THE  ABCVE  DATA  SETS  IN  TPFA57  70 
1HE  V  r  C  T  0  R  S  X  ( I  *  I  I  «  YCl(  II*  ANO  SIGY01II/15F  IS  AS  F  OLLOV.S  :  /  IHO ,  30X  ,  PHA5  7  80 
28HCATA  SET.3X.5H1NCSX//1  IH  ,24X,2I10II 

72oo,:°^a; 

24HT(’ll  ,9X,BHOBSERVrO,8x,10HCALCULATEC,6X,10HtIFFERENCE,4X,7HPERCEN"FA5820 


3T«2X«14HCHS-CALC/SIGY0I  p^5B40 

?400  forma!  (55H0THE2ApiTF2ETlC  MEAN  CF  THE  EFRCRS  FCR  THIS  SET  ARE  --PFA5850 

1  - - 3F 1 3 .4/ 59H0T  FF  STANDARD  DEVIATIONS  ABUIT  THc  MEAN  ARc -  PFA5860 

2  - -  3F13.4/31H0THE  NUMBER  OF  OBSERVATIONS  IS  ,151  PM5B80 

7500  FORMAT  (IFO.lOAd/lH  1  p(.,cag0 

7600  FORMAT  (  1  F  0 ,  T  3  5 , 2  1HT  F  R  PFRA  T  UPE  (KELVINSII 

7700  FORMAT  (IFO/IH  ,  A  E/1H0I  .  dhabqio 

7800  FCRMAT  (IF  ,  1  5X  ,  A  E  ,  8  ) ,  A  8,  a  X  ,  A  8,  8  X ,  A  8 , 8X  ,  A  e ,  E  X ,  A  e  ,  8  X  ,  A  8/ 1 H  I  PL.JS970 

ssss  F,gs:s? 

BIOO^CRMAT1!  -ISH^YOUR^IDWILY  COMPUTER,  OC  NCR  CUIT./9H01N  YCLP 

146FTF  DATA  SET  YOt  HAVE  MISSTYPEO  TFE  PHASE  NAME. I  p^5970 

8200  FCRMAT  (  l  F6C 12 . 5 , A  8 )  PHA59B0 

8300  FORMAT  (  I P ID  12. 5 , 52X ,2A 81  PHA5990 

8400  FORMAT  (7A8,8X,2Afil  phafqoo 

8500  FORMAT  ( 7 A8 . 8HC0R FLETE , 2A8  1  PHA6010- 

ENO 
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l 


c 

c 

c 

c 

c 

c 

c 

c 


fllCCK  DAT*  0LK 

IMPLICIT  RSAL*8<A-H,C-ZI  BLK 

L0GICU.*1  UCC, LABEL  8LK 

01  ML  NS  UN  SCINVt2I,STCCEFt 10I,YCSN0C2»,TKI19,2I,TYPEI  16 » ,NSC «LC < 5  I PLK 
1  tLABCL ( 50)  BLK 


10 

20 

30 

AO 

50 


—  ELK 

60 

COMMON  BLOCKS 

BLK 

70 

NAME 

ROUT INES 

0 1  K 

80 

EARTH 

MAIN  PROGRAM  , 

CRCLS2, 

c A  Fw20 , 

putout 

PLK 

90 

AIR 

MAIN  PROGRAM, 

OR  CL  S  2 , 

TCST20, 

EAFW20 

BLK 

100 

PUTLUT  ANC  UNICLE. 

eLK 

no 

FIRE 

MAIN  PROGRAM, 

7RCLS2, 

TEST29, 

EAFW20 

BLK 

120 

RATER 

MAIN  PRUG5  AM  f 

EAFW23, 

YCf.RIV, 

ODER i V , 

BLK 

130 

PUTOUT,  r-i,  01 

=  ,  AND 

PLCCK  OATA. 

PLK 

1  AO 

7  !  M  c 

MAIN  PROGRAM, 

PLTCUT 

PLK 

153 

SPACE 

MAIN  PROGRAM, 

-AFVJ20, 

YCEC.IV, 

ODER  IV, 

BLK 

160 

AND  FUTFLT 

eLK 

170 

PAN 

MAIN  PROGRAM, 

PLTCLT 

olk 

—eLK 

180 

190 

COMMON  /JCULES/UMTS 

CCWMCN  /NATFR/ZFRC,  ONE,  TrfC,  THREE,  FOUR,  SIX,  R,  F, 

1  SCINV,  TRL  F ,  STCOjF,  OlfcO,  A,  8,  TEETA,  YESHC,  TK, 

2  ASTAR,  TYPE,  NL  ,  N  SC  ALE ,  LABEL,  BCC 


P  LK 
“LK 


2  JO 
210 


THE  FOLLOWING  VALLES  CF  »R*  ANC  *F» 
UAT»  IN  CALORIES. 

CATa  R,f/6. 5 75616 50 0,2. 3 06 056606/ 

AP.- 

TC 

USnC 

FOR 

ENERGY 

THE  FOLLOWING  VALUES  CF  'R>  AND  »F* 
DATA  IN  JOULES. 

CATA  R,F/ 19.1663  7500,5.66870006/ 

ARE 

TC 

BE 

USSC 

FOR 

ENERGY 

BLK  220 
PLK  230 
BLK  2 AO 
PLK  250 
eLK  260 
PLK  270 
PLK  280 


CATA  L'IIT5/UF  JOLLFS  / 

OATA  R,F/  19. 166 3 7500,9.66  8  70006/ 

CATA  ZERO ,CNE,T«C, THREE, FOUR, SI  X/O.OCO,  1 . 00 C ,2 . OOC , 3 . 003 , A . 000 , 

I  fc.COO/ 

CATA  SCINV/-1.0D0.1.CCO/ 

CATA  TREF/2. 981502/ 

OATA  STCCEF  /I  .OCO,  2 . OCO , 3 . 000, A .  0C0.5 . OCO , 0. 500, 0. AOO, 0 . 3 2 3333 3 
133C0,0.25C0,C.2  10/ 

OATA  DItC  ,A,8,THET A/3. 05702,  1. 8 75C-2 1 . 2 7A  101,2.1902/ 

CATA  YF.SNr  /  3HNG  ,3HYF5/ 

CATA  TYPE  /  8HF  S  AT  CAP,  8KNTR0PY  ,  8H  ENTHALPY,  8HGI  8BS  EN, 

1  eHUlG  K  ,  8HC  ELL  EMF , H HFT2-HT 1  ,7*8HSPECIAL  / 

OATA  ASTAR/8H********/ 

OATA  TK  / 273. 1500,298. 15J0,  323 . 1 500, 3A3 .  1500 , 3 7 3.  1 500 , 398.  1 5C 0 , 

1  A23. 150  0,666.  1500, A  73. 15D0.A98. 15D0,523. I5D3,5AR. 15D0.573. 1500, 

2  5S3.1500,623.15CC,C. COO, 0.000, 0.030,0.000, 2 72. I 5CC , 2 98 . 15C 0 , 

3  A.  J02, 5.CC2.6.0C2,  7 . 002 , 8 . 00'. ,  9.002 , 1.  OC  2  ,  1 .  1C  3,  l .  203 , 1 . 303  , 

A  1. AC  3, 1.503,  1.60  3, 1.703,1. 803,1.90  2,2.00  2/ 

CATA  LABEL  /30HOeS  LESS  CALC  OVLR  UNCERTAINTY/ 

CATA  NL/3C/ 

OATA  BCO/lt-O/ 

CATA  NSCALE/5*0/ 

ENC 


PLK 

PLK 

eLK 

BLK 

PLK 

BLK 

BLK 


290 

300 

310 

320 

330 

3A0 

350 


BLK  360 
BLK  370 
IBLK  380 
BLK  390 
BLK  AOO 
BLK  A10 
BLK  A20 
BLK  A 30 
8LK  AAO 
BLK  A50 
BLK  AtO 
PLK  A70 
BLK  ABO 
BLK  A90 
BLK  500 
BLK  510 
BLK  520 
BLK  530 
BLK  560- 
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SUBROUTINE  0RGLS2  (CALC  I 


A  GENERAL  FCRTRAN  LEAST  SQUARES  PROGRAM 
3Y 

mILLIAM  R.  BUS  I NG  ANC  FENR I  A.  LEVY 
CHEMISTRY  DIVISICN 
DAK  RIDGE  NATIONAL  LAECRATCRY 
CAK  RIDGE,  TENNESSEE 

REFERENCE —  BUSING,  R.  R.  AND  LEVY,  F.  A.,  1962,  CR  GLS,  A  GENFR 
FCRTRAN  LEAST  SCLARES  PRCGRAM,  CAK  RIDGE 
NATIONAL  LABORATORY,  CAK  RIDGE,  TENN., 
CPNL-TM-271. 


SU8RCUTINE  CRGLS2 

DRGLS2  HAS  BECN  ACAPTtC  FROM  BUSING  ANC  LEVY'S  UR  GLS 
PROGRAM  TC  ALLDw  LSc  GIF  THE  STATISTICAL  ROUTINES  KITHCUT  BEING 
REST"  ICT  EC  TC  THE  RIGID  INPUT-OUTPUT  FEATLRF.  S  OF  TH.-.  PPOGRA". 
THE  CCNVFRSJCN  HAS  PERMITTED  CHANGES  KHICF  CLARIFY  THE  OUTPUT 
FDR  THIS  PECULIAR  PROBLEM. 

ADAPTATION  BY  HAAS,  B/20/72. 


IMPLICIT  REAL*8  (A-H.D-Z) 

REAL*'.  TlTLf  .PCTERR.CIT 
OIM;;NSIUN  OUT  (  71 

DIMENSION  SQSIGC 2), AM< 10013). V( 143 l,CV( 14C) .DIAGI 140) ,RCW( 1401 
DIMENSION  CCSFC6,73),PNAMfc(20), TINV<6,70,4> ,  I PHA SC- ( 6, 70 1 , 

1  NPHAiE ( 70) ,  IKOLNTI 7C  I, l GDI  701 , I  STATE! 6, 70 1 .NINVER 1 6 , 701 , 

2  INSTAT(6,70>  ,  INV»H C6 , 7 J, 5 1 , 1 NVSCI 6 , 70) 

DIMENSION  X( 2,  1200) , YC(  12 0 0 ) , S I GYC 1 1200 1  ,P(  14  0  I , KI 1 1 40  I , 

1  CC 1 14  0  )  ,PD( 140),TITLEI20I 
DIMENSION  REF< 10, 70) ,CRRP< 140 ) , I HCL C (40 ) 


CCMMCN  CLOCKS 
NAME 


ROUTINES 


EARTH 

MAIN 

PROGRAM, 

DRCLS2, 

EAFW20. 

PUTDUT 

MR 

MAIN 

PROGRAM, 

DPCLS2, 

TCST20, 

EAFH20 

PUTCUT  AND  UNIQUE. 

FIRE 

MA  IN 

PRCGRAM, 

DRCLS2, 

TEST20, 

EAFR20 

MATER 

MAIN 

PROGRAM, 

EAFH20, 

YCl R  IV, 

ODER  I V 

PUT'JUT,  FN ,  DIE,  AND 

PLCCK  DATA. 

TIME 

MAIN 

PROGRAM, 

PUTOUT 

SPACE 

MAIN 

PRCGRAM, 

EAFW20, 

YCFRIV, 

ODER  IV 

AND  PUTPUT 

MAN 

MAIN 

PRCGRAM, 

PUTCUT 

COMMON  /EARTH/  CCEF , PNAME.T  1NV , I PHA SE , NPE AS E , I KCLNT ,  I  GO ,N SETS , 

I  ISTATc, NINVER, INST  AT, I NV PH, I NV SC , L I  ST P 
COMMON  /AIR/  X,  P,  DC,  TITLE,  Y3,  SIGYD,  PD,  KI,  NC,  NV,  NX,  IW 
1  NP,  NO,  ISING,  ISTCP,  IL,  JDFLAG 
COMMCN  /FIRE/  REF,  RPRP,  IHRITE,  ICY,  IICY,  IREDC,  IRFG,  NHCLD , 
1  IEOLO 


CRG  10 
— ORG  20 
CRG  30 
C«G  40 
CRG  50 
CRC-  60 
rPG  70 
CRG  80 
CRC  SO 
DRG  100 
ORG  110 
CRG  120 
ALCRC  130 
CPG  140 
C»C  150 
DPG  160 
CRG  170 
CRG  180 
CRG  ISC 
CRG  200 
CRG  210 
CRG  220 
CRG  230 
CRG  240 
CPG  250 
CRG  260 
DRG  27C 
CRG  280 
CRG  2S0 
—  CRG  300 
CRG  310 
CRG  320 
CRG  330 
ORG  340 
DRG  350 
CRG  360 
CRG  370 
DRC-  380 
CRG  3S0 
ORG  400 
—CRG  410 
CRG  420 
CRC  430 
CRG  440 
CRG  450 
CRG  460 
CRG  470 
DRG  480 
CRG  4S0 
DRG  500 
CRG  510 
CRG  520 
CRG  530 

- ORG  540 

CRG  550 
DRG  560 
,  CRG  570 
ORG  580 
CRG  590 
CRG  600 
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C  FORMAT  STATEMENTS  ORG  610 

1J0  FORMAT  (1F120AAI  nRG  620 

200  FORMAT  (32H0NUMeER  OF  CYCLES  IN  THIS  JOB  I S I  2/ 3 7HQNUM8FR  OF  PARAMcORG  630 
ITERS  TO  BE  VARIEC  I SI3/51HQNUMBER  CF  INCEFENCENT  VARIABLES  PER  CESCRG  6A0 
2ERVAT ION  IS  12  >  CRG  650 

3J0  F CRM  AT  ( A  6 HOCFR I V  AT  I VE  S  PROGSAHMEC  IN  SU8RCUTINE  EAFW20.  1  CRG  660 

ACO  FORMAT  (  31H0Wt:lGHTS  TC  BE  SU.'PLIEO  eY  USED)  ORG  670 

50J  rCRMAT  1 3  AF3UN IT  WEIGHTS  TO  BE  SET  EY  PRCCRAM)  CRG  68C 

6J3  FORMAT  ( 2  SHONUMBCR  OF  PARAMETERS  REAO  ISIA)  ORG  690 

73J  FORMAT  ( 3 1H0NUM 3E R  OF  OBSERVATIONS  READ  ISI51  CRC  700 

800  FORMAT  (A  AHOCALCULATEO  Y  8AS70  ON  PARAMETERS  BEFORE  CYCLEI2)  CRG  7  10 

900  FORMAT  (  1F0, 10X,5FINCEX,  7X,AHT<  2  1 , 8X.  AHT  (  1 )  ,  10X  ,£HY  (08S)  .  10X,7HY(CDOG  720 
1ALCI  ,8X,  8PC8S-CALC,5X,6HSIG(0I,6X,  12Hi' C-C  )/SIG(C>  )  CRG  730 

100.1  FORMAT  (IF  ,U5,F12.3,12X,1PA016.3,0PF16.A)  CRG  7A0 

1100  FORMAT  (IF  » 1 1 5 , 2F 12 .3, IPAQ  16. 3, OPF 16. A  )  TRG  750 

1200  FCOMAT  (IF  ,I15,F12.3,12X,5F16.A)  CRG  760 

1300  FORMAT  (IF  ,I15,2F12.3,5F16.AI  CRG  770 

1A00  FORMAT  (  5 lHOAGRGEMtNT  FACTORS  BASEO  CN  PARAMETERS  BEFORE  CYCLE  1 2/2CRG  780 
lOHOSUMIrf* (O-CI **2 1  IS  C11.3/35H0SG»TF!SUM|W*(O-C 1**2) /( NP-NV1  I  IS  CRG  790 
2C10.AI  QPG  0 Q Q 

1500  FORMAT  ( 6CF0CSTIMATE0  AGRCFHENT  FACTORS  E  AS  EO  CN  PARAMETERS  AFTFR  CPC  810 
1CYCLCI2/2CH0SUM(W*(0-CI**2I  IS  0 11.  3/ 35HCSCR TF ( SUM  I w* (0-C )**2 I /(  NOCPG  820 
2-NVI)  IS  CIO. A)  ORG  830 

1600  FORMAT  ( 6  2H  MATRIX  HAS  A  ZERO  0 1  AGONAL  ELEMENT  CORRESPONDING  TC  PACRG  8A0 
l»AMtTERI3,16H  OF  TFCSF  VARIFOI  nRG  850 

1700  FORMAT  (ACH  SINGULARITY  RETURN  FROM  MATRIX  INVERTER)  CRC-  860 

1800  FORMAT  (37HOPARAMETCRS  AFTER  LF/ST  SOLA  RES  C YC LC I  2/ 1H0  ,T 16, 'OLD* , TORG  870 
1  35  ,  'C'-'NNGF*  ,T57,  •  NEW*  ,T77,  ‘ERRCR  •  ,T95,  '  PCT.  CHANGE *,T1 18, 'PCT.  F.FPOPG  880 
2CR '/ 1H  I  CRG  890 

1900  FCRMAT  (1F0.A8)  CRG  900 

2000  FORMAT  (IF  I3,2(AX,1P017.10I)  CRG  910 

2100  FORMAT  (IF  ,  13,  IP  502  1 . 10, OPF 1 7.  A  I  CRG  920 

2200  FORMAT  ( 66HJSUBRCLT  IN'E  TEST  INDICATES  THAT  JC8  IS  TO  BC  TFRNINATEOORG  930 


1  FOR  RGASCNI 2  I 


CRG  9A0 


2300  FORMAT  (16H0TRIAL  CCNSTANTS/  1H0, 2X ,  1H 1 , 9X  ,AHP  ( II  ,6X  ,  5HK I  (  I )  /  1H  I  ORG  950 

2A00  FORMAT  (1F0.9X.A8I  0RG  960 

2500  FORMAT  (IF  ,  1 3 , 5 X  ,  IRC  1  l.A , 6X,  I  1 1  OaG  970 

2600  FORMAT  (  1 SHOCORR E IAT I CN  MATRIXI  CRG  980 

2700  FORMAT  ( 1 FOI 3 , 100  12. A/( 1H  3X, 10012. All  qrG  990 


2803  FORMAT  (  1 AH1  ***W  ARNI  NG*«*/9  3H0THL  TERM  SIG/INC-NVl  IS  NEGATIVE.  TCRG1000 


1HE  ABSOLUTE  VALUE  IS  TAKEN  ANO  THE  REGRESSION  CONTINUES.) 
CO  2900  !*1,NP 
ERRP  (11*0.000 
2900  CONTINUE 

WRITE  (6,1001  (TITLE(1I,I*1,20I 
WRITS  (6,200)  NC.NV.NX 
WRITE  (6,200) 

IF  < I W  »  3  100,3000,3100 
3000  WF.ITE  (6.A00I 
GO  TO  320C 
3100  WP.ITE  16,500) 

3200  CONTINUE 

WRITE  (6, £00)  NP 
WP.ITE  (6,700)  NC 
IF  (NC)  3200,3200,3500 
3  300  CO  3 AOO  I*1,NP 
3A00  K I  (  I  1  =  0 

C  INITIALIZE  PROBLEM  ANC  ENTER  SUBROUTINE  PRELIM  IF  PROVIDED 
3500  NM=(NV»(NVM)l/2 
SQSIGI 1 ) =C.OOO 
CALL  PRELIM 


ORG1010 
CRC 1020 
0RG1030 
CRG10A0 
CRGI050 
CRG  1060 
CRG  1070 
CRG 1080 
O' G  1090 
CRG1100 
0RGU10 
CRG  1120 
CRG  1 130 
0RG11A0 
0RG1150 
CRG  1 160 
CRG1170 
CRGU80 
CRG  1190 
DRG1200 
CRG  12 10 


PUT  OUT  TPIAL  PARAMETFRS,  AN0  K£ V- INTEGE R  S . 

ORG  1220 

WRITE  (6,2300) 

CRG1230 

J0«1 

CAG1240 

CC  36J0  I=1,L!STP 

CRG1250 

WRITE  (6,2400)  PNAME(I) 

CSC  1260 

JN  *  J  C*6 

CRG1270 

WRITE  (  6,2500)  (  J,P(  JI,KI(JI,J*JC,JM 

CRG1280 

J0=JC*7 

C°G 1290 

3600 

CONTINUE 

CRG1300 

START  LOUP  TO  PERFORM  NC  CYCLES  AN0  ONF  FINAL  CALCULATION  OP  Y 

CRC 13  IQ 

NCY*NC»1 

CRG1320 

CO  /500  I C  =  1 ,NC Y 

OR  G 1 330 

CLEAR  ARRAYS  AM  ANC  V  EXCfcPT  CN  LAST  CYCLE 

CRC1340 

IF  (IC-NCY)  3  700,4000,4000 

CRG1350 

3 /JO 

00  3630  1*1, NM 

0RG1360 

3000 

AM ( I )  =  3.0C0 

CR5137C 

00  35JJ  1*1, NV 

ORG  1380 

3900 

V( 1) = J.ODC 

CRC1350 

IMTIALIZF  FOR  CYCLu  IC  AND  PUT  OL  T  CAPTION  FCR  LIST  CF  Y  (CALC )  OG  1400 

4000 

SQ  S1GI2 I  =  5CS  IG( 1 ) 

C 1G 14  1C 

SIC  =  0.  JCO 

CRC 1420 

WRITE  (6,100)  (TITLFt II , 1=1 ,20) 

CRG1430 

IF  (  IL.NE.O)  GO  TC  4100 

CRC 1440 

IF  (  IC.NE.NCY)  GO  TO  4200 

CRG1450 

4100 

CONTINUE 

CRG1460 

WRITE  (6,E0J)  IC 

CRC 14  70 

WRITF  (6,«03) 

CRG14K0 

4200 

CONTINUE 

CR51450 

START  inoo  THROUGH  NO  OBSERVATIONS 

C5C1500 

C')  5530  1*1  , NO 

0RG1510 

ENTFR  USE0  S  SUBROUTINE  TO  CCVFUTE  Y(CALC)  ANC  DERIVATIVES 

GAG  1520 

CALL  CALC  (YC, I  ) 

0RGI530 

OBTAIN  WEIGHT  ANC  CALCULATE  CUANTITIES  FRCF  Y I  OB S )-Y ( C ALC ) 

ORG  1540 

IF  (IW)  4400,4300,4400 

CRC- 15  50 

4300 

SORTW’l.CCC/SIGYCtI) 

ORG  1560 

GO  TC  450C 

CRG1570 

4400 

SIGYOI  I  )  *  I. 000 

GRC1580 

SCRTW=1.0CO 

CRC,  1590 

4500 

OY  =  YC(  I  )- YC 

ORG 1600 

WDY*SCRTW«CY 

CRGI6I0 

SIG=SIG»WCY*WDY 

ORG 1620 

PUT  OUT  YICALCI  ANO  OTHER  INFORMATION  FOR  CNE  OBSERVATION 

CRC  1630 

IF  (IL.NE.O)  GO  TC  4600 

CRC 1640 

IF  (  IC.NE.NCY)  GO  TO  5100 

CPG1650 

4600 

CO  NT  I.NUE 

CRGI660 

OUT (  1) *X(  1,1) 

0RG1670 

0UT(2)*XI2,II 

CRGI680 

CUT ( 3 ) * YQ ( I ) 

CRG1690 

CUT(  4  I  *  YC 

0  'G  1700 

OUT( 5 ) =D Y 

CRGI7I0 

OUT ( 6) =S I CYO( I) 

0RGI720 

OUT(  71  =WDY 

ORGi/JO 

IF  ( X( 2, 1  1.EQ.C.0C01  GO  TO  4800 

CRG1740 

if  ( (DAUSIYCI I ) ) .LT.  1.00-3) .AND.(YO(  I) .NE.O.OOO)  1  GO  TO  4700 

0RGI750 

WRITF  (  6,  13001  I  ,  (OLT(K)  ,K*I,7) 

CRC 1 760 

GO  TO  500C 

CRG1770 

4  700 

CONTINUE 

CRG1780 

WRITE  (6,1100)  I , CUT (l),OUT(2l,YO(I)fYC,CY,SIGYC(I),OUT(7) 

CRC  1 790 

GO  TC  5J0C 

CRG1B00 

o 

o 

CO 

*r 

CONTINUE 

CRGI8I0 

IF  I (DABS(YO(II).LT. 1 .00-3 ) . AND. I YO(II.NE.O.OCOI)  GO  TO  4900 

0RGI820 
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WRITE  (6.  12001  I  ,CUT  (1 1  ,  (OUKKI  ,K=d,7  I  0RG1830 

GO  TC  500C  CRG1840 

CONTINUE  CRG1850 

WRITE  (6.10031  I.CUTdl.YOdl  ,  YC ,0Y , S IGYC (I  I  ,0UT(7I  0RG1860 

CONTINUE  CRG18  70 

CONTINUE  CRG1880 

BY-PASS  DERIVATIVE  ANO  MATRIX  SET-LP  ON  FINAL  CALC  OF  Y  0RG1890 
IF  (IC-NCYl  5200. 5400 . 5400  C5G1900 

START  LOCP  TC  STORE  AN  ARRAY  OF  NV  CERIVATIVES  0RG1910 

J*1  C»G1920 

CO  5400  K *  1 , NP  ORG 1930 

IF  (XKKII  5400,3400,5300  QRGI940 

CONTINUE  CRC1950 

OBTAIN  DERIVATIVE  FRCP  THCSE  PRCGRA HHEC  BY  USER  CRG1960 
CV( J I  *  SQR TW*OC ( K I  CRGI970 

J“  J  ♦  1  CRG1980 

CONTINUE  CRG1990 

ENO  LCCP  TO  OBTAIN  OFRIVATIVES  CRG2000 

START  LCCP  TC  STORE  MATRIX  ANO  VECTOR.  r9G2010 

1604  CR  CL  S  STORAGE  SCHEME  IS  REVERSE  OF  7090  CR  GLS  CRG2020 

JK*1  CRC2030 

CO  5  8 JO  J  *1 , NV  CRG2040 

TfcMP=Ov(JI  CRG2J50 

IF  (TEMPI  5600, 55CJ, 5600  CRG2060 

3Y-PASS  IF  DERIVATIVE  IS  ZERO  CRC2070 

JK  =  JKfrNV»  1-J  CRC2080 

GO  TO  580C  CRG2090 

CO  5  700  K»J,NV  0RG2I00 

AM(JK|>AH(JKI*TEMF*OV(KI  ORG2IIO 

JK=JK*1  CRG2I20 

CONTINUE  CRG2I30 

V(  JI=V(  JMTfcHPAWDY  CRG2140 

CONTINUE  0RG2150 

ENC  LCCP  TO  STORE  HATRIX  AND  VFCTCR  CRG2160 

CONTINUE  CRG2170 

ENO  LOOP  THROUGH  NO  CBSERVATICNS  0RG2180 

COMPUTE  ANO  PLT  OUT  AGREEMENT  FACTORS  C°G2190 

CNVN=NO-NV  CRG2200 

SQ  S I  G ( 1I  =  CSQRT(SIG/(CNVNI I  CRG2210 

WRITE  (6,1400  IC.SIG.SCSIGdl  CRG2220 

oY-PttSS  MATRIX  INVERSION  ANC  PARAMETER  CUTPUT  ON  FINAL  CYCLE  0RG2230 
IF  (IC-NCYl  6000, 80G0, 8000  CRG2240 

START  LOCP  TO  TEST  FOR  ZERO  DIAGONAL  ELEMENT  0RG2250 

ISIMG»0  CRG2260 

d  =  l  CRG2270 

1 1 Q*N V  CRG2280 

DO  6300  1*1, NV  0RC2290 

IF  (AM(IIII  6200,6100,6200  CRG2300 

ISING*1  ORG2310 

WRITE  (6,16001  I  CRG2320 

1 1  * I IfIIO  CRG2330 

IIC*II0-1  0RG2340 

CONTINUE  CRG2350 

ENO  LOCP  TC  TEST  FOR  ZERO  CIAGCNAL  ELEMENT  CRG2360 

TERMINATE  JCB  IF  ZERO  OIAGCNAL  ELE“ENT  WAS  FCUNO  CRG2370 

IF  (ISINGI  8600, 64G0,  8600  CRG2380 

FNTER  SLEROLTINfc  TO  REPLACE  MATRIX  WITH  INVERSE  0RG2390 

CALL  MINV2G  ( AM ,NV, I  SING  I  CRG2400 

IF  (ISINGI  6500,6600,6500  0RG2410 

TERMINATE  JOB  IF  SINGULAR  MATRIX  WAS  FOUND  CRG2420 

WRITF  (6,17031  0RG2420 


133 


o  o  or 


GO  T C  860C 

C  START  LOOP  FOR  MATRIX  VECTOR  MULTIPLICATION  FOR 

C  PARAMETER  CHANGES 

6600  00  7 10 J  I«1.NV 
PCI  =  0.000 
IJ«I 

IJOaNV-l 
CC  /COO  J  a  1  *  NV 
POIaPOI*AM(IJ)*V( J) 

IF  *  J- 1  I  67C0,680C,65G0 
6700  I J  =  I  J*I  JC 
IJC=IJO-l 
GO  TC  70CC 

C  SAVE  01  AGONAL  ELEMENTS  CF  INVERSE  MATRIX 

6800  DIAG(II*AM(IJ> 

64 JO  IJaIJ*l 
7000  CONTINUE 
P0< I l=P 01 
S IG=S IG-P  C I *V(  II 
7103  CONTINUE 

C  l:NC  LOOP  FOR  MATRIX  VECTOR  MULT  I  PL  I C  A  T  I CN 

C  RcCOMPUTF  AGREEMENT  FACTOR  USING  “OCIFIEO  SIG 


IF  (SIG) 

7200, )30C, 7300 

7200 

CONTINUE 

WRITE  (6 

,2800) 

SIGaCAdS 

(SIG) 

7300 

CON'T  INUE 

SCSI G(  11 

=CSCRT(SIC/(NC-NV)I 

C  PUT  QLT  CAPTION  FOR  LIST  CF  CCFRECTtC  PARAMETERS 

WRITE  (6.1C01  (TITLF(l), 1-1,201 
WRITE  (6,  1800)  IC 

C  START  LOCP  TO  CORRECT  ANC  PUT  CUT  PAPA-STEPS 

Jal 
J0  =  1 

CO  7700  I C-l »L  ISTP 
WRITE  (6,14001  PN  A  ME  (  1 0  I 
JN-JC»6 

CO  7600  I  *  JO  ,  JN 
IF  ( K I ( I )|  7400,7400,7500 
7403  WRITE  (6,20001  I .  F ( I  I  *  P ( 1 1 
GO  TC  760C 
7500  POLO  =P (  I  I 

P(I laPOLCNPC(JI 
SIGPaOSJF  TICIAGI J  ll»SCSIG(l I 
PCTCt-Gal.C2*P0(JI/P(I  I 
PCTE  PR  =0A  tS (  1.002«SIGP/P(I  I  I 
FRRP ( I ) a  PC  TERR 

WRITE  (6,21001  I,FOLO,PO(JI,P(U,SIGP,rCTCHG,PCTERR 
JaJ*l 

7600  CONTINUE 
J0»JC(7 
7700  CONTINUE 

END  LO.P  TO  CORRECT  ANO  PUT  CUT  PARAMETERS 
PUT  OUT  ESTIMATEO  AGREEMENT  FACTORS 
WRITE  (6,1500)  IC  ,  SIG  ,  SCSIG  ( 1 1 

ENTER  USERS  SU8R0UT I  NS  TO  TEST  ANO  MODIFY  PARAMETERS 

OR  ENO  JOB 

ISTOP«0 
CALL  TEST 20 

C  TERMINATE  J08  IF  INOICATEO  BY  USERS  SUBROUTINE  TEST 

IF  IISTOP)  7800,7500,7800 


CRG244Q 
0RC2450 
C»G2460 
0RG2470 
CRG2480 
CRG2450 
CPC2500 
CRG2510 
0RG252O 
C’C2520 
CRG2540 
ORf  2550 
C°G2560 
CRG257C 
0RG2580 
CRG2550 
CRG2600 
C°G2610 
C°G262u 
0RC263O 
CRC264G 
CRG2650 
0RC2660 
CRG2670 
OP G 26 80 
CRC2640 
C»G27CO 
r  1 C  2  7 1 C 
CPG2720 
CRG273C 
COC2/40 
CPG2750 
0RG2760 
CSC2770 
0®G2780 
0RG2750 
CRC2800 
CRG2810 
CRG2820 
CRG2B30 
0PG284O 
C  RG  28  50 
CRG2860 
0PG2B70 
CRC2B80 
0RG2890 
CRG2900 
0RG2910 
CRG2920 
CFG2930 
0RG2940 
CRG2950 
0RG296O 
CRG2970 
0RG2980 
CRG2990 
0RG3000 
CRG3010 
0RG3020 
CRG3030 
0RG304O 
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7800  WRITF  (  6,2200)  ISTOP 
GO  TC  800C 

^7900  CONTINUE  £NC  LQCp  THROUGH  NC  CYCLES  ANO  FINAL  CALC  CF  Y 

C  TERMINATE  JUB 

8000  IF  (NC)  86G0, 8600, 8100 

C  CALCULATE  ANO  PUT  OUT  CORRELATION  MATRIX 

8100  WRITF  (6, ICO)  ( T  ITLE <  I  1 . 1 *1 1 201 
WRIT*  (o,2600) 

CT  82 JJ  1*1, NV 

0IAG(I)  =  1  .COO/OSORT(OIAG(I  I  ) 

8200  CONTINUE 
IJ*1 

00  350J  1*1, NV 
CO  9  200  J  *  1 , NV 
RC»( J)= J.CCO 
8 J JO  CONTINUE 

CO  8  AO  0  J*  I ,NV 

ROH( J)*AM (I JI'DIkGI I  I *01  AG (  J) 

I J  =  1 J*l 
8 AGO  CONTINUE 

WRITF  (6,2700)  l , (RCW ( J I . J’ l ,NV ) 

8500  CONTINUE 
0600  CONTINUE 
RETURN 
ENC 
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ORG  3050 

CRG3060 

CRG3070 

CRG3080 

CRG3090 

ORG  2 100 

CPG3110 

CRG2120 

ORG  2  l  30 

0RG31AC 

GRG3150 

CRG3160 

CRG2170 

CRG310O 

3RC-3193 

CRT, 3200 

CRG3210 

CRG3220 

0  R  G  3  ?  3  0 

0  R  G  3  2  A  0 

CRG2250 

CAG  3260 

CRG3270 

CRG2200 

0RG3290 

ORE  2300" 


■»  '-■’-a*" 


uuuuuoouuuuuu  juuuouuujjuudu 


SLBRCUT I  N£  TEST20 


TEST2C 

TEST2 3  TERMINATES  THE  REGRESSION  IF  THE  CHANGE  IN  ALL 
PARAMETERS  IS  LESS  THAN  l.JC-8. 

TEST23  WILL  ALSC  CAUSE  IRcOO  TO  TAKE  THE  INDEX  CF  THE 

PhRAMFTER  WHICH  HAS  THE  GREATEST  FRRCR,  PROVIDED  ICY  HAS 
AN  NCN-ZERO  VALUE  AND  PROVICSO  THE  EPRCR  CN  CNE  OR  MCPE 
PARAMETEPS  IS  GREATER  THAN  ID  PERCENT. 


IMPLICIT  PFAL*fl  (A-H.C-ZI 
R3AL»R  TITLE 

C I  ME  NS  I CN  X<2, 12DOI . YC< 12 JDI ,  SIGYCX  12001 ,P( 1 RC  I  . K 1  <  1R 3 1  , 
1  CCUAJ».P0(14CI,TITLCI201 
QI“'.NSIUN  REF(10,70I,CRRP(  1R0I,  IHJLC<R0I 


COMMCN  BLCCKS 

NAME  ROUTINE  S 


EARTH 

MAIN 

P  3CGR AV , 

0RCLS2, 

FAFW20, 

FUTOUT 

AIR 

MAIN 

G  n  C  L  S  2  , 

TEST20, 

H  A  c  w  20 

PUTCUT  A  NO  LMCLE. 

PIPE 

M4IN 

P*CG'< 

CRCLS2, 

1 1 S 1 2  0 , 

EAFW20 

HATER 

MAIN 

GS iM, 

EAF«2 J, 

YCER  IV, 

OOFR  IV, 

POTJJT ,  ft.  GIF,  AND 

PLCCK  U A  T A . 

TIME 

MAIN 

P3CGP 

PL  TOUT 

SPACE 

MAIN 

EAFR20, 

YCfcRIV, 

COER  IV, 

AND  PUT»UT 

HAN 

MAIN 

PRCG»AH, 

PLTCLT 

CC-HCN  / A I P /  X,  P,  CC,  TITLE,  YC,  SIGYO,  FC,  Kl,  NC ,  NV,  NX,  1W, 
l  NP,  NO,  ISING,  I  STOP ,  IL.  JOFLAG 
CUWMCN  /FIRE/  REF,  ERRP,  IWRITE,  ICY,  IICY,  IPEUC,  IREG,  NHCLO, 

1  IHJLO 
IREUC»0 
ISTAY*3 
J»1 

CC  R03  l-l.NP 

IF  (KKII  .EC. 01  GC  TC  ROQ 

IF  (DABSIPOIJI/Plin.GT.  1.00-81  ISTAY*! 

J*JH 

IF  ( ICY. EC. 01  GO  TC  ROO 
06RRP*ERPP(I  I 

IF  (CERRP.LE.1.0C1I  GC  TO  ROO 
IF  ( IREUC.tC.OI  GC  TC  100 
IF  ( CERRP.LE.ERRP (1REC0II  GO  TO  ROO 
100  CONTINUE 

IF  ( NH0LC.tQ.0l  GC  TO  300 
CO  200  K=  1 , NHOLO 
IF  I  l.EO.  IHOLDIKI I  GC  TO  ROO 
200  CONTINUE 
303  CONTINUE 
IREOC*t 
R03  CONTINUE 

IF  (  ISTAY. EQ.l)  GC  TO  500 
I STJ  P*1 
500  RETURN 


IES 

10 

US 

20 

TPS 

30 

U  S 

RO 

ns 

50 

TES 

60 

T'S 

70 

US 

EO 

Y  r  S 

90 

TCS 

100 

TFS 

no 

Tr  S 

120 

TFS 

130 

•TFS 

1<  0 

Tc  S 

150 

TFS 

160 

Tc  S 

170 

HS 

180 

Tc  S 

150 

■TFS 

200 

US 

210 

TFS 

220 

ns 

230 

US 

2Ru 

t  :  s 

250 

TFS 

26C 

Y'S 

270 

TFS 

2B0 

TFS 

290 

TES 

300 

TFS 

310 

TES 

320. 

-TFS 

330 

TES 

3R0 

TES 

350 

TES 

360 

TES 

370 

TES 

380 

TFS 

390 

Tr  S 

ROO 

TES 

RIO 

TES 

R20 

TFS 

R30 

TES 

RRO 

TES 

R50 

TFS 

R60 

TES 

R70 

"ES 

R80 

TFS 

R90 

TFS 

500 

TES 

510 

TES 

520 

TES 

530 

TFS 

5R0 

TFS 

550 

TES 

560 

TFS 

5  70 

TES 

5  BO 

TFS 

590 

TES 

600 
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SU610UT  I  f.  f  MINV20  I  A  M  ,  N  ,  NFA  IL  I 
IMPLICIT  afAL*a(A-H,C-Z) 

OIMcNSMN  AMI  1001C) 

**********  SEGMENT  1  CF  CHOLESKI  INVERSION  ********* 
*****  FACTOR  MATRIX  INTO  LORCR  TRIANGLE  X  TRANSPCSP 
K*1 

IF  (N-II  ICO  i  2  30  t  3C0 
J  NFAIL*K 
GO  TO  20 OC 
)  AM  1 1 1  *  I  . O/A  M 1 1 | 

GO  TO  19CC 

*****  LOCF  M  OF  AIL, Ml  ***** 

)  OC  1 20 J  M»1,N 
IMAX=M-1 

*****  LOCF  L  OF  AIL, Ml  ***** 

00  1I0J  L*M ,  N 

SU"A=0.0 

KLI--L 

Kw  1  =  V 

IF  I  IMAXI  6Q0,60C,A00 

«****SdM  CVER  1  =  1, M-l  A(L,II*A(M,I)  »•*»» 

I  CO  50  3  1  =  1,  I  MAX 
SUMa  =  SUmA«A.m(kLI  I  *AMIKMI  I 
J-N-I 
KLI=KLI*J 
I  KMI=K«I*J 

*****TERMsC(  L,  M  |-SUM  ***** 

I  T=RM«AM(K  l-SUMA 
IF  (L-M)  700,700, 1000 
l  IF  I  TERM  |  S0C,93C,90C 
*****  A(M,m|=SORT ITERM)  ***** 

I  OENOM=OSC«T1TERHI 
AM (K I »DfcNCM 
GO  TO  HOC 
I  NFAI L*K 
GO  TO  200C 

•**••  A(L.M|=TERM/AIM,MI  «•»•• 

AM(kI»TERF/CENCM 

K*K*  1 

CONTINUE 

**********  SEGMENT  2  OF  CHOLESKI  INVERSION  •»•**•*** 
••♦♦•INVERSION  OF  TRIANGULAR  MATRIX***** 

«M< 11*1. O/AMI 1  | 

K0M»1 

*****  STEP  L  OF  BIL.MI  •♦••• 

CC  1500  L*2,N 
K0M=K0M*N-L*2 

*****  RFC  IPRCCAL  CF  C I  AGONAL  TERM  ••••* 
TERR*1.0/AM(KOH) 

AM|KCMI=TERM 
KM  l»0 
KLI*L 
IMAX=L-1 

*****  STEF  M  OF  81L,M|  ***** 

00  IAJO  M«1  ,  IMAX 
K=KLl 

*****  SUM  TERMS  ***** 

SU“A=3.0 

00  1300  I  =  M,  IMAX 
l I *KMI*I 


MIN  AO 
MIN  50 


MIN  100 

min  no 

MIN  120 
MIN  130 
MIN  I«0 
MIN  150 
MIN  160 
MIN  170 
min  iao 

“IN  IRQ 
MIN  200 
MIN  210 
MIN  220 
MIN  230 
MIN  2A0 
•IN  250 
MIN  260 
“IN  27C 
MIN  280 
MIN  2-JO 
MIN  330 
MIN  310 
MIN  320 
MIN  330 
MIN  3A0 
MIN  350 
MIN  360 
MIN  370 
“IN  380 
MIN  3 <30 
MIN  AOO 
MIN  A10 
MIN  A20 
MIN  A30 
MIN  AAO 
MIN  A50 
“IN  A60 
MIN  A70 
»'IN  A 80 
MIN  A <30 
MIN  500 
MIN  510 
MIN  520 
MIN  530 
MIN  5A0 
MIN  550 
MIN  560 
MIN  570 
MIN  580 
MIN  590 
MIN  600 
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SU“A»SUMA-/SM(KLI  1* AMI  1 1 1 
13 JO  KL I  =  KL I *N- I 

-  **♦*»  MULT  SUM  *  RECIP  OF  DIAGONAL  ***** 

AMR  l=SUMA*TERM 

J«N-M 

KLl=K*J 

1400  KM  I  *  KM  |  *  J 
1  5 JO  CONTINUE 

'  **********  SEGMENT  3  CF  CHOLESKI  INVERSION  ********** 

***.*PKcHUTIPLY  LOhEP  TRIANGLE  8Y  TR  AMSFCSE  ♦**♦  * 

K  3  1 

CO  1800  M*1,N 
KL  I  *K 

CO  1700  L  *M , N 
KMI  =  K 

IMAX=N-L* 1 
SlMA=0.0 

CO  1600  1=1, IMAX 
SUM;=SUMA*Am<klII*AM(kMI| 

KLI=KLI*1 
1600  KMI=KMI*l 
AM(k  l«SUMA 
1700  K  =  K ♦ 1 
18 JO  CONTINUE 
1900  NFAIL'O 
2000  RETURN 
FNO 


MIN  610 
“IN  620 
MIN  630 
MIN  640 
MIN  650 
MIN  660 
MIN  670 
MIN  680 
MIN  690 
“IN  700 
MIN  710 
MIN  720 
“IN  730 
“IN  740 
“IN  750 
“IN  760 
MR  770 
“IN  780 
MR  790 
“IN  800 
"IN  810 
“IN  820 
MIN  830 
“IN  840 
“IN  850 
MR  860 
"IN  870- 
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StBrtOuriivi  SAFW20  (  YC  ,  I  I 


C 

C 

C 

C 

c 

L 

c 

c 

c 

c 

c 

c 

c 

c 


EAFW20 

r.-IIS  Vc'RSICN  WRITER  3/09/72  BY  HAAS. 
LAST  MCCI FISC  8/05/72  BY  HAAS. 


implicit  rfal»8(a-h,c-zi 
lcgical*i  ecc, label 
RLAL«a  title 
CIMlNSIDN  SGNI  2 ) 

01  FUSION  CDEF(6,70),PNAWEI20).TINVI6  ,70,9)  ,  I  PH  A  SF  (  6,  70  I , 

1  NPHAs2<701,lKDi:NT(7C»,IGQ{7J),ISTATE(fc,7G),NINVLM6,70), 

2  I  ISTAT(t»7GI,lNVPF(6,7j,5I»lNVSC(fc»70J 
01  P|  AS  I  2  I  X(  2, 12JGI  ,  Yf  (  I20J  ),  S1SYCU200)  ,P(  190)  ,KI  (  190)  , 

1  C.t  1'.  )J  ,P;i(  19J),T!TLr-<20) 

01*'-.  AS  I  JiS  '!^F(1J,7J),FFA0(  190) , IH0LCI90I 
0  I  Mc  A  S  I  J  I  3CINV(2).$TCCr.F(lj)#YiiSNH2),TK(I9,2)*TYPuI19)»NSCALE(5I 
1  ,LA3„t.(E0) 

CIHENSI ON  COCCI  7) ,CYCC( 71 


EAF 

FAF 

EAR 

EAF 

EAF 

OAF 

R(F 

CAR 

£AF 


10 

20 

30 

90 

50 

60 

70 

80 

90 


C C N u C N  PLCCKS 

NAPE 

ROUTINES 

EARTH 

MAIN  PROGRAM , 

CAGLS2, 

RAF, 20,  PUTCLT 

AIR 

MAIN  PROGRAM, 

CRCLS2, 

Tt  ST2  3 ,  FAFW23 

FIRE 

PUTUUT  ANC  UN  I  CLi  E  • 

MAIN  PROGRA*, 

0RGLS2, 

TFST20,  L  AF  H  20  1 

fcATf-R 

MA IN  PROGRAM, 

EAFW20, 

YCEHV,  DC  ER  I  V  , 

TIME 

PUTUUT,  FN,  C I 

F  ,  AND 

PLCCK  DATA.  i 

MAIN  PROGRAM, 

PLTOUT 

1 

SPACE 

main  program. 

EAFW23, 

YCERIV,  CDERIV,  1 

MAN 

AND  PUTPUT 

| 

main  prcgram, 

PL  TOUT 

1 

CCXNCN  /EARTH/  CCEF,PNAME,TINV,IPHASE,NPHASE,IKCLNT,IGC,NSETS. 

1  I  STATE  i  MNVF.R  ,  INST  AT ,  INVPH,  INVSC*LISTP 
COKMC'l  /AIR/  X,  P,  CC,  TITLE,  YC,  SIGYO,  PD,  KI,  NC ,  NV,  NX,  1H, 
1  NP.  NJ,  ISIHG,  ISTOP,  IL,  JOFLAG 

1CIH0LC  /F,Rf/  'tEF’  ERRP’  •lHR,TE*  IC¥'  HCY,  IRF.CC,  IREG,  NHCLD, 

CCPMCN  / m /T ER/ ZERO ,  CNE,  THC,  THREE,  FOUR,  SIX,  R,  F, 

1  SCINV,  TRFF,  STCDSF,  01E0,  ADI E ,  SCI E ,  THETA,  YESNO,  TK, 

2  ASTAR,  TYPE,  Nl  ,  N SCALE ,  LABEL,  BCC 
COMMON  / S  FACE/  DCCC,  OYCC,  SC,  TO 
YC*ZERO 

CO  103  K=1,NP 
CC IK  I - Z ERC 
100  CONTINUE 


DETERMINF  J  BY  LOCATING  I  IN  THE  VECTOR  IKCUNT.  IF  NOT 
LOCATED,  PRINT  ERROR  MESSAGE  ANC  STOP 

CO  200  JM.NSETS 
IF  <  I.LO.  IKCUNTIJII  GC  TO  300 
200  CONTINUE 

WRITE  (6,20001 
STOP 

300  CONTINUE 


EAR  100 
c AF  lio 
Sip  120 
EAF  130 
-'AF  190 
RAF  150 
EAF  160 
EAR  173 
RAR  130 
EAR  190 
R-AF  2JJ 
CAF  210 
EAF  220 
FAF  230 
RAF  290 
CAF  2  50 
SAF  260 
PAR  270 
EAR  280 
EAF  290 
EAF  300 
EAF  310 
EAF  320 
EAF  330 
EAF  390 
EAF  350 
FAR  360 
EAF  370 
EAF  380 
EAF  390 
EAR  900 
EAF  910 
CAF  920 
EAF  930 
EAF  990 
FAF  950 
EAF  960 
EAF  970 
EAF  980 
EAF  990 
EAF  500 
EAF  510 
EAF  520 
EAF  530 
■RAF  590 
EAF  550 
F.AF  560 
EAF  570 
EAF  580 
EAF  590 
EAF  600 


O  O  O  U  O  J  O  O  U  O  o  o  O  JU-JU  O  U  O  O  u  O  O  U  u  O  U  O  u  u  u 


WF  I T  c  REFFRFNCE  IF  THF  1-TH  OBSERVATION  IS  TFE  FIRST  IN  THF 
J-TH  DATA  SET. 

IF  I  IWRITE.FO.l.CR.I.FO.ll  GC  TC  400 
FC  TC  500 
400  CONTINUE 

WRITE  (6.2100)  IR£F(  H  ,  J  I ,  I  I  =■  l  ,  1 0  I 
INKITSU 

500  CONTINUE 

IF  (  I. SO.  IKGUNTI  J  )  I  I  WHITE”). 

IF  (  ISO ( J  I.LE.7)  CC  TC  600 

IF  ICC(J)  IS  GREATER  THAN  7,  CALL  UMCUE  10  CALCLLATC 

YC  (ANC  THE  OEF IVATI ViS  UC  (  III  IF  JCFLAG  (JN  IN  PHAS20)  IS 
ZTR  J. 

CALL  UNIOLF  (YC.I  ,J. JCFLAG) 

?sru°N 

60C  CONTINUE 

FOR  EACH  PHASE  IN  THE  J-TH  OATA  SET,  CALCLLATF  THE  SEVEN 
OcR I  VAT  IVES • 


r= x i  i,i) 

TC=X(2,  I) 

L A  ST  =  NPHASE ( J ) 

CC  18JJ  LU.LAST 


LCCATE  THE  FIRST  CONSTANT  FUR  THE  L-TW  PHAS=  IN  THE  VECTOR  P. 
INITIALIZE  SC,  THf  STC ICH1CMFTRIC  COEFFICIENT. 


IN0uX=1h7»(IPHASE (L, Jl-l) 
SO  =C  CEF  I  L  ,  J  ) 


CALL  YOfcK I V  TO  CALCULATE  THE  TEMPERATURE  TERM  IN  THE  DERIVATIVE. 
IF  THE  DATA  SET  IS  FOR  ENTROPY  ANC  HEAT  CAPACITY,  GC  TO 
STATEMENT  13C. 


IGCES’IGO (J ) 

IF  ( I  ST AT  E (L.JI.EC.-1I  IGOES  =  IGOCS»7 

CALL  YOnRIV  IT.IGCESI 

IF  (  I  GO  I  JI.LE.2)  CO  TC  1600 


CALL  OUE  K IV  TO  CALCULATE  THE  CLRIVATIVe  WITH  RESPECT  TO  -C-  AT 

2)8.15  KELVIN  IF  COMPONENT  IS  AN  ELEMENT.  PROCEEO  TO  ACJUST 
OC  ACCCROINGLY. 

IF  (  ISTATEIL.J  I.NF.l)  GO  TO  COO 
CALL  OOcR  IV  (TREF  I 
CC  ZOO  K«l,7 

CCU ♦INDEX-1  l*OC<K*INCEX-l I+SC*0C0CIK)*0VCC14> 

700  CONTINUE 
GO  TC  16QC 
BOO  CONTINUE 

CHECK  FOR  INVERSIONS.  IF  LCH  TEMPERATURE  PHASES  ARE  TO  eE 
CONS  ICCREO ,  INITIALIZE  SGN.  IF  LCW  TE'/0LRATURE  PHASES 
ARE  NOT  TO  BE  CONSIOEREC,  GO  TO  STATEMENT  HO. 


IF  (NINVERI L ,J) .LE.OI  GC  TO  1600 


•EAF  610 
E  A  F  620 
EAF  630 
EAF  640 
S Ac  650 
FAF  660 
EAF  670 
EAF  680 
EAF  690 
FAF  700 
EAF  710 
EAF  720 
■EAF  730 
EAF  740 
EAF  750 
EAF  760 
■EAF  770 
FAF  700 
FAF  740 
faf  aoo 
■caf  an 

EAF  B20 
EAF  030 
•EAF  840 
f*.F  050 
FAF  860 
FAF  070 
FAF  BRO 
•EAF  890 
EAF  900 
EAF  910 
•EAF  920 
EAF  930 
EAF  940 
•EAF  950 
EAF  960 
EAF  970 
EAF  980 
■EAF  990 
E  A  F 1000 
EA  F  10 10 
EAF  1020 
EAF  1030 
•  E  AF 1040 
EAF  1050 
FAF 1060 
EAF 1070 
-EAF1O0O 
EAF 1090 
EAF 1100 
F  AF 1 1 10 
EAF 1120 
EAF1130 
EAF1140 
EAF  1150 
-EAFH60 
EAF 1 l 70 
F  A  F 1 180 
E  AF 1 190 
■EAF 1200 
E  A  F 12 10 
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ooo  O  O  o  o  O  o  o  o  o  r-  o  oooor^  ooooo 


n  SCO  LL l *1 ,2 
SGN(  LLL  «  MCIKVtLLU 
900  CONTINUE 

IF  (INVSCU.JI.EC.OI  CC  TO  1000 
SGN(  U»SGM  11*STCCEF(  INVSC  (  L,  J)  I 
SGN(2I=3CM2I*STCCEF<  INVSCtL,  J  I  I 
1JO0  CONTINUE 


IF  THE  LOWEST  TEMPER  ATURF  PHASE  IS  AN  ELEMENT,  CALL  DOEP.IV 
TH;  ARGUMENT  IS  298.  15  ,  ANO  PACCEEC  TO  ADJUST  TC 
ACCORDINGLY. 


IF  (  INSTATCL.JI.FC.OI  GC  TO  1200 
1 1 NVP«1 *7»  C INVPhILt J, 11-11 
CALL  OJEPIV  (TRFF1 
00  1  1  J  J  K*l,7 

DC.  (<  »l  INVP-ll>OC(  MI  IKvK-1  )  »OQdC  (K  I  ♦SCN(  2  I* SC»DYDC <  4 1 
1  100  COM  INud 
1200  CONTINUE 


r.AF  1220 
c  A  F 12  30 
F  A  F 1240 
F  AF  12  50 
F  A  F 1260 

E  AF 12  70 
EAF 1280 
“  A  F  1 2  9  0 
EAF 1 300 
F  A  c  1 3  1 0 
0  A  F  1 3  2  0 
F  A  F 1330 
0  Ac 1 340 
" A  f 1 3  5Q 

FAF 1360 
F  A  F 1 3  70 
OAF  1300 
FAF 139J 
rAF 1400 


- -  j  e  m  j 

CALL  00  E3  IV  TO  CALCLLAT  F  THF  OHRIVATIVc  u  I  T  A-  A  f  SDF  C  T  TO  D  FOP  FAF1420 

EIOE  ANO  LOW  TEvflFATUK:  PHASE  AT  T.-EPC  =  ATLPF  of  f  c  f  14  30 

INVERSION.  FRCCEtO  TJ  AOJUST  CC  ACCCUINGLY.  EAF1440 


c 


LSTI  NV  =  NI  AVEML,  J  I 

CC  1 50 j  LL=t,L'TIM 

CALL  OJERIV  ( T  I  'I  V  ( L  ,  J  ,  LL  I  i 


IF  (INVSCIL.JI.ee. 01  GC  TO  1300 
IF  (LL.NE.LSTINvl  CC  TC  1300 
SGN(2I=SGM2»/STCCSF<  INVSCIL,  JI » 

1300  CONTINUE 

CC  1 40 J  L  LL  - 1 , 2 

IINVF=1*7MINVPH(L.J,(LL»LLL-11  )-l  ) 

CC  L  40  J  K*l,7 

0L<MIINVP-1I=0C(MIINVR-1  I  ♦0QaC(K»*SCNIlLL)*SC*DYCC<4» 
1400  CONTINUE 
1  500  CONTINUE 
1600  CONTINUE 


COMPLETE  THE  CALCULATION  OF  OC  BY  ACCINC  TEt  DERIVATIVE  AT 


CO  1700  K*l,7 

OC(K  M.N0EX-1I»DC(  KMNOEX-l 1  ♦SC*CYCC(K> 
1700  CONTINUE 
1800  CONTINUE 


ALL  CERIVATIVcS  HAVE  BEEN  CALCULATED.  NCW  CALCULATE  YC  BY 
SUMMING  THE  FRCOUCT  CC»P. 


CO  1900  K=1,NP 
YC  =  YC*P(K  I  ♦DC ( K I 
CONTINUE 


RETURN 


RETURN 


- E  A  F  14  50 

FAC  1460 
E  A  F 14  70 
E  A  f  14  8  0 

- F  £  F  14  90 

FAF 1500 
EAF 1510 

- FAF 1520 

FAF 1530 
EAF 1540 
EAF 1550 
EAF1560 
E  A  F 15  70 
E A  F 15  80 
E  AF 1 590 
F  A  F  1 6  0  0 
EAF161Q 
EAF  16  20 
FAF  1630 

- EAF1640 

TIOBS  I  .E'AF  1650 

- FAF 1660 

FAF167G 
EAF16B0 
E  AF  1690 
EAF 1700 
- EAF  1710 

F  A  F 1 720 
FAF 1730 

- EAF 1740 

EAF1750 
CAP 1760 
E  A  F 1770 

- FAF  1780 

EA  F  1790 

- E  AF 1800 

c  A  F 18 10 
- F  AF 1 8 20 


f 4F1Q30 
— -EARIBRO 
RVATEAF 1850 
IS  FAF I860 
E  AF 1870 
FAF  1880 
EAF 1890- 


P 


SUBROUTINE  YCERIV  <T,J) 


YOSRIV 

THIS  VERSION  WRITFN  3/09/72  BV  HAAS. 
LAST  Rf  VI  SEC  8/C5/72  BY  HAAS. 


IMPLICIT  REAL*8(A-H,C-7I 

Cl  MENS  I JN  CYCCI7) , SC INV I  2 ) , CCDC I  7 ) , STCOeF I  1C) 


CALCULATE  THE  APPROPRIATE  TEMPERATURE  TfPH  IN  THE  DeRIVATIVF 
AMO  RETURN.  —  APPROPP.I  ATE  AS  INDICATEC  PY  J. 


YCS 

-YCE 

YCF 

YCE 

YCE 

YCE 

YCP 

•YCE 

YCE 

YD! 


10 

20 

30 

AO 

50 

60 

70 

80 

50 

100 


c 

■ — YDE 

110 

c 

CCMMCN  BLOCKS 

YCE 

120 

c 

NAPE 

ROUT  INFS 

YCE 

130 

c 

EARTH 

MAIN  PPCGRAm,  CPCLS2, 

EAFR2Q, 

PLTOLT 

YCE 

140 

c 

AIR 

D4CGpAMt  C3GIS2» 

TCST2J, 

E  flc  W2C 

Y  Cr 

150 

c 

PUrCUT  ANC  UNICLE. 

YDr 

160 

c 

FIRE 

MAIN  PROGRAM,  C  c  C  L  S  2, 

TFST2J, 

EAFH20 

YCE 

170 

t 

RAT  LR 

main  program,  eafa20. 

YD*MV, 

DCC» IV, 

YC' 

180 

c 

PUTCUT,  FN,  C I c ,  AND 

PLCCK  DATA. 

YCE 

190 

c 

TIME 

MAIN  PR  JG°  A  w ,  PLYCUT 

YCE 

2  DO 

c 

SPACE 

MAIN  PROGRAM,  CAF«2D, 

vri-RIVt 

DDFB IV, 

YCE 

210 

c 

ANO  PLTPL'T 

YCE 

220 

c 

c 

PAN 

MAIN  PRCGPA«,  PLTCUT 

YDS 

230 

■—YCE 

240 

CO“MCN  /RATFR/ZfRC, 

CNF ,  thc,  three,  fclr,  six 

.  R,  F, 

YCE 

250 

1  SC  INV,  TRCF,  STCCFF.CILC, A,B, ThETA 

YCE 

260 

c 

CCMMCN  ,'SFACC/  CCCC 

.  CYCC,  SC,  TC 

YCE 

—  YCE 

270 

280 

YDE  29D 
YCE  300 

- -  3l(J 

GC  TC  HOC,  200. 3CC,40D.  500. 600,  700.  SCO.  1CC0,11»0,  1200, 1 3CC  *  1A00  *  1 5  YCE  320 


ICO),  J 
100  CONTINUE 


THE  CFRIVATIVE  FOR  HEAT  CAPACITY  CF  NCN-ICMC  SPFCIES. 

DYCC (I )*CNF 
CYCC I2)*TV0«T 
CYCC(3)=0NE/(T«T) 

DYCC I4)*ZER0 
CY  CC 1 5 )  *7  FRO 
DYCC I  61 *T*T 
DYCCITI'ONE/DSQRTIT) 

BETUPN 

200  CONTINUE 

THE  CFRIVATIVE  FCR  ENTRCPY  GF  NCN-ICMC  SPECIES. 

DY  CC(l)«ClCG!T) 

DYCC ( 21  *  T VC*  T 

CYCC I31=-CNF/ITW0«T«T) 

CYCC  14)  *Z  ERG 

DYCC (5)»0NE 

DYCC(6)=T«T/TH0 

DYCC  (71=-TV«0/DSCRT(T1 

RETURN 

300  CONTINUE 

THE  CFRIVATIVE  FCP  THE  ENTHALPY  Cf  NCN-IONIC  SPECIES. 


YCE  330 
YCE  3A0 
-YCE  350 
YCE  360 
-YDE  370 
YCE  380 
YDE  390 
YDE  AOO 
YCE  A10 
YCF  A20 
YCE  A30 
YCE  AAO 
YDE  A50 
YCE  A60 
-YCF  A70 
YCE  A80 
-YCE  A90 
YCF  500 
YDE  510 
YCE  520 
YDF  530 
YCE  5A0 
YDF  550 
YDE  560 
YCE  570 
YDE  580 
-YCE  590 
YCE  600 
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o  o  o 
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CYCCIll’TC 
rycc(2i=TC*TC 
CYCC(3I=-CNE/TC 
CYCC ( 4  I =CNE 
CYCCI5I* ZERO 
DY0C(6I=(TC**3)/THRE£ 
CYCC(7l=Tk0*0SQRT(TQI 
ocrjPN 

<503  CONTINUE 


THE  CEP.IVATIVf  FOR  THE  HEAT  CAPACITY  CF  ICMC  SPECIES. 

CYCC ( 1 1 *CNF 
C  YCC  121=  H3*T 
CYC-I3»=>CA6/IT*TI 
DYCC(4l*ZERO 
CYCC(5I=ZFRC 
CYCC 1 61  *  T*T 
FNT*  FNI T ) 

CYOCI  M=-T*CIEIT*»I  4*ent*a«fnt*a*fnt*a*a*fnt*twc/th;:ta*jnc/i 

ITH'TAI I 
RETURN 

1JG0  CCNTINUE 

THE  CFRIVATIVE  FCP  THE  LNTECPY  CF  ICMC  SPECIES. 

CYCC ( 1 1 -C ICG ( T ) 

PYRCIZIMkOM 

CYCU31*-CNC/|TWC*T*T» 

CYCC ( 41  =  Z  ERO 
CYCC(5I=CNE 
CYCC 1 6 1  3  T*T /TWO 

CYCC  1 7 1  =-C  I G  t  T I«(A*FN1 T I  PONE/THETA) 

RETURN 
CCNTINUE 


1100 


THE  CERIVATIVF  FCP  THF  ENTHALPY  CF  ICMC  SPFCIES. 


CYCC 1 1 1 “ T 

DYCC(2I=T«T 

CYCC 1 3 I =-CNE/T 

CYCC I  4  I =CNE 

DYCC(5I»ZERC 

CYCC(6l=(T**3l/THRcE 

CYCC I  7 » =  C IF  I  TI*1CNG-A«T*FN<TI-T/TFETA» 
RETURN 

12J0  CDNT1NUL 


THE  CERIVATIVF  F CR  THE  GI3BS  ENERGY  CF  ICMC  SPECIES. 

CYCCI1I* T-T*DLCG(T» 

CYCCI2I— 7»T 

DYCC(31=-CNE/(TW0*TI 

CYCC(4I*CNE 

CYCCI51 «-T 

CYCC16I=-(T*»3I/SIX 

CYCC  (  7I  =  C  IEIT  I 

RETURN 

1300  CONTINUE 

;  THE  DERIVATIVE  FCF  THE  ECU11IBR1U"  CONSTANT  CF  IONIC 


YCE  1220 
YOC  1230 
YOG1240 
YCE1250 
YCE1260 
YDE 1270 
YCE128Q 
YCE  1290 
YCE1300 

- YDE1310 

YCF1320 

- YDF1330 

YC.31340 
YCM350 
YOG  1 3fc  0 
Y06 1370 
YCE 1380 
YTE  1390 
YCr 1400 
THCTA*Y0“1410 
YCE  1420 
YCE 1430 
YCE  1440 

- Y  D  E 1 4  5  0 

YCL  1460 

- Y  C  “  1 4  7  3 

YCE  1480 
YOP 1490 
YCC  1500 
YCE 15 1 0 
YOF  1520 
YCE  1530 
YOG  1 54  0 
YCE  1550 
YCE  1560 

- - YCE1570 

YCE15eo 

- Y  C  E 1590 

YOE1600 
YCE1610 
YDE1620 
YCE  1630 
YC  G 1640 
YCE  1650 
YCE  1660 
YCE  1670 
YCF1680 

- YCE 1690 

YC  f  1 70  0 

- YCE  1710 

YDE1720 
YCE  1730 
YCC 1740 
YCE 1 750 
YC  L  1 760 
YCE  1770 
YC  E  1780 
YCE  1790 
YCE  1800 

- YCE  1810 

YCE  1820 


k 
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1 1  .. 


o  c  o  o 


C  SPECIES.  YrE1830 

c  _ YCF18A0 

DYCC11I=1CL0G(T)-CNP)/R  YC'1850 

CYCCI2IM/P  YCE1860 

CYCC1  3l=CNE/(  TWC*P*T*T)  YHE  It  70 

0y:C(4I*-CNF/(8*I)  YH*:  1880 

CYCC  (5)  - CNE / R  YC:i890 

CYc;i8»an*n/i  s  i  x*r  »  yc«1900 

DYCCI  7M-CIE  IT  )/«R*T)  YCC 1910 

5  ?  T JPN  YCS1920 

19J0  CCNTINUE  YC  0 19  30 

c  _ YCE19A0 

C  THE  CEPIVATIVf  FCR  THE  CELL  PCrENTIAl  CF  ICMC  SPECIES.  YC“1950 

c  _ YCP1960 

CYCC11MIT*CL0G(  Tl-TI/F  YC197C 

CYCCI2MM*T)/F  YC*  1980 

CY:cm  =  CSE/(Tw3«F*T)  VCC1990 

CYCC (9  I =-CN£/F  Y^C 2000 

CYCl!5I=T/F  YC'2J10 

CYCCI6I= IT*«3)/I S IX*F  i  VC?  2020 

DY3C 1 7 ) *-C I i ( T ) / F  YD- 20 30 

RETURN  YCc  2090 

1500  CCNTINJE  YCE2050 

_ _ 

THE  CFRIVATIVE  FOR  THE  RELATIVE  FEAT  CONTENT  CF  ICMC  Y  C c  2  0  70 

SPECIES.  YCE2080 

- - - YCF2990 

IF  (SC.LT.O)  GO  TC  1 6C0  YCE2100 

CYCC 1 1 1  *T  Y0E2110 

CYCC  I  2 1  =*  T  ♦  T  YCE2120 

CY  CC  (  3  t  =-CNE/T  YCE2130 

OYCCIAMCNE  YC:2lA0 

OYCC 1 5 ) *  7  SRO  YCE2150 

DYCC(G)=<T**3)/TFPCE  YCE2160 

UYCC(7l=niE( T)*ICNE-A*T*FNI TI-T/TFETA)  YCE2170 

RETURN  YCE2180 

1600  CONTINUE  YCE2190 

CYCC ( 1 1 *TC  YCE2200 

CYCC  ( 2 )  =  T  C*TO  YCE2210 

OYCC  (31—CNE/TC  Y0E2220 

OYCC(A)=UNE  Y0E2230 

CYCC ( 5 )  =  7  f RO  YCE2290 

CYCC(6)=(TC**3)/TFREE  YCE2250 

CYCC I 7)»C IE  (  TO) *ICNE-A*TO»FNI TOJ-TQ/THETAI  YCE2260 

RETURN  YC52270 

c  _ _ 

EMC  YCE2290- 
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I 


c 

c 

c 

c 

c 

L 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

t 


c 

c 

c 

c 

c 

c 


SUBROUTINE  CCER  I  V 

(  T  1 

CCE 

10 

CCE 

30 

CCER  IV 

OPE 

60 

CCS 

SO 

THIS  VERSION  WAS 

VRITEN  3/09/72  BY  HAAS 

• 

CCE 

60 

LAiT  REV  I  Sc  C 

8/ 05/72  BY  HAAS. 

CCE 

70 

CCE 

60 

IMPLICIT  RcAL«8(A 

-H.0-2I 

CCC 

100 

01  He  NSI ON  OCCC 1 7) 

tSC INVI2 1 . OYCCI 71 

CCF 

110 

COMMCK  BLOCKS 

OCE 

130 

N  AMH 

ROUTINES 

CCE 

160 

EARTH 

MAIN  PRCG^AV, 

CRGLS2,  EAFW20, 

PLTOUT 

CCC 

150 

AIR 

PAIN  PPCG’JAP, 

CPCLS2,  TEST2J, 

EAFW20 

CCE 

160 

PUTCUT  ANT  UNICUF. 

ere 

170 

f  IRe 

MAIN  PACG^Av, 

0®CLS2,  TEST2J, 

K AFw2C 

180 

WATER 

MAIN  P ^CG P A ^ , 

cJFw20,  YCrRIV, 

C OCR  I  V i 

cr.E 

190 

PUTDur,  FN,  Cl 

F.  AND  PLCCK  OATA. 

ere 

203 

TIME 

MAIN  PRCGpax, 

PLTCUT 

CCS 

210 

SPACE 

MAIN  PRCGPAM, 

E<FW20,  YCERIV, 

DOER  IV, 

CCE 

220 

AND  PUTPLT 

CCE 

230 

^  AN 

MAIN  P^CGTA", 

PLTCUT 

COE 
—  CCP 

260 

250 

CHXMCV  /WATIV/ERC,  ONE,  TwC,  THRCc,  FOUR,  SIX,  R,  F,  CCC  260 

1  SCINV,  TCCF,  STC17F,  Die?,  4.  P,  THETA  COE  270 

cj"i:CN  /space/  occc,  cycc,  sc.  to  ere  200 

- CCC  290 

CALCULATE  THr  DERIVATIVE  OF  -0-  WITH  RESFECT  TC  THE  CTHCR  OCE  300 

CONSTANTS  WHERE  THE  GIBBS  ENERGY  IS  C.O  AS  IS  TH=  OCE  310 

CASE  FOR  ELEMENTS  AT  29B.15  KELVIN  CR  AT  THE  TFWP'RATuPc  OF  CCE  320 


nvEPSICN.  OCE  330 

- CCE  34Q 

CCCC<1I=T«0LCG<T l-T  COE  350 

C0CC<2I*T«T  CCE  360 

C0CC<3I=CKE/ITW0*T»  CCE  370 

OCCC (61 a-CNF  DCE  3B0 

OOCC  <  5 1  »T  COE  390 

C0CC(6I=(T«*3I/SIX  CCE  600 

CCCC(7I=-FCUR«0SCPT(T»  OCE  610 

RETURN  CCE  620 

END 
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t.W^.  .....  - -  ■  . . .. - -  .. .  ' - — 


*  'vKSW-  t 


C0U3LE  PRECISION  FUNCTION  0IE(TI 

Cl* 

10 

IMPLICIT  P*AL*a(A-H, 

1 1 

Cl* 

2C 

CIPcNSION  STCOFF ( 101 

,  SC1NVI2) 

CIS 

30 

40 

c 

COMMON  BLOCKS 

Cl* 

50 

c 

NAME 

ROUTINES 

Cl* 

60 

t 

Earth 

MAIN  FRCG“AM, 

CRGLS2, 

C  A  F  w  2  0 , 

PU (OUT 

CIS 

70 

c 

AIR 

MAIN  PROGRAM, 

CRCL  S 2, 

TCST20, 

EAFW20 

Cl* 

80 

c 

PUTGUT  ANO  UNI  CUE. 

Cl* 

SO 

c 

FIRE 

MAIN  PROGRAM, 

C  0  CL  S  2, 

TCST20, 

6  AF  W 2  0 

DIE 

100 

c 

RATCR 

MAIN  PROGRAM, 

E  A  F  M  2  0  , 

TCERIV, 

DCFRIV, 

C  1  P 

no 

c 

PUTUUT ,  PI,  Cl 

S,  A  NO 

PLCCK  DATA. 

C  I  * 

120 

c 

TIME 

MAIN  PRCGRAM, 

PIT  JUT 

CIS 

no 

c 

SPACE 

MAIN  P^CGPft*, 

EAFrEO, 

VCER IV, 

CCER  IV, 

Cl* 

l'.O 

c 

AND  OUTPUT 

CIO 

150 

c 

c 

MAN 

MAIN  P&CG^i*, 

PLTCUT 

Cl* 
—  Cl  * 

160 

170 

common  /-.ater/ecrc,  cn-  ,  t«c,  -hafc,  pout,  six, 

1  SCINV,  TR'F,  STCCCF,Ol.iO,A,e,  THE^A 

ci  F=r:xp(  fmt)  ♦t/th?T4  i/oi'iu 

F.  F  T  U  0  N 
CNC 


3i*  iuo 
Cl'  ISO 


o i r  200 

Cl*  210 
CIS  220- 


O  o  ‘ J  OJOOOOOOOOO 


DOUBLE  PRECISION  FUNCT ION  FN I T  J  FN  10 

IMPLICIT  «CAL»R(  A-F,  C-2 (  FN  20 

Cl  PENS  ION  STC')EFtlO»,SCINVt2l  FN  20 

- -  'g 

COMMON  BLOCKS  FN  50 

NAME  ROUTINES  FN  60 

EARTH  MAIN  POCGFAM,  0RCLS2,  bA-.Cj,  PUT"LT  FN  70 

"AIN  PROGRAM,  C=GLS2,  TE5T0J,  EAFW20  FN  80 

PUTGUT  ANC  LMCU.  FN  90 

EIRE  MAIN  PRCORA*,  0FGLS2,  TES'GJ,  EttF»20  FN  100 

t-ATER  MAIN  PROGRAM,  EAFW20,  YCERIV,  CCERIV,  FN  110 

PUTGUT,  FN,  C I E  t  AND  FLCC \  0  A  T  A •  FN  120 

TI"E  MAIN  PCCGPAM,  PLTCUT  Fn  130 

SPACE  MAIN  PROG  PA v  ,  f  A  F  W  20  »  YC  I  V  >  CD£RIV,  FN  no 

ANC  PUTPUT  CN  150 

NAN  MAIN  PRCGRA-,  PLTCL  T  FN  160 

- -  17() 

COMMON  /  n  A  T  V  R  /  2  F  R  C  »  CNc,  T„C,  TH°LC,  FCUP  ,  SIX,  ft,  =,  FN  180 

1  SCI'JV,  T'CF,  STCCEF  ,C1L0,A,  B,  THETA  FN  190 

fn =3"XP I B  •  A*T 1  FN  200 

R-TUPN  fN  210 

ENC  FN  220- 


onooooor»oc'Oc'<on 


IW'.n  p: -r 


SLBRCUTINf  UMOUE  ( YC,  1 ,0  ,  JCFLAG  ) 


C'JVHY  SUBROUTINE  LNILLE 


SCTLRNS  ERRCP  STATEMENT  IP  CALLEC. 


UNIQUE  is  TJ  Be  PPCGRAPMEO  BY  THC  u"P  imcip  k 

T1  5UPPLY  YC  ANC  CC  (K  I  (  NP  ,  .^O  *  FCP  X?  ,'c4|  A  \^V 

C0V..PC0  BY  THE  PCLTIMES  E  AF  w  20  ,  YCETIV,  AAOCol  U. 


™E  f UNICUf  WILL  BE  CALLEC  WHEN  TH  = 
GREATER  THAN  7. 


INCEd  IGQ(J)  IS 


ARGUMENT  S- 
YC  — 


UNI  30 
LM  A  0 


LNI  100 
LM  110 
UM  120 


‘  VC  ”  '“THE0  pPAMP,Mp;sVp<nCU-  FP0M  C°^ENT  ViLUES  CF^>  130 

l  \  "  [(N”X  POP  T  H  *  I-TH  CBSEPVATICN.  !  }?“ 

j  —  piQPX  F(J i*  IN'!  J-TH  SFT  ^ 

■  JCFLAG  -  l*  UNIOUC  IS  ciLLfC  »T%EG»««5ir%  TIP«  n! 

=  ,F  is'lnl  ‘I  nLr  JCFLiG  «N  «S 

:  _ _ ncVbe  0t",V4T,v6s  NEED  UM  210 

IMPLICIT  PEAL»B(A-H,c-2l  UN  I  230 

C I  «i  NS  I  Qf,  Cn:p(6,7:),PNAME(20),Ti;,V(C,7''  41  IPI-PPI^  7'.  tJNI  2<*a 

NPHASP(7j),IK:uNT(70t,IGoi7o;,ISTAU.;  Jj!  N  ,Kvr  ‘  ?  ^  I  UM  250 

i  I  I  ST  AT  ( l ,  70  )  >  I N  VFH  (  fc  ,  70, 0  I  ,  I  vv  SC  (  t ,  7  0  I  '  '  '  UM  260 

- - * - - - - - 1 _  UM  270 

CUHPCN  BLCCKS  UN  I  280 

N«PE  ROUTINES  UNI  290 

mTH  “mw  '  £««*.  E*F-2J.  mnur  lm  310 

t£m  Tesi?0*  r,Fh?0  ®  ™ 

WATER  HA,N  PRCCPAM  l.ftH  SUffc  ft  { 

asrwsfirfttSi  puck  °*ti-  ™  -° 

m%?88a 64FW20' 00£fi,v' 

_ _  NAIN^PhCGPAA,  PLTCLT  ft*  ”0 

if f j;;p^sss?^^nn;  s;{  :is 

WRITE  (6,1001  J ,  J ,  I GC  ( J  )  *(NV5e,LISTP  UNI  430 

STCP  UNI  440 

1U0  POQMAT  ( 6F0YQUR  .I3.39H-TH  DATA  SCT  CALLEC  iMnic  e  ,r»  r.r.™  ,  UNI  *50 
1 1  GO  (  (12,1 7HI  HAS  A  VALLE  0*  .16. IBM  n.iT  ,  *  .  rui  ?!  I  C  EiFW™./5h  LNI  460 
2RAM  TtO,  )  *  1 6 , 38M  OUT  LM  CUE  IS  CURRENTLY  UNPPCGUM  470 

ENC  UM  A80 

UM  490- 


N»P£  ROUTINES 

aE?rPTH  £!!  •,CLS?‘  E*F-2J*  "UTOUT 

1  M!  p'fnRA».  CRCL«2,  T'STTO,  F  AFW20 

PUTnur  A«.U  LMCLE. 

f!JE  *'4IN  *■*«***•  dCLSE,  TEST 20 ,  EAFW20 

h  TER  PRCGAAP,  EAPW20,  VCERIV,  ODEPIV, 

_  PUT  IJT,  PN,  CIM  AND  PLCCK  DATA. 

**»*•*  PRCG5  AP  ,  PLTCLT 

SPACE  RRCGPA-.  6APW20,  rct-fllv,  ODEPIV, 

AND  PUTPUT 

,'AN  NAIN  PRCGRAP,  PLTCLT 

WRITE  (6,1001  3,J,IGC(J) 

STCP 
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t&vmv-KjTPfr** 


IJI X i.  jj, .  .pi*  WJf.K*  WM^.WJ  J. 


,u'|,ujs,i,'.i.'.i^nw 


SUfiRCUT I NC  PUTCUT  IJI 

PUT 

10 

IMPLICIT  REAL*8(A-H,C-ZI 

PUT 

20 

LCC I  CAL* 1  ecc, LABEL 

PUT 

30 

r.EAL*9  TITLE, ACCFF 

PLT 

90 

CIMfcNSIJN  SGNI2I ,A<7,7) 

PUT 

50 

DIMENSION  C3FFI6,70l,PNAMEI20l,TINVIt,7 

:.9I,IFHASEI6,T0I, 

PUT 

60 

1  NPHAjrl 70 1, I  MOUNT ( 7CI 

,  I  GO  17.11  , 1  STATE  (t  ,  7C  1  ,NIN  Vf  » 16  , 

70), 

PUT 

70 

2  I  ISTAT  It, 73),  INVFHI6, 

70,5 1 .INVSCI t , 7  0 

I 

PUT 

80 

DI  Me  NS  I  IN  XI  2 , 12  CCI  ,  TCI  1200  1 ,  SIGYC  I  12C0I  .FI19CI 

.  K  I  <  190), 

FLT 

90 

1  CCI19J1, PDIl 901, TITLE! 201 

PUT 

100 

Cl  MENS  UN  SCINVI2I  .STCQtFI  10 1 ,  V  E  S.NO  I  2  1 , 

TK (19, 21 

,TVPE( 19  I .NSCALEI 51  PUT 

111 

1  .lABEUECI 

PUT 

120 

DIMENSION  CDCCI 7I.CVCCI7I 

PUT 

130 

DIMENSION  AA  (1901 

PUT 

1 A  0 

L 

C 

COMMON  BLOCKS 

PLT 

160 

C 

\a*:- 

RCUT INFS 

PLT 

170 

c 

EARTH 

MAIN  PRCO=AM, 

:rgls2. 

S AFR2  J, 

PUTGLT 

PUT 

180 

c 

AIR 

MAIN  PR^GnAV, 

~  =  GLS2, 

TEST2J, 

cAF*20 

PUT 

190 

c 

PUTCUT  AND  UNI 

our. 

PUT 

200 

c 

EIRE 

MAIN  PRCGFAw, 

"=CLS2, 

TF.ST2J, 

CAFW20 

PUT 

210 

c 

RATER 

MAIN  PRCGRAM, 

C  A  E n  20 , 

VCERIV, 

DOSFIV, 

PUT 

220 

c 

PUTQUT,  FN,  Cl 

E  ,  AND  PLCCK  DATA. 

PUT 

230 

c 

TIME 

MAIN  PR CO MAP# 

=  L  TCUT 

PUT 

290 

c 

SPACE 

MA  IN  P«CG 

5AFJ2G, 

YCt.R  IV, 

CDER IV, 

RUT 

250 

c 

ANC  FLTFUT 

PUT 

260 

c 

c 

MN 

MAIN  0CCGCAV* 

OLTCUT 

PUT 
—  PUT 

270 

280 

CCMMCN  /F/RTH/  CC£F,PNA“F,TINV, I FHA  S  E  ,  N  P F A SE , I KCUNT , I GC ,NSET S , 

1  ISTAT2.  MNVER,  INSTAT,  INVPH, INVSC  ,  L  1ST0 
CCM.MCN  /AIR/  X,  F,  CC ,  TITLE,  VO,  SICVC,  FD,  KI,  NC.  NV,  NX,  IVt, 
1  NP,  NO.  ISING,  ISTPP,  IL,  JOFLAG 
COMMON  /WATFR/2ERC,  GNc,  ThC,  THREE  ■  ECUS,  SIX,  A,  F, 

1  SCINV,  TRE F ,  STCCEE,  DIEO,  AOIE,  ECIE,  THETA,  YESNC,  T K , 

2  astar,  type,  nl,  nscale,  label,  bcc 


TO 


CC«MCN  /TIKE/  CATE 
COMMON  /SPACE/  QCCC,  QVQC ,  SC, 

CO“MCN  /MAN/  AA 

CP  RIME’* EEC 

LSTPFA  =  NPFASE(JI 

CO  1 2 JO  L  *1 , U5TPHA 

INOEX*l-  /  *(  I  PHASE  (L  ,  J  I  - 1 1 

C*C.SOJ 

IF  ( ISTATEIL, JI.KE.lt  GC  TC  ZOO 
CALL  DOER  IV  (TREF) 

CO  ICO  K*  1 , 7 
D=D*P( INOEX*K-ll*CCOC(Ki 
100  CONTINUE 
GO  TC  900 
ZOO  CJNTINUS 

IF  (MNVEPIL.JI.EC.OI  GC  TO  900 
CO  300  1*1, Z 
SGNi I l-SC  INV( 1 1 
300  CONTINUE 

IF  (INVSCIL.JI .EC. 01  GO  TO  500 
CO  90 J  LLL'l.Z 

SGNI  LLU  *  SGN  (LLL I  *STCCEF  I  INVSCIL.JI  I 
900  CONTINUE 
500  CONTINUE 

IF  I  INSTATE L, J) .EC.O I  GC  TO  700 
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PUT  290 
PUT  300 
PUT  310 
PUT  320 
PLT  330 
PUT  390 
PUT  350 
PUT  360 
PUT  370 
PUT  380 
PUT  390 
PUT  900 
PUT  910 
PUT  920 
PUT  930 
PUT  990 
PUT  950 
PUT  960 
PUT  970 
PUT  980 
PUT  990 
PUT  50 0 
PUT  510 
PUT  520 
PUT  530 
PUT  590 
PUT  550 
PUT  560 
PUT  570 
PUT  580 
PUT  590 
PUT  600 


InCcX2=  I  ♦  7*(  INVPHU,  J.ll-11 
CALL  DJERIV  (TREFI 
CJ  6CO  K»I,7 

C  =  C»P( INCfX2*K-I  I *OCCC ( K ) *SGN <21 
600  CONTINUE 
700  CONTINUE 

LSTINV=NINVF R(L,  Jl 
CJ  dCC  L S  1  =  1 , L  ST  I A  V 


PLT  610 
PLT  620 
PUT  63 C 
PLT  660 
PLT  650 
PLT  660 
PUT  670 
PLT  68C 


IP  (  (LST.E0.LST)NVI.ANU.(1NVSC<L,JI  .NE.OI  I  SGM  2  l=SGN(  2 1 /STC95F  (INPUT  690 


IVSCIL.JI)  PLT  70C 

CALL  CUdT.IV  (TINVIL,  J.LST1  I  PLT  71C 

CO  800  LLL* 1 i 2  PUT  720 

INRcX2si*7*IINVPh(L,  J,  LSTALLL- 1 1- 1 1  PLT  73C 

CO  SCO  K=  1 , 7  PUT  760 

C*C»P( IN0EX2*K-I |*CCCC(KI*SGN(LLLI  PLT  75C 

800  CONTINUE  PUT  76C 

600  CONTINUE  PUT  770 

00  1 30J  K ■ I »  7  PUT  780 

A( L ,  K I =  p  (  INDEX* K- II  PLT  7SC 

1000  CONTINUE  PUT  800 

A(L»6I’A(L»6I*D  PUT  810 

CO  UJO  K « 1 1  7  PUT  820 

AM  )NDEX«K-1)3«(L|KI  PUT  830 

1100  CONTINUE  PUT  860 

1200  CONTINUE  PUT  850 

LST»LSTPH»1  PUT  360 

CO  1300  K=l,7  PUT  870 

A< LbT.K I =C  .000  PUT  88C 

1300  CCNTINUL  PUT  890 

PO  I  500  L*I,LSTPHA  PUT  900 

CO  1 5 0 J  K«l, 7  PUT  910 

IF  I  (K.EC.7l.AN0.USTATt(L.JI.L0.-ll  I  GC  TC  1600  PUT  92C 

A(L3T,KI=A(LST,KI«A(L,KI*C0EFIL,J1  PLT  930 

GO  TC  I  5 C C  PUT  960 

1600  CPPIPE=CPPIPE»AIL,KI*COCF|L,JI  PUT  950 

1500  CONTINUE  POT  960 

w°ITF  (h,  ISCCI  OATE  PUT  97C 

CO  ltJO  L= I , LSTPHA  PUT  980 

ACCiE  FsClIE  FI  L  •  J I  PUT  990 

WRITE  (6,20001  PNAVF.  ( I  PHASE  ( L,J)1,ACCEF,(A(L,K|,K=1,7I  PUT1000 

1600  CONTINUE  PLTI010 

IF  (  (NPHASE1JI.SC. II. CR.IJGOIJl.GE. 711  GC  TC  1800  PLT1020 

IF  ( CPR  I PE.E0.2ERCI  GC  TO  1700  PUT1030 

NFITE  (  6,22001  PUT  1060 

WRITE  (  6,23001  ( A (L S T,K 1 ,  K* 1 , 7 1 ,CPRl  PE  PLT1050 

GO  TO  1 8CC  PUT  1060 

1700  CONTINUE  PUT  1070 

WRITE  (6,21001  PUT  1000 

WRITE  (6,23001  (A(LST,KI,K»1,7I  PLT1090 

1800  CONTINUE  PUT  1100 

RETURN  PUT  11 10 

1900  c JRMAT  (  1F0/1H0,  Afi,  1CH - /6H0PHA  S  E  , 8X  ,  6HC CFF  ,9X,  7H  1 A  » /  »G •  •  1PUT1120 

13X, 3H'0» ,  I5X,  3H'C ', I5X,3H»0*,  I5X,  3H'E“ , IEX, 3h' F • I  PUT  1130 

20  JO  Foo.NST  (  1H0,AR,6X,F6.3,  IX,  1P6U18.  7/lW  ,  1 9  X ,  IPO  I  8. 71  PUT  11  60 

2100  FORMAT  (  1F0/ 1HC, 2tX,  7H'  A' /'G*  ,  13X, 3t 'B» , IEX.OH'C' ,  15X,3H'0' , 15X, 3HPUT  1150 
l'F',  I5X, 3F'F' I  PUT  1 1 60 

2  200  FORMAT  ( I  FO/ IHO,  26X,  7F  »A  •  /  •  G»  ,9X  ,  1 2h  *0  »  /  •  Cl  IOM  • ,  I  OX,  3H»C  • ,  15X  ,3H  *PL  T 1 1  70 
ID* ,15X,3H»r» ,I5X,3H'F*I  PUT  1180 

2300  FORMAT  ( 20PORE ACT  ION  CCNSTANTS  ,  1P6C 18. 7/ (  IE  , 1SX, 1P2018. 71 1  PL  T 1 190 

ENC  PUT  1200— 


GO  TO  18CC 
1700  CONTINUE 

WRITE  (6,21001 

WRITE  (6,23001  ( A (LST ,KI  ,K>1 ,71 
1800  CONTINUE 
RETURN 

1900  c JRMAT  ( 1F0/1H0, AC, 1CH - 
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sue* cum e  prplct 

l^PU  !C1T  ICGICAL*1(I»),LCG1CAL*1CK) 

Pnr,;-?.,nN  NFCjlL-(51-  A8N0S<26),  X(l),  Y(ll 

S  ss:1. 

1NTi-GEK*2  VCR 

n4T„  FJX3/M1H0F  ,  F  ....  .  ,, 

C4T4  VCR/ 7AF00/  ’  ' 

OAT  A  KPLrTl/.FALSt./,KPLCT2/.FALSF./ 
t  CATA  <AaSC,KO»D,KPCTGL/3*.  FALSE./ 

JTl»LTPU2l'SCAL^,NhL,KSBH,,'VL,NSBV' 

K  PL  JT2-.  FALSE. 

NFWA8SCNFL  I 
NSF=  IASS ( 0S8HJ 
NV  =  iA3S(MiL) 

NSV= I A3S  (  KSHVI 
N  SCL  =NSC  A  LF ( 1 1 

ip.'moldor5''^-0'  G°  TQ  200 
10J  KPLuJIJIfs;!^  fLCTl  ARG*  "> 

RETURN 

200  KPL J  T- . TPIF. 

IP  (NV.Lfc.25l  GO  IC  AOO 
RRITE  (6,2001 
KPLJ  T». F A  I  St • 

300  o^uSn  (T!,,N0*  0F  VERnCAt  LINES  >25') 

AOO  CONTINUE 
NV*»=NV-I 
NVP=Nv»l 
NDh=NH*NSfc 
NOF'PjNOH*-  I 
N0V=NV«NSV 
NDVP-NQV* I 
M“G  =  (NOFP»NCVP| 
ip  (NCV.LE.120I  GC  TO  600 
KPLUT*. FALSE. 

WRITE  (6, EOOI 

5U°  PE  TURN  (T5',',ICTH  CF  GRAPH  >121*1 
600  CONTINUE 

IF  (NSCL.EQ.OI  GO  TO  70G 
F  SY= 10. *»NSCALE( 2  ) 

PS>=IO.«*NSCALP(AI 
IV«MINJ( 1A8S(NSCALE( 3  >> , 71*1 
I  X - M I N  J ( I ABS(NSCALE(5II ,9)f I 
GO  TC  800 
TOO  FSY*1. 

FSX»I. 

IY’A 

IX*A 


PRP 

10 

PRO 

2  C 

pop 

30 

PRO 

AO 

PRP 

50 

PRP 

60 

PPP 

7  C 

PRP 

80 

PRP 

90 

PRP 

100 

PRP 

no 

PRp 

120 

PRP 

130 

PRO 

1A0 

PRP 

150 

PRP 

160 

OTp 

1  7  C 

POD 

1HC 

po  P 

190 

POD 

2JC 

PRP 

210 

PRP  220 
PPP  230 
PPP  2A0 
PRP  250 

PPP  260 
PRP  270 
PRP  280 
PRP  290 
PPP  300 
PRP  310 
PRP  320 
PRP  330 
PRP  3A0 
PRP  350 
PRP  360 
PRP  370 
PRP  380 
PRP  390 
PRP  AOO 
PRP  A10 
PRP  A20 
PRP  A30 
PRP  AAO 
PRP  A50 
PRP  A60 
PRP  A  70 
PRP  A80 
PRP  A90 
PRP  500 
PRP  510 
PRP  520 
PRP  530 
PRP  5A0 
PRP  550 
P»P  560 
PRP  570 
PRP  580 
PPP  590 
PRP  600 
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800  FCA1(1J)*NCS(1Y) 

NA  =  M IN  J I I >  »  NSV I - 1 
N  S  =  N  A-y  I  NC(  N  A  *  1 20-N0  V ) 

NB  =  1  1-NS  *NA 
Il*iIE/10 
I2=NP-I  1*  10 
F0nj(6)=NCSl 11*1 ) 

F a n  3  <  7 ) =  NC  S  ( 12*1) 

F0R3(9)=NCS(NA*1  ) 

If  (NV.3T.0)  G2  TC  1000 
CO  9C0  J- 11  *  18 
900  FniJ(J|=ei 

cc  tc  hoc 

000  1 1=  J  V/ 10 

I  2=  9V- I  1  »  10 

FOR  J  ( 1 1 ) =  NCS ( 1 1 ♦ 1 1 

FC°  3 ( 12 )  =  NCS  f 1 2*1 1 

FCRj(l3)=NF 

I1=NSV/1CC 

T  3  =1 S  y-  I  1  MCO 

12=13/10 

13  =  12-12*10 

F0A3 <19)*N0S( 11*1 ) 

f  3R  3  ( 1 5  ) =NCS( 12*1) 

FOi(lo)  =  NCSI  13*1  I 
F  0  A  3  ( 1  7  )  *  F  F  1 
FCR3 (16! -FCF3(9) 
lluO  IF  (KPLJI1)  RETURN 
KPLJT1 =. TPUE. 

FNTK  Y  PLCT2I  I M AGE  ,  XPA X , XPl N, YMAX , YP IN  I 
KPL0T2=.TPUE. 

IF  ( KPLOX  1)  GO  TC  12C0 
NSCL=0 
NH  =  5 
N  SH= 10 
NV  =  1 0 
N  SV= 10 
CO  TC  200 
1203  CONTINUE 

IF  ( .N JT.KPLCTI  RETURN 
YPX* YMAX 

OH=CYPAX-YNIN)/FLCAT(NDHI 
CV  =  (  X'1AX->MIM/FLCATINCV) 

CO  1303  1=1, NVP 

1300  ABNOSI  I  )  =  (XPINfFLCATI  (I-ll  *NSV)*CV)  *F$X 
DO  1900  I = 1 i NI MG 
1900  IMAGC(1I»EL 

00  1600  1  *  1 «  NOHP 
12*1 «NOVP 
I1=I2-N0V 

KNF0S=Ma0CI-l,NSH).N6.0 
IF  I KNHOQ  |  GO  TO  1600 
00  1500  J • 1 1 i 12 
1500  IPAGEIJIoK 
1603  CONTINUE 

DO  1 80 J  J=I1,12,NSV 
IF  ( KNHOR  |  GO  TO  1700 
INAGFC ji =nc 
GO  TO  180C 
1700  1  PAGE  I J )  =  VC 


PRP  610 
POP  620 
PRO  630 
POP  690 
BAP  650 
PRP  660 
PAP  670 
PRP  680 
PAP  690 
PAP  700 
PRP  710 
PRP  720 
PAP  730 
POP  790 
PRP  750 
PAP  76C 
poo  770 
PAP  78C 

Poo  793 
FRO  800 
PAP  810 
PAP  820 
PA F  830 
Pop  890 
P»P  850 
PRP  860 
PAP  870 
PRP  B80 
PRP  890 
PRP  900 
PAP  910 
PAP  920 
PAP  930 
PAP  990 
PRP  950 
PAP  960 
PRP  970 
PRP  980 
PRP  990 
PRP 1000 
P  A  P  1 0 1  0 
PAP1020 
PRP 1030 
PRP1090 
PAP1050 
PRP 1060 
PAP  10  70 
PRP 1080 
PRP 1090 
PRP 1 100 
PAP  1 1 10 
PRP1120 
PRP1130 
PRP 1 1 90 
PRP 1 1 50 
PR  0 1 160 
PRP 1170 
PRP1180 
PRP 1190 
PR  P 12  JO 
PRP1210 
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18J0  CONTINUF 

xmi:ji*amin-cv/2. 

Yw  IN  1  =  YM 1  N-CH/2 • 

return 

c 

ENTRY  PLCT3(CH,X,Y,N3I 
IF  (KPLQT2)  GC  TC  21C0 
1903  wPITF  ( 6  ,  2C30) 

2000  C0R3AT  (  T5  ,  *  PLCT2  WLST  BE  CALLEC’I 
2100  COST IN'Ji 

IF  (  .NOT.KPLCTI  RETURN 
IF  (N3.GT.0I  GO  TC  2300 
KPLQT*. FALSE. 

RRITC  (  6,22001 

2200  FORM  AT  (T'.'PLOTS,  APG2  I  O'l 
PFTUPN 

23  jJ  PC  3  C 0 J  I  *  1 , N  3 

IF  (CV)  2300,2600,2500 
2600  CUM1=0 

GO  TO  2600 
25J3  CONTINUE 

0UP1  =  <  X(  I  1-XMlNl  I  /CV 
2600  IF  (CHI  2EOO, 2700, 2800 
2703  CUM2=0 

GC  TC  290C 
2800  CONTINUE 

CUM2  -  (  Y  (  I  l-YMINll/OH 
2000  CCNTINUL- 

IF  (  CUU.LT.J., CR.CUP2.LT. 0.1  GO  T  J  300C 
IF  (  DU11 . GC .NDVP.CP .CUM2.GE .NDHP I  CC  TC  2000 
NX  = l *  I  NT ( CUM  1 1 
NY* 1  *  I  NT (  rL'P2l 
J  =  ( (CHP-NYI*NOVP*AX 
1HAGF(JI*CF 
3 0 JO  CONTINUE 
RETURN 
C 

ENTRY  PICT6(NL,LAEELI 
ENTRY  FPLCT6(NL, LABEL! 

IF  (  .N  3T . KPLCT I  RETURN 
IF  (  .N3T.KPLCT2I  CC  TC  1900 
CO  3203  1*1, NDHP 

IF  (  I.EO.NCHP.ANC.KBCTGLI  GOTO  3200 
RL*BL 

IF  (  I.LE.NU  HL  =  LA  BEL(  1 1 
12*1 *NUVF 
1 1  =  1 2-NDV 

IF  (P3iJ(  I-UNSHI  .EQ.C.ANO..NOT.KCRDI  GO  TC  3100 
HRITf  (6.FCR2I  WL  ,  (  I  "AGcU  I  ,  J*  1 1, 1  2  ) 

GC  TC  320C 
3100  CONTINUE 

CRC  JC  =  (  Y“X-FL0AT( 1-1 l*DHI*FSY 
IF  (  I.EO.NDFP)  CRCNO=YMJN 

HRITF  (6, FORI)  WL .CRCNC, ( I  MAGE (  J I  ,  J* 1 1 , 1 2  ) 

3200  CONTINUE 

IF  (KABSCI  GO  TC  2300 
WRITF  (6.FCF3I  ( A  ENOS ( J I , J* 1 ,  NVP  I 
3300  RETURN 
C 

ENTRY  OMIKLSWI 
KAPSC=N0D(LSW,2I.EQ.l 


PPP1220 
PRP1220 
PRP1260 
PRP12  50 
PRP1260 
PRP1270 
PBP1280 
PR  P  1290 
PR°1300 
PR  P 13  10 
PRP1320 
PRP1330 
PR  P 1360 
PRP1350 
PR  P 1360 
PR  P 13  70 
RRP1380 
PRP  1290 
RR°16 JO 
PR  P 16  1 0 
PR  0  16  2  0 
PR  P  16  30 
PRP1660 
PRP 1650 
PRP166C 
PRP1670 
PRP 16  80 
PR  P 169  0 
PPP1530 
PPP1510 
PPP1520 
PR  P 15  30 
PRP 1 560 
PPP1550 
PRP 1560 
PRP 1570 
PRP 1580 
PRP 1590 
PRP160C 
PR  P 16 1 0 
PRP 1620 
PRP1630 
PRP 1660 
PPP1650 
PRP1660 
PR  p 670 
PRP1680 
PRR1690 
PRp  1700 
PPP1710 
PR  P 1 72  0 
PPP1730 
PRP 1760 
PRP 17  50 
PPP1760 
PRP1770 
PRP1780 
PRP1790 
PRP180C 
PRP1810 
PRP 1820 
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